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ABSTRACT 

We study various properties of the Z{3) interface which forms between 
the different ordered phases of pure SU{3) gauge theory above a critical 
temperature. We simulate the theory on a computational lattice of two space 
and one time dimensions, at various temperatures above this critical point. 
At high temperatures, we perform an accurate measurement of the interface 
tension, which shows good agreement with the prediction of perturbation 
theory. Near the critical temperature, we show that it is possible to observe 
the behaviour of the interface amidst the phase turbulence near its collapse, 
and we characterise its fluctuations in terms of various displacement moments. 
We show that the interface behaves like a fluctuating string, and model its 
behaviour in terms of an interacting scalar field theory. We also examine other 
properties of the interface, including its intrinsic and Debye electric screening 
widths, its formation and movement on the lattice, and its wetting. 
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Finally Cutangle spoke, very slowly and carefully. 

'I look at it like this, ' he said. 'Before I heard him talk, I was 

like everyone else. You know what I mean? I was confused 

and uncertain about all the little details of life. But now, ' he 

brightened up, 'while I'm still confused and uncertain it's 

on a much higher plane, d'you see, and at least I know I'm 

bewildered about the really fundamental and important facts 

of the universe. ' 

Treatle nodded. 'I hadn't looked at it like that, ' he said, 

'but you're absolutely right. He's really pushed back the 

boundaries of ignorance. There 's so much about the universe 

we don't know. ' 

They both savoured the strange warm glow of being much 

more ignorant than ordinary people, who were ignorant of 

only ordinary things. 

Terry Pratchett — "Equal Rites" (1987) 
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NESCME QUAEBAM MAGNA PARS SAPIENTAE EST 

Grotius 



For in much wisdom is much grief: 
and he that increaseth knowledge increaseth sorrow. 

Ecclesiastes 1.18 
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i neoretical r reliminaries 



1.1. Prologue: Phases in Physics 

:ii::::Ni|i:' ANY PHYSICAL systems possess assorted phases, undergoing transitions between 
;;;;l;\||:;iN: the different phases under certain conditions, for instance at a particular 
■■■■^■■^^■■^- ■■■■■■ temperature. Obvious examples include H2O, with phases as ice, water and 
steam and transitions at 0°C and 100°(7, and ferromagnets, with ferro- and paramagnetic 
phases. Phase transitions in systems such as these have been studied extensively in the 
past. Where different phases coexist, interfaces can form between them, and the interfaces 
have also been objects of previous study [1]. 

Quantum chromodynamics (QCD) is now generally accepted as the theory responsible 
for the "strong force" of particle physics. The gluon sector of this theory is described by the 
SU{3) gauge theory. This gauge theory itself exhibits several phases at finite temperatures, 
as will be detailed in the following sections. Interfaces can form between the different 
phases, and this thesis is concerned with the study of several properties of one particular 
type of interface, which may be of relevance to cosmology and to the heavy ion collision 
experiments of nuclear physicists at the Relativistic Heavy Ion Collider. 

This chapter will begin with some fundamentals of field theories at finite temperature, 
in particular SU{3) gauge theory. It will then continue with a discussion of the types of 
interface present in the theory, the mathematical framework for their description, and some 
of the predicted properties of one of the types. Subsequent chapters will proceed to consider 
this one particular type of interface, first at high temperatures, where it is most stable, and 
then at lower temperatures, where it is close to collapse. 

1 
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1.2. Introduction: Phases in SU(3) Theory 

T T has been known for some time that pure gauge theories have a non-confining high- 
temperature phase [2]. This was first shown non-perturbatively, in the strong-couphng 
limit of Hamiltonian lattice gauge theory [3]. In this situation, for instance, Abelian gauge 
theory can be shown to reduce to the "Villain model". This is a well studied statistical- 
mechanical system which is known to be disordered at high "Villain temperature", with 
a pair of opposite charges experiencing a confining linear potential; and ordered at low 
temperature, with a normal Coulomb potential. The Villain temperature is actually 
inversely proportional to the temperature, T, defined in gauge theory, which is used from 
now on. Therefore, any theory in a confining phase at T = must have a phase transition 
at some finite temperature, called the critical temperature {Tc)-, separating the confining 
and non-confining phases. These ideas are supported by numerical simulations [4]. 

A consideration of QCD at finite temperatures reveals that the pure gauge SU{2>) 
sector has the two expected phases. At low temperatures, there is the "disordered" phase 
where the colour charge of QCD is confined, and the vacuum is symmetric under the group 
Z(3), the "centre" of the SU{3) group (as discussed further in section 1.7). This is why 
no free quarks are seen in the relatively cool universe of today. At high temperatures, 
though, one finds a non-confining, Z(3)-breaking "ordered" phase, corresponding to the 
free quark-gluon "plasma" believed to exist in the hot early universe. In fact, owing to the 
Z{3) breaking, there exist three different high-temperature (high-T) phases, namely three 
degenerate vacua, corresponding to the three members of Z{3), the cube roots of unity. 
Fermionic matter breaks the vacuum degeneracy, albeit on a small scale, so we consider 
only the pure gauge theory in what follows. 

It is apparent that two different types of interface are possible in the theory. First, 
there is one between the ordered and disordered phases; this type is only stable at the 
critical temperature, since only then do the temperature and pressure match on both of 
its sides. Second, an interface can form between two of the ordered phases with different 
Z{3) vacua, and this is usually called an "order-order interface", or "Z(3) interface". All 
thermodynamic quantities are the same on both sides of a Z(3) interface, but they dip 
slightly across the object itself. Clearly, this second type of interface can only exist above 
the critical temperature, as only there does one find the ordered phases themselves. 

The aim of this thesis is to study some properties of the Z{3) interfaces, by means of 
Monte-Carlo simulations on a Euclidean lattice, in two very different temperature regimes: 
first, in the high-T limit, where analytic predictions have been made concerning properties 
such as the surface tension and Debye screening mass of the interface; and second, just 
above the critical temperature, where the interfaces are in the process of collapse. 
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1.3. Cosmological Implications 

"POSSIBLE cosmological significance has been claimed for these interfaces in the past, in 
the role of domain walls in the early universe, which could be responsible for matter 
formation[5] [6] . The suggested process is that of bubble nucleation, forming hadronic matter 
during the phase transition from the high-T quark-gluon plasma to the low-T hadronic 
phase. It has been believed for some time[7], and has been shown in simulations [8] [9], that 
the interfaces of S'f/(3) show something called "complete wetting". The meaning of this 
is that a Z{2>) interface between two ordered phases really consists of two order-disorder 
interfaces stuck back to back, leaving a very thin slice of disorder between two ordered 
domains. One possible scenario for the formation of hadronic matter in the early universe 
envisages many different domains of ordered phase at high temperature, separated by Z{'i) 
interfaces. As the temperature drops towards the critical value, the interfaces start to widen; 
and when the temperature hits the critical value, the constituent order-disorder interfaces 
split apart, spreading the slice of disordered hadronic phase outwards. 

However, it has more recently been emphasised[10][ll][12] that the order parameter 
distinguishing the SU{2>) phases is fundamentally a Euclidean object, with no counterpart 
in Minkowski space. Thus, whilst the Z{2>) symmetry exists, and hence Z{'i) domains and 
interfaces exist, in the Euclidean path integral, it is less certain that the interfaces can 
exist as physical objects in the universe. The claim has even been made[13] that only one 
true physical phase exists, even in Euclidean space, at high temperatures, and we shall 
address this when discussing our results later. In Minkowski space, the interfaces have been 
shown [11], on current understanding, to have unphysical thermodynamic quantities when 
a certain number of fermion families is considered in addition to the gauge fields, resulting 
from the unusual fact that the ordered phases occur at high temperatures. The free energy, 
F, is proportional to the fourth power of the temperature with a positive coefficient for many 
SU{N) systems with varying numbers of fermion families present: F = |7|T'*. Thus, the 
pressure, P = —\^\T^; the internal energy, E = — 3|7|T'*; the specific heat, C = — 12|7|r'^; 
and the entropy, S* = — 4|7|T'^. 

As pointed out in [12], no physical systems of this type can exist with positive 
temperature, but metastable states can exist with inverse population, i.e. negative 
temperature, flipping the sign of the speciflc heat and entropy as desired, and forming 
a type of quark-gluon "laser" . The problem may be assuaged by embedding the interface 
within a system whose positive entropy exceeds the negative entropy of the interface itself, 
and such a system has recently been postulated[14], with the interface contained within a 
larger volume of metastable vacuum in the early universe. However, it is noted in [12] that 
any statistical description of a system containing even a subsystem of negative entropy must 
fail, since no states are available to the subsystem. For instance, the standard electroweak 
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model well above the QCD phase transition has negative free energy density [5] [6], thanks 
to the contributions of the leptons, Higgs and weak gauge bosons; but if we consider the 
case where the quark- gluon subsystem has Fqg = \^\T^, then the weak coupling of the 
other particles cannot rectify the situation. Even setting this problem aside, the problem 
of interpretation of the interface remains, as the order parameter is non-local in Euclidean 
time, corresponding to an imaginary time-like gauge field in Minkowski space, and thus 
an imaginary chemical potential for the colour charge[15]. Adhuc sub judice lis est... In 
any event, by virtue of their contribution to the partition function and to expectation 
values calculated using the Euclidean path integral, the interfaces must be included in a 
non-perturbative analysis of the thermodynamics of QCD, if only (like the instanton) as a 
calculational device without physical form in our universe. 

In order to make analytical predictions about the interfaces, a mathematical description 
of them is necessary. For this reason, the next sections introduce some elements of quantum 
field theory at finite temperatures which will be of use in future chapters. 



1.4. The Partition Function and Thermodynamics 

/^NE now introduces a Schrodinger field operator at time t, (f){x,t), and its conjugate 
momentum operator, 7r(x, t)[16]. For convenience, one labels the eigenstates of the 
field operator at t = as \(j) >, satisfying 

<^(x, 0) I (/.>=(/.(x) !</.>, 

with eigenvalue 0(x). The amplitude for a transition from initial state {(pi > at t = ti to 
final state \(j)f > at t = tf is then given, as usual, by the following functional integral in 
(d — 1) space and one time {d or "(d — 1) -|- 1") dimensions: 

<cbf\e-'^^'f-'^^\cbi>= VirV^e -'^^ J ^ ^ ". (1.1) 

The integral over vr is without restriction, but the integral over (p is constrained by the 
requirements that (j){t = ti) = (pi and (j){t = tf) = (f>f. H is the Hamiltonian, and TC the 
Hamiltonian density. 

Now, the partition function of a statistical mechanical system in thermodynamic 
equilibrium is given by 

Z = Tt e-^^" = f d(t)< 0|e-'^^^|0 >, (1.2) 
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where /3t is the inverse temperature, T"-*^, and the integral sums over all states. From Z, 
all thermodynamic properties can be calculated: 

etc... (1.3) 



1 




dF 






dF 


\nZ 


P = 






S = 




(3t 




dV 


T 




dT 



Also, the expectation value of any operator is simply given by: 

<6>= ^Tr{de-^^^). (1.4) 

Zi 

We now generalise (1.1) and (1.2) to the logical conclusion of a path integral over gauge 
fields, before turning to consider how to distinguish between the different phases of a gauge 
theory at finite temperatures. 



1.5. Finite- Temperature Field Theory 

/'"^E can make (1.1) look more suggestive for Z by introducing an imaginary time coordinate 
T = it and putting ti = and itf = (3x, giving 



where previous functions of t are now functions of r, and (p = -J^. This enables one to write 

f^"^ dr J d''-^x{iTi<j>-n{Tr,4,)) 



Z = Vtt P0 e 

J -/periodic 

where the periodic boundary conditions on (p are that 0(x, 0) = (f>{x,(3T), arising from the 
trace in (1.2). This periodicity restricts the energy values associated with the fields to a set 
of discrete values. 

For a general renormalisable scalar field Hamiltonian, 

with C being the Lagrangian. One can perform the integration over vr to give 

-/periodic 

where the (irrelevant) normalisation constant has been discarded, and the Lagrangian is 
Euclidean. 

This technique can be extended straightforwardly to gauge fields. When applied to 
quantum electrodynamics (QED), for instance, one obtains the following functional integral 
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Z= I VA'^SiJ') det (^^ eJ 

Aperiodic \ ^^ / 



for the electromagnetic field alone, given for an arbitrary gauge specified by ^ = 0, where 
^ is a general function of the vector potential, A^^, and its derivatives: 

^^^ j;-d./d-^.(-iF,.F-)_ 

'periodic 

The electromagnetic field tensor F^"^ = d^A'^ — d'^A^, no colour indices being contained 
within A^ for QED, and a(x, t) specifies possible gauge transformations. 

For a general gauge group Q, with generators T"" and structure constants fate, the usual 
vector potential acquires a colour index, A^^, and thence the field strength becomes 

Fr = d>^A:-d''A>^-gU,A>^A-^^, 

introducing dimensionless coupling g. For a general gauge group SU{N), a runs from 1 to 
N"^ — 1. In a convenient shorthand, we shall later write the product over colour indices as 
A^^ = A^ • T. The infinitesimal gauge transformations are now given by aa{x,t), and there 
are A^^ — 1 gauge- fixing functions specified by J-"", so 

Z= [ V A^JCF") det (^) X^'^^''-'-'--Vl^'\ (1.5) 

•/periodic \ ^^ / 

To work in a physical gauge, one can simply choose 8{J-^^ to pick such a gauge within the 
generating functional above. However, it is often more convenient to work in an unphysical 
gauge such as the covariant gauge, where ghost fields ria,ria are introduced to give 

Aperiodic 



1.6. The Polyakov Line and Confinement 

T^'HE phases of a thermal gauge theory are characterised by the free energies of static 
configurations of quarks and antiquarks[17]. To calculate such free energies, one needs 
to consider creation and annihilation operators, ^J(xi,t) and ipa{xi,t), for static quarks 
with colour a at position x^ and time t, together with charge conjugates, tpl!^ &iid ip^, for 
antiquarks. All equal-time anticommutators of these fields vanish with the exception of 

{^a{^i,t),'>pl{^j,t)} = 6ij6ab (1.6) 

and its charge conjugate. 

For the general gauge group Q of the previous section, using the shorthand defined 
there for A^^, the quark fields obey the static time-evolution equation 



f -i^ -5^°(xi,r)j^(xi,r) = 0, 
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whence integration yields 



(1.7) 



Consider states | ^ > containing heavy quarks at xi , . . . , x^r and antiquaries at x'^ , . . . , x^_. 
The free energy of a collection of Nq quarks and Nq antiquarks is given by a sum over all 
such states, from (1.2) and (1.3): 

^-/3,F(x„...,x.,,x;,...,xW_) ^ ^ ^ ^|e-^-^|.& > . (1.8) 

Using the quark field operators, one can rewrite this expression in terms of new states |^ > 
without heavy quarks: 



-/3t-Fjv Ar_ 



Y.1^ Y. ^a,(xi,0)...Va.^(x^,,0)Vg^(x;,0)...^,V(^W'0) 



xe-'^-^<(xi,0)...^t (x^,,0)<^(x;,0)...VlL(x^^,0) 



(1.9) 



^ . 



Remembering that e~^'^^ will evolve an operator 0{t) in Euclidean time {0{t + Pt) = 
e'^^^0(r)e~'^^^), one now introduces pairs of exponentials e,-PTH^pTH ^^ ^j^g ^eH of 
each destruction operator in (1.9). Evolving these operators to time (3t and using the 
anticommutation properties to group operators into pairs at the same location, one obtains 

—I3tFn„n- 



*> 



E(^ E e-^-^Va,(xi,/3T)Vl,(xi,0)...Va.^(x^,,/3T)C(xiv,,0) 



{a,b} 



xV,^(x;,/3t)CK'0) • • • V'6._(x'iv^'/3T)C (^W'O) 
Using (1.6) and (1.7), one can finally write 

,-/5.F.,.^ = Tr(e-^-^L(x?) . . . L(x^ )Lt(xf) . . . Lt(x^J), 



(1.10) 
fl.ll) 



where the trace is taken over states of the pure gluon theory, and L(x) is the "Polyakov 
line", sometimes called the "Wilson line" (T denotes time ordering): 



L(x) = TrO(x), with 0(x) ~ Te'^/o ^ ^"^"("-") 



In particular, the normalisation in (1.8) for gauge group SU{N) is — nI+n- , and 



(1.12) 
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Recalling (1.4), one divides (1-11) by the vacuum, viz. the equivalent with A''^ = Nq = 0, 
to give 

^-^.AF.,.^ =< L(x?) . . . L(x^^^)Lt(xf) . . . Lt(xy >, 

where AF = -Fat, at- - Fqq. 

In particular, the expectation value of a single Polyakov line gives the self-energy of a 
single quark in the gluon medium [Fq = AFiq): 

Since these quarks are static and infinitely heavy, Fq is meaningless in itself. However, the 
correlation function of two Polyakov lines, F, gives the static quark-antiquark free energy 

{Fqq=AFuy. 

g-/3TF,-(x-y) ^^ L(x)Lt(y) >= r(y - x). 

Static quarks will be confined if this correlation function vanishes as |x — y| -^ oo, but 
cluster decomposition demands that 

r(x) =< L(0)Lt(x) > — > < L(0) >< Lt(x) >= I < L(0) > ^ 

x| — ^oo 

so that Fqq{x.) — > '^Fq, and so this condition is equivalent to < L(x) >= 0, showing that 

x| — ^oo 

the insertion of a single quark requires infinite energy, i.e. the theory is in a confining phase. 
In contrast, if < L(x) >^ then Fqq{x) will be finite and the colour charge deconfined. 
Thus, the Polyakov line gives an effective test for confinement. 



1.7. The Centre Symmetry and its Breaking 

'"puRNiNG again to the expression for the generating functional in (1.5), one is now led to 
consider the symmetries of the Euclidean action, 

drd^-'xF^^Fr, 

and the behaviour of the Polyakov line under the symmetry-related transformations. 

Returning to the shorthand notation of section 1.5, one sees that under a general 
non-Abelian gauge transformation 

Af" -^ UAf'U-^ + iUd''U-\ C/(x, t) e SU{N) (1.13) 

the Polyakov line will transform in the following manner: 

L(x) ^ Tr(t/(x,0)f7(x)C/t(x,/3T)). (1-14) 
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This means that L is invariant if U is periodic in time. However, Se is actually invariant 
under a larger group than merely the periodic gauge transformations; the only physical 
restraint is that the gauge potentials remain periodic in time. Let C/(x, t + [3t) = zU{x, t) 
in (1.13); since A^ must be single- valued, this tells one that z must be independent of x, 
satisfying z~'^A^z = A^ . This almost always implies that z is an element of C, the "centre" 
of the gauge group (the set of elements which commute with all members of the group): 
C = {z a G\zgz~^ = g for all g G G}. This global symmetry is an invariance under 
gauge transformations which are only periodic up to an element of the centre; it leaves 
local observables unchanged, and it remains after all gauge fixing. However, (1.14) reveals 
that the topologically non-trivial Polyakov line, which wraps around the periodic boundary 
condition in the time direction, is not invariant under these transformations, but rather is 
rotated by an element of the centre: 

L(x) ^zL(x). 

For the SU{N) groups, the centre is Z{N), whose elements are the A^th roots of unity, 
^2-!Tiq/N ^q _ 0, 1,...,A^— 1). This Z{N) symmetry may not be dynamically realised, 
however: if it is impossible to interpolate from one Z{N) vacuum to another through field 
configurations of finite energy and measure, the symmetry will be spontaneously broken. 
In the absence of such breaking, < L{x.) > must vanish, corresponding to Fq = oo. This 
will lead to an exponential decay of the correlation function at large distances: 

r(x) ~ e-'^i^'/^. 

From the expression for the free energy of a quark-antiquark pair given earlier, we obtain a 
linear potential with string tension cr{T): a confining phase. However, the symmetry could 
be broken at some finite temperature, with N possible broken-symmetry ground states 
(vacua) for the pure gauge system. 

How can the symmetry be broken at high temperatures? Consider the continuum 
generating functional again, for a moment, absorbing the dimensionless coupling through 

Z= f P^Me-^/o'^^^/^'-^-i^^.^. (1.15) 

Aperiodic 

If one rescales r -^ I^t'^-, ^^ -^ A^ /Pt^ this can be rewritten as 

/■p^^e-^r"-/'^'"^^^(^'^'+^'\ with A^{t = 1) = A'{t = 0), (1.16) 

where the electromagnetic tensor has been written explicitly in terms of its electric and 
magnetic components. As T ^ oo, the only non-zero contributions will come from field 
configurations with E = 0, implying that A*, the spatial components of the gauge potential. 
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must be static up to a gauge transformation. This, together with the boundary conditions on 
A^ ^ imphes that r2(x) must be an element of the centre group[18]. Thus, as -L(x) = Trr2(x), 
one sees that 

<L>=e2^^«/^/0, (^ = 0,l,...,iV-l), (1.17) 

drawing the conclusion that the global centre symmetry is broken at high temperatures, 
corresponding to a deconfined phase, with the N degenerate vacua being labelled by the A^ 
possible distinct expectation values of L, < L >= e^'^^^'^ Lq. This spontaneous breaking of 
the discrete global centre symmetry [<L>^ ^ Fq finite) will lead[19] to a correlator 

r(x)~ |<L>|2(l+/3re2e-'^l^l/|x|^), (1.18) 

giving a short-range static potential Fgg ~ — e^e"^'^'/ |x|P: the Debye screening of an electric 
plasma. 

The point to emphasise is that < L{x.) > is acting as an order parameter for the centre 
symmetry, distinguishing between broken and unbroken phases, and between the different 
broken phases. This has been seen perturbatively[20][21] and by simulation[17]. 

It should also be noted, from (1.16), that this high-T behaviour of a theory in d 
dimensions, with bare coupling g'^, is equivalent to that of a T = theory (where the 
electric field component of the action will also vanish) in (d — 1) dimensions (because of 
the static configurations) with effective coupling g^T. This effective coupling will become 
infinitely strong as T ^ oo, leading to strong-coupling behaviour. Thus, the Wilson loops 
mentioned in the next chapter show confining behaviour at all temperatures, and do not 
provide a valid test for confinement. 



1.8. The Background Field Method 

"\^7HEN we move beyond the classical level and seek to include quantum corrections, the 
manifest gauge invariance of a gauge field theory is usually broken by the gauge-fixing 
procedure. The "background field" [22] [23] quantises the field theory in the presence of a 
background field whose transformation properties are chosen to restore certain symmetries 
of the generating functional, and we shall rely on this method in the next chapter. 

The procedure is to write the gauge field from the classical action as A^ = A^ + a'^, 
where A^ is the background field and a^ the quantum field over which the functional 
integral is taken. The "background field gauge" is then used, breaking the gauge invariance 
of a^, but retaining that of A^. External sources are coupled only to a^, thus allowing 
quantum calculations without losing explicit gauge invariance. 
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First, conventional quantisation: adding a source J in the conventional fashion to the 
expression for the generating functional in (1.5), and writing x = (r, x) so that df^x = 
drdx = drd'^'^x, we have 

Z[J] = f VA^^det (^^^) el"'^ (/:(A)-^(^„(A))^-j»a^)_ (^_^g) 

./periodic \ ^^ / 

Here, 

1 



£(A) = --F;^Fr and J^a{A) = d^At^ (covariant gauge), 



4 
and ( "^"^ ^f'' ) is the derivative under the infinitesimal transformation 



dad 



M^(x) = -Df'aix), (1.20) 

where we define covariant derivatives for the standard and background fields as 

Dt" = df" - ig[A'' , ] and L>^ = 0^ - ^^[A^, ]. (1.21) 

The functional derivatives of Z[J] with respect to J are the disconnected Green's functions 
of the theory, with the connected functions generated by W[J] = lnZ[J]. We define the 
effective action by the Legendre transformation 

T[A] = W[J] - I d'^x J^Al, where ^^ = |^, 

and generate the one-particle-irreducible Green's functions by taking derivatives of the 
effective action with respect to A. 

Now, the background field: we define Z, W and T by analogy with Z, W and T: 

Z[J,Ab] = f Pa^det f^^^pA^) J ^"^ (£(A.+a)-^(^4a,A.))^-j>^)^ ^^ 22) 

Aperiodic \ ^^ / 



W\J,Ab\ = lnZ[J, As], r[5,^B] = W\J,Ab\ - I d'^x J^a^, a\ 



sw_ 



The single transformation of A in (1.20) now corresponds to connected transformations 
of Ab and a. Invariance of the action under (1.20) implies invariance under "quantum" 
transformations, where 5A^ = 0, and "background" transformations, where 5A^ = 
—D^a{x), with 5a^ being chosen to give (1.20) in each case. To maintain explicit 
gauge invariance of the action with respect to the background transformations, we choose 
J-a{a,AB) to transform covariantly under them, for example 

.Fa(a,ylB) = (Z)^a^),. (1.23) 

Now, ( Qd-^^) is the derivative under the quantum transformations, for which 5a^ = 

—D^a{x), so a consideration of J-'a{a + 6a, Ab) reveals that 

dJ-''^(a,AB) . „ ^^ , ^ 

Q^a = -{D%D,r,. (1.24) 
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By changing variables a ^ a — Ab in (1.22), we find the connection between F and T: 

f[a,AB] = T[a + AB\ => f[0,AB]=T[AB] 

The latter is invariant under the background transformation 6A'^ = —D'^a{x), so is a 
gauge-invariant functional of A'^ . 

We can now use (1-23) and (1.24) to rewrite (1.22), introducing complex Grassmann 
ghost fields r/ in a relation of the form 



so that 



Z[J,Ab] = [ Va^VrjVr] J""'' (^(^B+a)-i^(D-a,r+r^i-D-D,)r,-.r^a^J_ ^^ 25) 

Aperiodic 



Cnapter Z. 

rredictions lor Zfoj Interlaces 
in oU(o) dauge Ineory 



2.1. The Z(N) Interface 



:iii|i:i:ll ^ THE deconfined phase, where the Z{N) symmetry has been spontaneously 
lllllll broken, we have seen that the system can exist in any one of A^ degenerate vacua. 
mmmm Jt is possible to arrange boundary conditions of the system so that different 
parts of it exist in different vacua, i.e. different Z{N) phases. The existence of these 
distinct domains forces the appearance of "domain walls", or ^^Z{N) interfaces", where 
they meet [24] [21]. Within these interfaces, the gauge fields interpolate between the different 
vacua, as does the expectation value of the Polyakov line. 

The calculation of the interface tension is clearly an instanton problem: the field is 
interpolating from one Z{3) vacuum to another. An effective one-dimensional theory will 
describe the interface profile. The instanton is the solution of the classical equations of 
motion, with action ~ 1/5^ as in (1.15), and so we might naively expect the interface 
tension to have the same behaviour. To test this guess, we can construct an effective 
action from the classical piece, acting as a kinetic term, and a quantum piece, a potential 
term formed by integrating out fiuctuations to one-loop order. The Z{N) instanton is the 
stationary point of this total action. 

In a case which will be relevant for future chapters, consider a system of finite extent 
at a temperature T, as illustrated in fig. 2.1. The Euclidean time dimension runs from 
to (St = 1/r, and let us suppose there are two space dimensions, one {z) much longer than 
the other (x), and both much longer than the time dimension: L^'^ L^'^ fix- We assume 
that the system is in one Z{N) phase at the z = end of the longitudinal dimension and 
another at the z = L^ end, forcing the appearance of a Z{N) interface somewhere along the 
z direction. Where the Polyakov line was previously parametrised by the discrete values 

13 
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1^1 



1=0 




x=L 



z=0" 



z=L 



Fig. 2.1: A system of finite extent, with two space and one Euclidean time 
dimensions. 



q = 0,1, . . . ,N — 1 m (1.17), we shall now make q a continuous function of longitudinal 
position: q -^ qiz). 

As we have seen, the order parameter for the Z{N) symmetry is the Polyakov line, 
given in the fundamental representation by (1.12), where the gauge coupling is now written 
explicitly: 

L{z) = ^Trre''!o^^''^'^'\ 



N 



(2.1) 



We assume that the gauge fields obey strictly periodic boundary conditions. Since the 
time direction is finite, states with constant A^^ are inequivalent. Whilst we could work 
in a gauge with A^ set to zero at zero temperature, the periodic boundary conditions in 
the time direction make this impossible for T 7^ 0: were we to set A^ = by a gauge 
transformation of the functional integral, the A^ would violate the boundary conditions. 
Therefore, following the method outlined in [24] [21], we choose to make a global colour 
rotation so that A^ = {A^ , 0), with A^ a time- independent, diagonal N x N matrix: 



A'Hz) 



2TTTq{z) 

gN 



In-1 




-(iV-1) 





/I 






\ 


2TiTq{z) 










gN 




1 








\ 




-{N- 


-1)/ 



(2.2) 



In the presence of this A^ field, the Polyakov line is given by 

L{z) = e2-*«(^)/^ j^i - 1 [1 _ e-2-9(.); 



One can now understand the Z{N) structure of the vacuum: the A^ degenerate vacua occur 
for q{z) = 0,1, . . . ,N — 1, corresponding to L{z) = e27rig(z)/Af_ rpj^g ordinary perturbative 
vacuum has A*^ = q{z) = 0, with L{z) = 1. It has been shown[21] that the path taken by a 
gauge field when interpolating between different vacua has minimal action for this form of 
A^{z), as desired. 
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Returning to the generating functional in (1.25), without source term, and extracting 
the dimensionless couphng exphcitly again, the Euchdean action is given by 

The first term is the (classical) kinetic part, and the second and third terms form a potential, 
V. 

For the A^ = ^4*^(2;) field of (2.2), with transverse space dimensions of volume Vtrans, 
the kinetic part reduces to 



9^ J g^N J \dz 

Classically, there is no potential for A^, and all values of A'^ 7^ are degenerate. 
However, by computing one- loop quantum effects about the background field 

A'^ = (A°, 0), where A^ is as given in (2.2), (2.3) 

we shall see that the degeneracy is broken. We treat the background A^ field as constant 
in z, justifying this later. Expanding V to quadratic order in the quantum fluctuations a^, 
we can then integrate out both these gauge fields and the ghost fields, to give (to one-loop 
order): 



V = < iTrln 



(i?B)'V+(l-0^B^ 



Mn- 



Trln( -{Def 



from gauge fields from ghost fields 

The potential is independent of ^, so we use the Feynman gauge, ^ = 1, in our d dimensions, 
to reduce the previous expression to 

y=(^-l)Trln(-(Z)B)'). 

Now, each matrix function a{x) on which D^ acts can be written as a product of a 
scalar function of space-time with a vector in an (iV^ — l)-dimensional matrix space, where 
the basis matrices are the (A^^ — 1) generators of SU{N). To find Trln(— (D^)^), we should 
diagonalise Db in function space and then sum over its eigenvalues. A suitable basis is 
given by the plane waves ep{x) = e^^^*^ , with 4-momentum p^ = {po,p); the coordinates 
of a{x) are its Fourier components a{p), and the action of Db on these comes from (1.21): 

D^^a{x) = d'^aix) - iff[A^, a(x)] => D^a{p) = ip^a{p) - ig[A^^,a{p)]. 

Our choice oi A^ earlier in (2.3) eases the eigenvalue calculation. For the space components, 
the eigenvalues of — (D^)^ are simply (p*)^. For the time component, the periodic boundary 
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conditions force po = 2imQ/ (5t = 2TTnQT for integer hq. A suitable set of SU{N) generators 
enables the eigenvalues of — (D^)^ to be calculated rapidly; for instance, we can split the 
(N'^ — 1) generators into (N — 1) diagonal generators similar to (2.2), 

-(iV-2) 1, (i = l,...,iV-l), 



I 



N-i 



and [N — 1)! pairs of off-diagonal ladder generators [25], T^ •, with all elements zero except 
for one, whose value is common to the pair, at the jth row and ith column of T^ • and the 
ith row and jth column of T~- (i = 2, . . . ,N and j = 1, . . . ,i — 1). When these generators 
are commuted with A^, as given by (2.3), it should be apparent that only the {N — 1) 
pairs of generators Tjy • will give a non-zero commutation bracket (proportional to ±T^ ■ 
respectively), with the other (iV — 2)! = i(iV — 2)(A^ — 1) off-diagonal pairs and (A^ — 1) 
diagonal generators commuting with A^. Thus, the eigenvalues of — (Z)^)^ are seen to be: 

[27rT(no ±9)] : (A — 1) copies of each (charged generators); 

[27rTno] : (A — 1)^ copies (neutral generators). 

The part of the trace over space-time now simply gives the total space-time volume, Vtot = 
PTLzVtrans', howevcr, this is conventionally omitted from the "potential", so we shall ignore 
it until we insert V{q) into the action. Hence, putting V{q) = V{q) — V(0) so that V{0) = 0, 

V{q) = (A- l)(f - 1) T j;y"^^^|(A- l)ln((2vrTno)2+p2) 

+ ln((2^r(no - q)f + p') + ln((27rr(no + q)f + p')|. 

Note that for finite temperatures, we have replaced the integral over pq with a discrete sum: 

d^ ^ f d'^-'p 



d-1 • 



(2ttY ^^ J (2tt) 

Since uq and p run over all negative and positive values, the (no ± q) terms give the same 
contribution, and the subtraction of V(0) yields 

V{q) = 2(A - l)(f - 1) T j; I ^^^[ln((27rr(no + q))' + p^) - ln((27rrno)2 + p^)]. 

Now, we can use the following two integral formulae to make progress: 
dx = \n.b — Ina and / e""^ 



'0 ^ J -00 V 1 

The first enables us to rewrite V{q) as 

-1 /"OO 



v^,)^nN-m-i)TY.j^l r 



tp^r (27rTno)^t _ ^{2TTT{na + q)fU 
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before the second gives the result of Gaussian integration over p: 

V{q) = 2{N-l){l-l)TY, 



"°° dt 1 



t U-Kt¥/'^ 

no " " ^ ' 



t^{2nTnayt _ ^{27TT{na + q)fU 



A combination of the Poisson summation formula and a standard exponential integral leads 
to a useful relation: 



e-(""'+^"+^)(ix = ^e^b^-^"^^/^" and Yl ^(") = Y. \ e^^''''F{x)dx \ 

-OO ^__,^ T— -.^ W-oo J 



n^- — oo r— — OO 



— 7r{r / a)—27vibr 



n r 

Along with a change of variable to t = AnT'^t, this leads to 

°° /"OO if 

V{q)=2{N-l){^-l)T'^ V] / -^£-'^/2e-^'-Vt(i _ e2^^'-9) 

_, Jo t 



r= — OO 



A change of variables t ^ u = 1/t, followed by exponential integration over u, gives the 
final result: 

where, re-expressing the sum over all r as a sum over positive r, 



V{q) = m - l)(f - l)T'^^B,{q), (2.4) 



OO 

Bd{q) = J2-[l-cos{27rrq)]. 



r=l 



For SU{3) theory in 2 + 1 dimensions, V{q) is illustrated in fig. 2.2. The potential is 
clearly periodic: V{q) = V{q mod 1). The effective potential is minimised, as mentioned 
before, when q takes an integer value, i.e. it vanishes when < L >G Z(N), and < L > is 
indeed an order parameter for the A^ degenerate vacua. 

It should be noted in passing that this picture changes with the addition of fermions, 
as mentioned in section 1.3. This is because a gauge transformation of the form (1.13) 
causes the fermion wavefunction to transform simply as V' ^ U^l^. Thus, whereas the gauge 
fields are invariant under a Z{N) transformation, the fermion wavefunctions are generally 
not: t/j -^ e'^^^'^'^ip. This means that the boundary conditions, which should be periodic 
for the gauge fields and anti-periodic for the fermions, actually correspond to "twisted" 
conditions for the fermions: tpiPr) = e~'^'^^'^'^'ijj{0). Thus, the Z{N) symmetry is lost, 
and the structure of fig. 2.2 changes[10], with the central maximum and minima rising as 
fermion families are added but the outer maxima and g = minimum falling, the latter 
becoming the true vacuum. For A = 3 and four or more fermion families, the metastable 
minima at g = 1 and q = 2 actually correspond to positive free energy density of the form 
mentioned in section 1.3, and the peculiar thermodynamic properties discussed there set in. 
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Fig. 2.2: The effective potential, witfi temperature dependence factored out, for 
d = 3 and N = 3. Tfiis quantity is equal to 2Bz{q)/-K. The solid line is the 
continuum version given here; the dotted one is the lattice version for Nt =2, as 
calculated in section 2.3. 



We now have an expression for the effective action in g, remembering to integrate 
V{q{z)) over the total space-time volume, Vtot -^ PrVtrans J dz, and using (3t = l/T: 



■Seffiq) = VtransT / dz 



47r2(iV- 1) (dq 



+ 4r'^-2(iv-i)(g-i; 



For d = 3, we convert to the dimensionless variable z' = \jNTgz : 



^^^^BM\- 



^eff 



u 



T^-^ A-K^N -I) 



dz' 



dq 



+ 



BM 

47r3 



9 VN 

Treating the integrand as a standard Lagrangian with "time" z' , we know that the equivalent 

Hamiltonian (the energy, £ = (^)^ i^) ^^ ^ constant of the motion. In the simplest 

case, we set g = at one end of the system and g = 1 at the other, with ^ = at both. 
With these boundary conditions, the instanton has iS = 0, so 



L, 
L, 



9 VN 

9 



N 



dq^/Bsiq). 



47r3 



dz' J 7r3/2 
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The interface tension, a, is defined to be the effective action per unit area of the 
interface: 

^eff = PtLxOI. 

This gives a prediction for the surface tension of 



2^2.5 

= ao , 

g 



(2.5) 



as the integral is just a number, 1.018. Evaluating (2.4) for d = 3, we see that V{q)/T'^ 
(N — l)i?3(g)/7r, so, in a result that we shall use again later, 



N - 1 



«o = 4^y ^^ J dq^V{q)/T\ (2.6) 

For the SU{3) system we shall go on to consider, this means that 

ao = 8.33. (2.7) 

It only remains to justify our original assumption that q could be treated as independent 
of z. In fact, the instanton must have width of order 1 in terms of the dimensionless 
coordinate z' , as there is no length scale, and this corresponds to l/gyT in terms of z. For 
weak coupling, this is much larger than the natural length scale of the system, 1/T. Thus, 
the instanton field varies slowly enough in z that it can be treated as constant to leading 
order. 



2.2. Lattice Formalism 

A LTHOUGH some of the above properties can be examined in the continuum, many 
calculations are performed on the lattice [26] [27]. The advantage of this approach 
is that the short-distance cut-off enters the definition of the field theory in a natural 
way, allowing non-perturbative calculations to be made. In particular, lattice theories 
lend themselves to strong-coupling expansions and Monte-Carlo simulations, allowing the 
statistical evaluation of path integrals to a specified accuracy. In particular, the confining 
and non-confining behaviour of quarks can be seen in the lattice formalism, together with 
the transition between these two regimes. Using numerical simulations on the lattice, one 
can also approach the continuum limit and compare results with analytic predictions. 
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Time 



Longitudinal space 



Polyakov 
line 



Transverse 
space 



Fig. 2.3: The Euclidean lattice in two space and one time ("2+1") dimensions. 
Link variables, U, on the lattice are illustrated by bold arrows, as indicated by the 
example to the top left. A plaquette and a Polyakov line are shown; the Polyakov 
line wraps around the boundary condition in the time direction, but this is not 
shown. 



In the Euclidean formalism, one considers a hypercubic lattice whose sites are separated 
by lattice spacing a. This is the only quantity on the lattice which has dimension, and acts 
as an ultraviolet cut-off for the theory on the lattice. If we specify Nf links in the time 
direction, this corresponds to a temperature 



T = 1/(3t = 1/Nta. 



(2.8) 



Similarly, we have N^ links in the longitudinal direction and N^ in the transverse, so 



Lz = Nztt and Lx = N^a. 



For some gauge group Q, we can define a gauge field 



Ua = e*"-^-^"^ G g, 



(2.9) 



residing on the link leaving site n in direction ^u. The exponential of the field strength, 
e*^'"' , is now equivalent to the product of fields around a plaquette (four connected links 
forming the smallest possible square, as illustrated in fig. 2.3): 






(2.10) 
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The usual procedure is now to consider the Wilson action, which only depends on plaquette 
variables (A^ is the order of symmetry, as in section 1.5): 

Sw = — Y.(^- ^KeTrf/D- (2.11) 

•y n 

The conventional bare (dimensionless) coupling, g, is related to the lattice parameter (3 by 

P = 2iV/o4-V (2.12) 

in d dimensions. For our future purposes, d = 3 and A^ = 3, so that (3 = 6/ag'^ . 

This action is, naturally, invariant under local gauge transformations, [/^ -^ VnUj^V^,-, 
for Vn & G, as well as the global centre symmetry discussed in section 1.7. Using the gauge 
fields on the links, we can define "Wilson loops", gauge-invariant observables given by the 
trace of an ordered product around a closed curve C: 

W[C] = Tr n U^. (2.13) 

These topologically trivial quantities are not aware of the centre symmetry and Z{3) phases. 
Owing to the periodicity of the lattice in time, topologically non-trivial loops also 
exist, winding around the lattice one or more times; these cannot be shrunk to a point. 
The Polyakov line is one such, defined, as in (2.1) for the continuum, by considering a path 
parallel to the time axis at some spatial location x: 

Nt 

L(x)=Trn^xV (2.14) 

n=l 

An example of a Polyakov line is shown in fig. 2.3. Many of these non-trivial quantities are 
aware of the Z{3) phase structure, as we saw for the Polyakov line in section 1.7. 



2.3. The Lattice Surface Tension 

Q INCE we shall be examining the properties of an interface in SU{3) theory in 2 -|- 1 
dimensions, with Nt = 2, it behoves us to repeat the procedure giving the interface 
tension, this time on the lattice. This simply involves replacing continuum quantities in the 
calculation with their lattice equivalents, following [28]. 

If we consider the space dimensions to be essentially infinite, taking only the time 
dimension to be small, then the continuum momenta are restricted as follows on the lattice: 

2 . f Pia\ IT -K 

Pi^-sm — - , <pi<-; 

a \ 2 J a a 

po^ -sinl-YJ, Po = 2TmoT, 27r(no ± q)T, no = 0,1, ... ,Nt - 1. 
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The techniques used are largely the same as for the continuum calculation, but we leave the 
integration over the momenta until last, as the restriction is now —tt<p^ < tt (setting a=l 
for convenience) rather than — oo < p^ < oo. These integrals over momenta give modified 
Bessel functions: 

I^{t) = — f dk e'^^^^k+ink (2.15) 

The effective potential on the lattice is finally written (C.F. (2.4)) as: 

oo 

Viatiq) = m - l)(f - 1) ^(1 - cos{27rr q))F{rNt,d), (2.16) 

r=l 

where the newly introduced function F is defined in terms of the modified Bessel functions. 

In-. 

F{n,d)= / e-'''Ut)[Io{t)f-\ 

Jo ^ 

For d = 3 and Nt = 2, the values of F{2r, 3) are given in Appendix A. For these values, the 

lattice quantity, Viat{<i)Nf ^ is plotted in fig. 2.2 alongside the continuum version, V{q)/T'^. 

The lattice expression for the surface tension is given by (2.6): 



ao = 4^Y ^V^ / dq^Viat{q)Nf. 

For A^t = 2 and A^ = 3, this gives the prediction (C.F. (2.7)) that 

ao = 9.82. (2.17) 

To check that the correct continuum limit is recovered as Nt ^ cx), we need to find the 
behaviour of F{n, d) as n ^ oo. We use the method of steepest descent, which determines 
the asymptotic (s -^ oo) behaviour for an integral of the form 

l{s) = j dzg{z)e^f^^\ 

as follows: 

• Find the saddle point of /(-z), i.e. Zq such that the first derivative, /'(zq) = 0; 

• Expand around this point using z = zq + (5e*", and choose a so that ?R.e{f{z)} is a 
maximum, giving the path of steepest descent; 

• Substitute into the following formula, the approximation sought for large s: 



^/2^ff(zo)e^/(^")e^° 
"-^''^ \sf"{zo)\^/^ ■ 

From (2.15), we can write 

I^{nt) = ^ f dk e"(*^°^'=+^'=), Io{nt) = ^ / dk e"(*™^'=) 
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The forms of In{nt) and Io{nt) enable us to find saddle points at A; = cos i( v t^+i ^ g^j^^^ 
k = respectively, and the steepest descent in each case is for a = 0, giving 

In(nt) ss , — , lo(nt) « , 

^ ^ ^/2i^(l + l/t2)l/4 ' "^ ^ ,/2^t 

Substituting into our formula for F{n, d), 

Fin,d)= r ^e-'^-'Unt)[Io{nt)]'^-' 
Jo nt 

oo ^^ gn[t(^l + l/t2-l)-sinh-i(l/t)] • 



/o t (27rnt)'^/2 (1 + 1/^2)1/4 

This exponent has its stationary point at t = +cx), and, since the dominant contribution to 
the integral will occur around this value of t, we can replace the integrand by its approximate 
form as t ^ +oo. A final transformation of variable, t ^ u = 1/t, then reveals that 

poo d/2-l 

SO, by the definition of the gamma function. 



^, „ r(d/2) 

Fin, a) -^ ,,^ , as n ^ cx). 



as desired from a comparison of (2.4) with (2.16). 



2.4. The Debye Mass 

("p HE electric ("Debye") mass m, or inverse Debye screening length, governs gauge- 
invariant correlation functions of the time-like component of the gluon field {Aq) at 
large distances and high temperatures. The free energy of a quark-antiquark pair, over and 
above the sum of their separate free energies, vanishes as the quarks become infinitely far 
apart: 

Fg^(x) - 2F, = V{\^\,T) ^0 for T > T,. 



X — >oo 



From perturbation theory, we know that quarks are screened at a distance of the order of 
the inverse electric mass, so rather than being logarithmic in x, the interaction potential 
takes, in the manner of (1.18), the form 

F(|x|,T) ~ -C7e-2"'l''l forT>rc. 

x| — ^oo 

The factor of two arises because gauge invariance leads to an exchange of two gluons 
being the lowest-order contribution in perturbation theory. The electric mass is one of the 
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fundamental parameters of a gauge theory at finite temperature. One can perform a self- 
consistent and gauge-invariant calculation of its value, through the subtracted correlation 
function 

«L(0)Lt(x)> = <L(0)Lt(x) >- <L(0) ><Lt(x) >, 

and the Feynman diagrams contributing to this correlation function have been calculated [29] 
[30], with the dimensional regularisation mass scale set to m. We shall quote the expression 
found for m? later in this section, but it may be instructive beforehand to examine screening 
in a simple model: a plasma of bosons, each with mass nig and charge ib(7[31]. 

A positive test charge inserted into the plasma will attract (repel) negatively 
(positively) charged bosons, and the net induced charge density will result in an electrostatic 
potential (f){r). The energy of a given boson in the absence of the test charge, E = 
\ V^ + nvi^ will change in its presence by ^gcj), so we can write a consistent expression 
for this induced charge density using the bosonic density of states: 

pr^ - I (f~^p{ „ ,^^ ,. + ^ - ~^^, V (2.18) 

Far away from the test charge, 0(r -^ oo) -^ 0, and we can expand e^^^*^ ~ 1 + (3g(j). 
Changing variables so as to pull factors of /3t out of the resulting integral, we find that 
P ~ ~g^T'^4' for d = 4. The inclusion of the test charge in p gives exponential screening 
with this expression for the Debye mass: 

p — mr 

(/) , m'^r^g'^T^. (2.19) 

r 

For d = 3, we find that p ~ —g'^T(j)ln —r, and shows Gaussian rather than exponential 
screening. In fact, the effective one-loop potential for small q shows the following behaviour: 

Veff ~ q^ In q, 

leading to a divergence in the Debye mass. However, the full loop calculation mentioned 
earlier, treating the Debye mass self-consistently, reveals that the logarithmic divergence is 
halted at an infrared cut-off, resulting in the return of exponential screening. In our simple 
plasma model, this can be seen by an approximate evaluation of (2.18) for non-zero mass 
TUg, putting e^^-^ ps 1 + i3tE for rUg < E < /3^^: 

, . -,^m( r' E,E^\ttf + r E,E. '*" 



WtEY y^-. (eSrE _ 1)2 



-g'^T(j)ln I 1 to leading order. 
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m r^ g Tin 



m. 



The prediction of the full loop calculation [29] [30] mentioned previously is that 



Veff ^ q m , 



with a self-consistent solution for the Debye mass: 



m 



2 _ 9'NT 

4tt 



HT/g^) + O 



g^THHT/g')) 



(2.20) 



Clearly, this mass is a non-analytic function of coupling constant and temperature. This 
formula is only strictly applicable for T/g^ ^ 1. 
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Fig. 2.4: The predicted lowest-order behaviour of the Debye mass, m. 



The prediction for m is shown in fig. 2.4. In terms of the lattice parameter, (2.20) 
becomes 



m^ _ 3 



ln(/?/12) + O 



ln(ln(/?/12)) 



(2.21) 
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Since the one-loop V{q) ~ q'^ln{l/q) rather than q^ for small q, we are not able to 
extract an expression for the Debye mass from — o'°2 I _r,- 



2.5. The Deconfinement Transition 

"\A7^ can re-express the d-dimensional functional integral over gauge fields, seen in section 
1.6, in terms of the Polyakov line, L, by inserting a delta function conforming to (2.1), 
and so construct an effective action in terms of L [19]: 

e-Se//[L] ^ fvA^e-'''^^-^ll5(Lix)-Trre'I'''^"^''A. 

Integrating out the space-like gauge fields, A^, their correlation length, ^, sets the range of 
a set of coupling- and temperature-dependent kernels, 5*2, 5*4, . . ., in 

Seff [L] = fd^ V{L{i^)) + Jd^dy Lt (x)52 (x - y)L(y) 

-|- / dx dy dz dw L' (x)L' (y)S'4(x, y,z, w)L(z)L(7i;) -|- . . . . 

It can be shown[19] that ^ is finite at low and high temperatures, with space-like Wilson 
loops at both showing an area-law behaviour. Assuming that this will also hold at the phase 
transition itself, the kernels should be short-ranged. 

The deconfinement-confinement transition now appears as a (d— l)-dimensional order- 
disorder transition in statistical mechanics, the time dimension having gone. Integrating 
out the short-range interactions up to length ^ leaves a simple Landau-Ginzburg effective 
action: 



Seff= I dx{{d,Lf + V{L) 



Seff retains the invariance under the centre symmetry L -^ zL, so that the potential V{L) 
must do the same. However, we know little else about V{L), so we need to study it in 
mean- field theory (MFT), or by Monte Carlo simulation, to find out details of the order of 
the phase transition, as seen in [32] for the strong-coupling limit. 

In MFT we neglect fiuctuations, so the equilibrium state of the system is represented 
by the absolute minimum of V. If the minimum is at L = 0, we have the symmetric phase; 
at the transition, the transition is shown to be second order by the minimum moving away 
from zero smoothly, and first by it moving discontinuously. If the transition is continuous, 
renormalisation group (RG) theory tells us that an RG fixed point will govern the critical 
behaviour. Whole "universality" classes of theories have their critical points mapped onto 
the same fixed point, with the same exponents. This can be very powerful: if, in the space of 
(d — l)-dimensional theories with short-range interactions invariant under the same centre 
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symmetry, there is only a single fixed point, it reveals to us that a d-dimensional gauge 
theory will have the same critical behaviour as a d-dimensional spin system with the same 
symmetry. 

MFT predicts that systems with a global Z{3) symmetry should have first-order 
transitions. The symmetry allows a general effective potential depending on the cubic 
invariant 3fteL^, which we can then expand around the origin, choosing L to be real: 

V{L) = /(|L|2)+5(|L|2,3fieL3) = al"^ + bL^ + cL^ + . . . 

A sufficiently large value of a will give a unique minimum at L = 0, with a metastable 
minimum at L 7^ 0. The true minimum will jump to this second value when a is reduced 
towards zero, so the cubic term has produced a first-order transition. This argument has 
been confirmed for d > 4, but is invalidated for d = 3 by topological excitations. For 
our interests, SU{3) in d = 3 dimensions should have the same critical behaviour as the 
Z(3) Potts spin system, viz. a second-order transition described by the following critical 
exponents: 



Table 2.1: Critical Exponents 


Exponent 


Definition 


Regime 


Value 


Magnetisation 


<L(x) >~ |r-Tc|^^ 


T>Tc 


13m = 0.11 


Susceptibility 


x = Er(x)~|T-Te|-'^ 


T <Tc 


7 = 1.44 


Specific heat 


C = Td^ In Z/dT^ ~ |T - Tel"" 


T <Tc 


a = 0.33 


Correlation length 


^ = T/a r^ IT-Tcl-" 


T <Tc 


V = 0.83 


Critical correlation 


r(x) ~ |x|-(''-2+'') 


T = T^ 


r/ = 0.27 



The critical coupling, for 5*^^(3) theory in 2 + 1 dimensions, has previously been 
estimated [33] as 

/3e = 8.175 ±0.002 (2.22) 

in the infinite- volume limit of a lattice with Nt = 2, from measurements of (3c ~ 
8.04,8.11,8.12 and 8.14 for L^ = L.j, = 20,30,40 and 60 respectively. The coupling is 
related to the temperature, as previously, by /3 = VlTjg^. 
It has also been predicted[30] that 



and, leading from this, that 



/9. 



M 



(iV2 - 1)^2 



6 



0.23 for 5[/(3), 



levrTc 7r/3c 

but this is based on an approximate analytic calculation predating universality arguments, 

and so cannot be treated as reliable. 
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Cnapter o. 

Monte-Carlo lecnniques &: 
Results at Hign lemperatures 



3.1. Computer Simulation (2+1 dimensions) 

lii:ii|i:i:l:i' N ORDER to test the predictions of the previous chapter for the properties 
:|||||||||| of the Z{3) interface, it is possible to simulate pure gauge SU{3) theory on a 
■mmM^Mfr, computational lattice. In what follows, we have simulated the theory in 2 + 1 
dimensions, i.e. two space and one Euclidean time dimensions, as in fig. 2.3. One would 
prefer to work in the real world of 3 + 1 dimensions, of course, but it should be remembered 
that the interface profile is a one-dimensional object, interpolating between one Z{3) vacuum 
at z = and another at z = L^. All other spatial dimensions are therefore transverse, simply 
specifying the dimensions of the interface itself: in 2 + 1 dimensions, the interface is one- 
dimensional (a string), and in 3-|-l dimensions, it is two-dimensional (a sheet). An interface 
in a (2 -|- l)-dimensional system, as well as being worthy of study in its own right, may tell 
us much about properties in the higher dimension. The most obvious reason to work in 2 + 1 
dimensions, apart from the benefit to be had from studying the simpler, one-dimensional 
interface, is that by working in three Euclidean dimensions rather than four, we are able to 
considerably reduce the computational time needed for results of a desired accuracy. 

The SU{3) system with 2 + 1 dimensions provides an extremely useful testing ground, 
and is interesting in its own right. However, it differs from a similar system in 3 + 1 
dimensions in at least one fundamental respect: the nature of the phase transition at the 
critical temperature. In 2 + 1 dimensions, the transition is second order, meaning that 
the Polyakov line diverges smoothly from zero in the disordered phase, at and below the 
critical temperature, to its non-zero values in the different Z{3) phases above Tc- Thus, an 
interface present in this theory above the critical temperature would be expected to shrink 
continuously to nothing as T ^ Tc. However, the phase transition in the pure glue theory 
has been confirmed[34] to be first order in 3 + 1 dimensions, so that there is a discontinuity 

29 
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in the Polyakov line at the critical temperature, i.e. a jump from zero to a non-zero value as 
the temperature increases through Tc, followed by a further smooth change away from zero 
as the T -^ oo. Hence, while many inferences may be drawn from results in 2+1 dimensions, 
it is not clear that properties such as those describing the collapse of the interface at Tc will 
apply equally to 3 + 1 dimensions. It should be borne in mind, though, that a second-order 
phase transition in 3 + 1 dimensions is seen in an SU{2) system, and so those properties 
which do not carry directly over to an SU{2>) system in 3 + 1 dimensions may, instead, be 
relevant to an SU{2) system. 

In the computer simulation of pure gauge S'f/(3) theory, the fundamental variables 
are simply the link variables of (2.9). The variable U^^ is attached to the link leaving 
site n in the ^ direction, where ^ = 0, 1,2, corresponding to our one time and two space 
dimensions, and now U G S'f/(3). In an obvious extension of notation, U^'^ = U~t^~. We 
can construct physical observables from these link variables, as in (2.13), and then use 
Monte-Carlo procedures to compute a statistical average over all possible configurations of 
the system. Leaving the details of Monte-Carlo techniques used for the next section, we 
first deal with the creation of an interface on the lattice. 



t 



^^ ^^ ^' ^' ^' ^^ x^ ^^ .^^^ "V* 



Fig. 3.1: The Euclidean lattice in 2 + 1 dimensions: there are Nt links in the time 
dimension, r; the longitudinal space dimension, z, is considerably longer than the 
transverse dimension, x. The shaded plaquettes indicate the position of the twist, 
when present. 



The type of lattice used is shown in fig. 3.1, with two space and one (Euclidean) 
time dimensions. To study the phases of the system, we construct a Polyakov line at 
each point in the space dimensions (the x — z plane), using (2.14). We impose periodic 
boundary conditions in all directions, and set the temperature, specified through the lattice 
parameter, to be /3 ^ /3c, where /3c corresponds to the temperature of the deconfinement 
phase transition of (2.22). Left to its own devices, the entire system will settle into just 
one of the Z{?>) phases, this being the configuration of lowest energy. An example of such 
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a configuration, as may be verified from (2.11), is for all U^ = 1. We need some way to 
force the system into producing a Z{?>) interface, and this is achieved by judicious use of a 

"twist" [8]. 



1 1 1 z z 

11111 



Fig. 3.2: The effect of the twist: part of the t — z plane of fig. 3.1 for fixed t, x. 
Polakov lines to the left of the shaded (twisted) plaquette show a different Z{'i) 
phase from those to the right. 



The twist is introduced by modifying the action (2.11) so that every plaquette in the 
T — z plane at some particular value of r and z is to be premultiplied by an element, z~^, 
of the centre group when it appears in the action. This is illustrated in fig. 3.1, with the 
twisted plaquettes shaded. The action becomes 



S 



6 



w 



E 

shaded Ir 



1 _ is,ie(z-iTrf/p)) + J^ (1 - i^eTrC/p)) 



other 



where P represents a plaquette (C/p is the product of fields around a plaquette, as in (2.10)), 
and z is a member of the centre group: z G ^(3). Whereas the original action favours all 
[/ = 1 at high temperatures, so that f7p = 1, the twisted part now favours f7p = z to 
compensate for the pre-factor. This means that if we order the system so that all C/ = 1 
to the left of the twist, giving Polyakov lines in the Z(3) phase given by < L >= 1, then 
the twist will cause the link variables around the twisted plaquettes to behave as in fig. 3.2, 
with the right-hand time-like link given hj U = z. All equivalent time-like links to the right 
of the shaded plaquette at the same r will also take this value, so that so that C/p = 1 
again, making < L >= z to the right of the twist. No extra energy is associated with 
the twist itself; it merely acts as a change of variables, transforming one Z{2>) vacuum into 
another. However, the periodic boundary conditions force the appearance of a real physical 
Z(3) interface somewhere else in the z direction, where the Polyakov lines interpolate from 
z back to 1. Thus, we have forced the creation of two Z{'i) domains on the lattice, with one 
Z(3) interface between them. We shall place the twist at the boundary in our simulations, 
so that the left-hand end of the system is in one phase and the right-hand end in another. 
One might ask why the interface should form along the z direction rather than the r. 
In our simulation, the fact that Nt = 2 precludes the formation of an interface in the r 
direction; the dimension is simply not long enough for any interpolation to occur. Thus, 
thermodynamics favour an interface in the z direction, as long as Pt ■^ L^. 
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We thus set the initial conditions of our simulation to be that all link variables are 
unity except those in {r = , < z < L2/2 , x}, where they are e^^^'^ . The system is 
then allowed to equilibrate before measurements start, using the Monte-Carlo techniques 
described below. 



3.2. Monte-Carlo Techniques 

'"pHE object of Monte-Carlo methods [26] is to evaluate the expectation value of a quantity 
not by means of a path integral of high dimension, 

fVUOe-^ 
< O >= 

JVUe-s ' 

but by a statistical average over gauge configurations {U}i, generated according to the 
Boltmann distribution e~^: 



<0>=- Yl ^(i^}^)- (3.1) 



Each {U}i represents a set of link variables, U, on all links of the lattice. This technique is 
one of "importance sampling" , where the effort has shifted from the evaluation of the path 
integral, with many configurations contributing little, to the generation of {U}i with the 
Boltzmann weight. The statistical fluctuations in the mean value for O fall as l/\^imax as 
imax ^ 00. In practice, many tens of thousands of configurations are generally needed for 
a good estimate. 

In the "heat-bath" algorithm, one starts with a configuration {U}, and then changes 
one of the link variables to generate a new configuration {U'}. This new configuration 
replaces the old if the rules of the algorithm are satisfied. One "sweep" of the lattice is 
accomplished by applying this procedure to each link in turn. The object of the exercise 
is to generate, first of all, a configuration at thermal equilibrium: {C/}i, a member of 
the Boltzmann distribution. One then sweeps the lattice again until a new, statistically 
independent member of the Boltzmann distribution is generated: {U}2- Repeating this 
procedure many times, and calculating 0{{U}) for each, one obtains < O > through (3.1). 

The heat-bath algorithm works by touching a heat bath to each link variable in turn. 
Consider a particular link / with link variable U: all the links which interact with / have 
fixed U matrices, providing a background. On /, one chooses U' from the gauge group Q 
with a probability density proportional to the Boltzmann factor: 

dP{U') = e-^^^^'^dU'. (3.2) 
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This algorithm satisfies detailed balance, ensuring that any ensemble will be brought to 
equilibrium eventually. Moreover, since any other algorithm which varies one link at a time 
will simulate the heat bath after repeated application, this method must be the fastest 
to attain equilibrium in terms of the number of iterations. The challenge, therefore, is to 
ensure that the computational time necessary for a heat-bath iteration is on a par with that 
of alternative algorithms. 

A suitably fast algorithm for SU{2) was originally suggested by Creutz[35], and then 
an alternative by Kennedy and Pendleton [36]. One can write the Wilson action (2.11) in 
the form 

(3Sw{U) = (constant) - /?'3f?eTr(S/i), (3.3) 

with (3' = P/3, S a complex 2x2 matrix equal to the sum of products of link variables 
neighbouring /, and h the SU{2) matrix (link variable) associated with the link I to be 
updated. Thus, from (3.2), the problem reduces to the generation of an SU{2) matrix 
obeying the distribution 

g(/i)(d/i)oce'^'^^^(^'^)(d/i), 

where (dh) is the group-invariant Haar measure. We can parametrise all 2 x 2 matrices in 
terms of the 2x2 unit matrix, I, and the vector of Pauli spin matrices, a: 

a = oqI + ia- a, (a^ G M, a^ = Oq -|- a • a = 1). 

In fact, we define a = uh and (uq = ^~^5?eSo, u = ,^~^3f?eS), where a,u £ SU{2). In this 
notation, 3f?eTr(S/i) = ^Tra. Furthermore, the Haar measure takes a simple form, which 
can be rewritten in Polar coordinates: 

{da) = (5(1 — a'^)d'^a = ^y 1 — Oq "^(^ ~ \^ ~ ^o) sin 9 daodrdOdcj). 

Thus, we need to generate a^ with the (normalised) probability distribution 

Piao)dao = ^^\A^ e""Mao, 

TTli[a) V 

where a = 2/3'^, the factor of two coming from TrI, and Ii{a) is a modified Bessel function. 
We then generate Oi uniformly on a two-sphere of radius y^l — Og. 

The algorithm of Creutz generates such oq values with an exponential distribution, and 
then imposes the square-root factor by means of an accept/reject procedure. However, for 
our purposes, a is generally large, so that the normalisation pre-factor is small; this leads 
to a high rate of rejection. The Kennedy-Pendleton algorithm improves on this by working 



in terms of a new variable 5 = ^/l — Oq, so that 

P'{6)d6 oc -^1 - i(52 s^e-'^^^dd, {0<6<V2). 
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This algorithm is faster than that of Creutz for a > 1.7. It can be implemented as follows 
[36]: 

• Generate four uniformly distributed pseudo-random numbers in the unit interval: 
rt, JX , -tt , 11 , 

• Set X = -^lnR, X' = -^lnR', C = cos^{2ttR")- 

• Set 5 = X' + XC; 

• If i?'"^ > 1 — i(5, start again; otherwise, set Qq = 1 — 5. 

The extension of this algorithm to produce N x N unitary matrices with distribution 
(3.2) becomes more and more complicated as N increases, reducing the speed of the 
procedure, and favouring alternative algorithms to the heat bath, such as the Metropolis 
method[37]. However, Cabibbo and Marinari[38] proposed a faster algorithm for SU{N) 
based on updating SU{2) subgroups. Their method consists of first selecting a set 
{F : SU{2)k, k = 1,...,N - 1} of SU{2) subgroups of SU{N) such that there is no 
left ideal, i.e. only the whole SU{N) group is invariant under left multiplication by F. The 
choice of F that we use has elements of the form 



flfc 



ttfc 



1 



V 1/ 

where ak is an SU(2) matrix located at the kth and {k + l)th rows and columns. In each 
step of the iteration, the new link variable is obtained by premultiplying the previous value 
by one matrix from each of the A^ — 1 subgroups: 



U' = aN-iaN-2---aiU, Uk e SU{2)k, (k = 1, 



.N-l) 



We obtain the matrices a^ by using the Kennedy-Pendleton algorithm to choose the 
appropriate ak- It was shown in [37] that U' emerges with a Boltzmann distribution if 
U has one; thus, the procedure is ergodic. 

The Monte-Carlo study of lattice gauge theory often suffers from slow evolution of 
physical variables, owing to large correlations between configurations separated only by a 
few sweeps. Near a second-order phase transition, "critical slowing down" means that 
the autocorrelation time of the Monte-Carlo evolution diverges as the critical point is 
approached and physical correlation lengths increase. One way of circumventing this 
problem is to use "over-relaxation" to accelerate the heat-bath evolution. 

This use of over-relaxation was first suggested by Adler[39] for fields with quadratic 
coupling in the action. Brown &: Woch[40] and Creutz[41] extended the procedure to more 
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general theories. The idea is to pick a trial change for a link variable so that the system goes 
to an area of phase space as far as possible from that of the original value without increasing 
the energy by a large amount. This is done by finding roughly the locus of minimum energy 
for the variable, and then picking a trial value on the "opposite" side of this minimum. 
Consider a trial change of link variable h ^ h' . If a group element Hq, independent of h, 
minimises the system energy, then we wish h' to be on the "opposite side" of ho to h, for 
instance 

h' = /io/i~ /iQ- 

One can use a Metropolis algorithm to accept the new value if the configuration energy is 
unchanged or reduced by the change h ^ h' , and accept it with probability e"^*^ if it is 
increased. 

For an action of the form (3.3), a natural choice for /iq is the inverse of the group element 
obtained by projecting S onto the group SU{2). For SU{2), T, is always proportional 
to a group element, S. The energy of the new configuration will be the same as that 
of the old for SU{2), so the Metropolis algorithm will always accept this change; the 
algorithm is deterministic and microcanonical. To apply this to our SU{3) system, we 
again use the SU{2) subgroup procedure of Cabibbo and Marinari, over-relaxing each of 
the subgroups of our SU{3) matrices in turn. A previous study[42] suggested that for 
SU{3) gauge theory in 3 + 1 dimensions below the critical temperature, a mixture of seven 
Brown-Woch over-relaxation steps followed by one heat-bath step appeared optimal for 
reducing the autocorrelation time. Above the critical temperature, their results did not 
indicate any particular optimal ratio. For our simulations, therefore, there appears to be 
no particular optimal ratio to choose; we use seven Brown-Woch steps for every one heat 
bath, after many initial heat-bath sweeps to equilibrate the system. The link matrices are 
renormalised after each sweep, to ensure that rounding errors do not cause violation of their 
unitarity. Measurements of physical operators take place only every four sweeps, to reduce 
the correlations caused by adjacent configurations. 



3.3. The Wandering Interface 

"\A/^ can see in fig. 3.3 an interface on a 2 x 24 x 96 lattice {Nf x N^ x N^). The longitudinal 
direction goes across the page, and the transverse direction into it. The "height" in 
each case is determined by the value of the Polyakov line, and has nothing to do with the 
time dimension, which has been integrated over to form the Polyakov line. The twist has 
been placed at the boundary of the z direction, as mentioned before, so the left-hand end 
of the lattice is in one Z{3) phase with < L >~ e^'^*''^, and the right-hand end in another 
with < L >~ 1. The temperature is fairly high: /? = 50 ^ /3c. Thus, the interface is 
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Fig. 3.3: The real (top) and imaginary (bottom) parts of Polyakov lines on a 
2 X 24 X 96 lattice with twist. The interface can clearly be seen between the 
< L >ss e^'^*" domain on the left and the < L >ss 1 domain on the right. The 
twist is at the boundary between z = and z = Lz, and should not be confused 
with the position of the interface. The temperature is fairly high (/? = 50), resulting 
in a high suppression of the fluctuations in phase. 
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Fig. 3.4: The real (top) and imaginary (bottom) parts of Polyakov lines on a 
2 X 24 X 96 lattice without twist, also with (3 — 50. 
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fairly well-defined and stable, since its free energy is correspondingly high, suppressing the 
formation of fluctuations in phase. In fig. 3.4, the same system is seen without a twist. 

In this picture, the interface is roughly in the centre of the lattice. However, the 
interface is actually translation invariant; there is nothing in the action to specify that it 
should be located at any particular value of z. We shall always start the interface in the 
centre of the lattice, as mentioned in section 2.1, but for this chapter, we shall allow it to 
wander freely along the lattice. We can keep track of the interface most easily by examining 
a transverse- averaged profile of the Polyakov lines. 
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Fig. 3.5: Profiles of the Polyakov lines on a 2 x 16 x 48 lattice for /? = 50, averaged 
across the transverse direction. The real (solid) and imaginary (dotted) parts of the 
Polyakov lines are shown. The first picture shows the initial configuration before 
Monte-Carlo evolution begins; the second, the configuration after 20k heat-bath 
sweeps; subsequent ones, 40k sweeps apart. 



We can see the interface walking along the lattice in fig. 3.5. Placed at the centre of 
the lattice before Monte-Carlo evolution begins (first picture), the interface has taken up its 
shape of least energy after twenty thousand ( "20k" ) equilibration sweeps (second picture) . 
It has also moved along the lattice, and continues to do so in subsequent pictures, which 
are 40k sweeps apart. In particular, the interface often approaches the twist. In the 200k 
picture, the interface is passing through the twist, and in the 240k picture, the system has 
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tunnelled into a different ^(3) vacuum: now, < L >= e^^*/'^ on the left and < L >= e~^^*/'^ 
on the right. After 260k sweeps, the system has tunnelled again. 
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Fig. 3.6: The movement of the interface with computer 'time', for a 2 x 24 x 96 
lattice and various /3. 



The movement of the centre of the interface as proceeds is shown in fig. 3.6, a plot of 
the square of the position of the interface centre against the iteration number. To make 
this particular measurement, we allow the interface to move each sweep, but then slide the 
whole gauge configuration along the lattice so that the interface is put back in the middle 
of the z axis; essentially, we have put the interface on a treadmill, so that it is allowed 
to move, without really going anywhere. This technique will be used throughout the next 
chapter, and will be described in more detail there. It is used here simply to cure problems 
that would arise in identifying the interface as it passed through the twist. 

In a random walk with step size a, the lattice spacing in our case, the probability that 
at time t the object will be at position z is 



implying that 

< z^{t) >= ca^t, 
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where the speed of movement is given by c. The interface, from fig. 3.6, clearly obeys this 
linear relation for all the values of fi measured, demonstrating that it is indeed executing 
a random walk, in accordance with translation invariance. As the temperature decreases, 
the speed of the walk is seen to increase slightly. This is also expected, as the energy 
penalty suppressing fluctuations in phase will decrease as the temperature drops, allowing 
the interface to move more swiftly by combining with, and budding off, fluctuations nearby. 
We shall look again at the random walk of the interface in the next chapter, at much lower 
temperatures. 



3.4. The Surface Tension 



T~iHE free energy of a system is -F = —TlnZ, as in (1.3), so 



d_fF\ 1 



a3^5.;-5<S(C/)>. 

This gives us the behaviour of the free energy in terms of the average plaquette action. We 
can perform simulations for the same lattice with and without twisted plaquettes, i.e. with 
and without a Z{3) interface, and the difference in free energy between the two gives us the 
free energy of the interface itself, from which we can measure its surface tension: 

^interface ^twisted ^untwisted CV/1, 

where a is the surface tension, and A is the 'area' of the interface (in our case, the transverse 
size of the lattice). Thus, 

Therefore, recalling (2.5) and substituting A = L^a, 

AS d faoT'^-^L^a\ L^ 






(3 d(3\ g T J 2,/6p Nt^ 



ao, (3.4) 



where we have used (2.8) and (2.12) to substitute for T and g. 

The prediction (2.17) that ckq = 9.82 for a system with Nt = 2 was made in the infinite- 
temperature limit. Hence, to test this prediction it is necessary to estimate Oq in this limit. 
The best way to do this is to carry out a survey of ao for various values of /3, and then to 
extrapolate the results to the infinite limit. 

To estimate oq for a particular /3, we need to measure the difference in total plaquette 
action between two simulations which are identical except for the presence or absence 
of a twist. Several tens of thousands of heat-bath sweeps were used to ensure complete 
equilibration of the lattice before measurements began. "Main sweep" refers to a sweep 
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after equilibration is complete. As mentioned in section 3.2, the ratio of over-relaxed to 
heat-bath sweeps was seven to one, and measurements of physical quantities were taken 
every four sweeps. A DEC 2100 A500MP machine was used for most of the computation, 
and the simulation details are listed in the following table. The CPU time listed was needed 
for runs with or without twist, and therefore needs to be doubled to give the total time 
required for an estimate of Uq. 



Table 3.1: Simulation Details 


/3 


Lattice Size 


Heat-Bath + Main Sweeps (OOO's) 


CPU Time (hrs) 


16 


2 X 12 X 30 


15 + 400 


40 


28.125 


2 X 12 X 36 


20 + 400 


60 


50 


2 X 16 X 48 


20 + 400 


130 


112.5 


2 X 24 X 72 


40 + 400 


400 


200 


2 X 32 X 96 


40 + 400 


800 



For each of these runs, with and without twist, we now list the results for the average 
action for a single plaquette, < S >. This has to be multiplied by the total number of 
plaquettes (3 x iV^ x A^^ x N^) to give the total action of the configuration. The difference 
between the two values is then used in (3.4) to obtain the measurements of ao ii^ the last 
column. 



Table 3.2: Simulation Results 


/3 


^ 'J -^twisted 


^ ^ ^untwisted 


ao 


16 


0.17510(1) 


0.17379(1) 


13.1(1) 


28.125 


0.097575(6) 


0.096832(6) 


11.79(13) 


50 


0.054310(3) 


0.053916(3) 


11.12(12) 


112.5 


0.0239740(7) 


0.0238074(7) 


10.58(6) 


200 


0.0134549(5) 


0.0133644(5) 


10.21(8) 



One would not expect the finite size of the lattices used to have a great effect on the 
untwisted results, as the whole lattice will be in one phase, but it is less clear how great the 
effect would be in the presence of the twist. The shape of the interface may not flatten out 
completely at opposite ends of the longitudinal dimension if the lattice is too short, skewing 
our measurement of the total interface free energy. The lattice sizes above were chosen to 
be large enough for the finite-size effects to be negligible, by studying results for /? = 50 on 
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various sizes of lattice and scaling the spatial dimensions with the interface width, which 
is proportional to y^, as we saw at the end of section 2.1, and which we choose as our 
measure of physical correlation lengths. The only exception to this was the measurement 
for (3 = 16, for which the lattice could not be safely made any smaller. The results are 
given in the following table, each from 2+20k sweeps. Similar runs without twist gave a 
result of < S > untwisted^ 0.05390(1), and this value was assumed for each lattice in the 
calculation of oq: 



Table 3.3: Finite-Size Survey: /3 = 50 


Lattice Size 


< >^ ^twisted 


ao 


2 X 16 X 32 


0.05449(1) 


11.1(3) 


2 X 16 X 40 


0.05438(1) 


11.3(3) 


2 X 16 X 48 


0.05432(1) 


11.9(3) 


2 X 16 X 64 


0.05422(1) 


12.0(3) 


2 X 24 X 48 


0.05432(1) 


11.9(3) 


2 X 24 X 60 


0.05423(1) 


11.6(3) 


2 X 24 X 72 


0.05418(1) 


11.9(3) 



This suggested, given the errors, that a lattice size of 2 x 16 x 48 was adequate for 
/3 = 50. 

Now, as mentioned above, the prediction of (2.17) is only valid in the limit (3 -^ oo, so 
the question arises of how to extrapolate our results to test this prediction. This is really 
a question of calculating infra-red divergences. In 3 + 1 dimensions, one would expect the 
one-loop correction to (2.17) to go like 1/(3 in this limit, since it was seen in (2.19) that 
the Debye mass, m ^ T. However, in 2 + 1 dimensions, it is conceivable that a dependence 
on ln/3 would need to be taken into account, as this appears in the infra-red correction for 
the Debye mass in this case (2.20). For this reason, it is difficult to predict the form of the 
/3-dependence. In any case, it would be quite possible for the logarithmic dependence to 
masquerade as a power law over a restricted range of (3. 

In the absence of a definite prediction for the form of the finite-/? correction, we perform 
the extrapolation using a best fit to a power of 1//3. A computer-generated fit to the 
data, using the IDL package and its "COMFIT" function and leaving the power as a free 
parameter, suggests a power of 0.85 with x^ = 2.39. This gives a (3 = oo limit of ao = 
9.91(1), as extrapolated in fig. 3.8. IDL fits for various imposed powers give a lowest x^ of 
2.34 for a power of 0.80, with oq = 9.83(1). However, x^ < 2.6 for powers between 0.75 
and 0.90. Therefore, in the absence of a definite prediction for the power of 1//3, we allow 
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0.15 



Fig. 3.7: Extrapolations of ap are shown as /3 — > oo for various powers of [3. For 
each power, the same five data points are plotted with a particular symbol, and 
the accompanying line of extrapolation is marked with a corresponding label for 
the X-axis. 



for a reasonable error in the fit, and use 0.85(10), the range for which extrapolations are 
illustrated in fig. 3.7. Thus, our infinite-temperature extrapolation for the interface tension 
is 

ao = 9.91(1)(14), 

in very good agreement with the instanton calculation result of (2.17). The first error is 
associated with the fit for l//?'^'^^, and the second with the uncertainty in that value for 
the power. 

It has been argued[13] that the different Z{'i) vacua distinguished by different values 
of the Polyakov line actually correspond to one and the same physical state. This claim 
comes from considering the role of infra-red divergences in the calculation of the surface 
tension, and, if true, would imply that the surface tension is actually zero (since no physical 
interface can exist). The agreement of our results for oq with perturbation theory tends 
to disprove this argument. It is, of course, conceivable that our measurements for oq 
would dramatically fall towards zero if we went to higher and higher /3, meaning that our 
extrapolation is simply an artefact for the range of (3 that we have considered. However, 
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Fig. 3.8: Extrapolation of ao for preferred power of 1/(3 as (3 ^ oo. 

our highest value is /3 ~ 25/3c, a temperature well above the critical point, and it seems 
unlikely that any dramatically different behaviour will set in at a temperature higher than 
this. A sceptic could also argue that a lattice with Nt = 2 could show behaviour quite 
distinct from that in the continuum; but quite apart from the agreement with the lattice 
prediction of Oq for this value of Nt, a similar survey carried out for SU{2) gauge theory 
for a range of A^t < 5 has found similar agreement with perturbation theory [43]. 



3.5. The Debye Mass 

T^iO measure the Debye mass on the computer, we need to study the correlations of 
fluctuations in the gauge fields. The behaviour of these correlations will give us 
information about the Debye screening of sources, as in (1.18). We need to find a "good" 
operator, viz. one whose correlator rapidly converges to a decaying exponential of only the 
lowest mass state. This will allow us to obtain, in a straightforward manner, an estimate for 
this lowest mass. A simple plaquette correlator does not fit our requirements, as it depends 
on flux loops between the plaquettes being considered, and thus has to be considered over a 
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long range of separation before a single exponential dominates. However, a good operator 
turns out to be given by a transverse average of Polakov lines, as shown in fig. 3.5. 
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Fig. 3.9: Estimation of the Debye mass from iJi{z) for untwisted data aX [3 — 200. 



We measure the correlation function of this operator along the untwisted lattice, as 
it happens to be easier to obtain accurate results from this system than from the twisted 
lattice, where the Z{2>) interface acts as a dominating source term and the Debye mass can 
be estimated from the shape of its tail. Owing to the periodic boundary conditions in the 
space directions, we expect a fluctuation at one value of z to influence another point further 
along the lattice by correlations going in the +z and —z direction. This gives correlations, 
at long distances, proportional to g"^™"^ (as in section 2.4) and e~2m.(Lj-z)^ Therefore, we 
estimate 1m by a cosh fit, an improvement on the following function which fits simply to a 
decaying exponential: 

■<Lt(z-l)L(0) >- <L>2~ 



/x(z) 



In 



<Lt(z)L(0) > - <L>2 

where each < L{^z~) > represents an average of L{z) over all configurations (sweeps), and 
< L > is a further average over all z. From the definition of the Debye mass, we expect 



/i(z 



oo 



2m. 
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In fig. 3.9, a typical set of data is shown for ii{z) at high (3. For each similar set of data, 
we extracted a value of ii{z) corresponding to the level at which the function flattened out 
before errors began to increase wildly. These values are given in the following table for the 
simulations listed in table 3.1, in units of the inverse lattice spacing. 



Table 3.4: Debye (Electric) Mass (units of a~^) 


a 


Twice Debye mass, /i = 2m 


16 


0.80(1) 


28.125 


0.690(5) 


50 


0.555(5) 


112.5 


0.410(5) 


200 


0.332(4) 
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Fig. 3.10: Results for the Debye mass. The prediction of (2.21) is plotted 
to leading-order only (dashes, C.F. fig. 2.4); with an additional correction of 
ln(ln(/3/12)) (dots); and with an alternative additional correction of 
0.14 ln(ln(/3/12)) +0.94 (dot-dashes). 
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The values are plotted in fig. 3.10, together with the theoretical prediction from fig. 2.4 
(dashed line), which is expected to be valid in the high-temperature continuum limit. The 
results show some qualitative agreement with the leading-order prediction of (2.20) and 
(2.21), as plotted previously in fig. 2.4. However, although the shape of the rise of the mass 
with temperature looks correct, the results differ from the prediction by a factor of two or 
so. The reason for this becomes clearer when we consider the size of the O ln(ln(/3/12)) 

correction compared to the leading-order ^ ln(/3/12) that we plot: for our range 16 < /3 < 
200, the leading-order term is between 0.1 and 0.7, but the correction is 0(— 1.2) to O(l.O)! 
Thus, we shall find it very difficult to draw any firm conclusions from the data to support 
or deny the prediction. 

To show the effect of the correction, the dotted curve in fig. 3.10 is ^ln(/3/12) -|- 
ln(ln(/3/12)). Our data certainly look more reasonable when this correction is included, but, 
of course, this estimate is only of the order of the correction, not of an exact expression. A 
free fit to a correction of the form "c^ ln(ln(/3/12)) -|- C2" gives the dot-dashed line in the 
fig. 3.10, with ci ~ 0.14 and C2 ~ 0.94, which shows extremely good agreement with our 
data. Although we have two free parameters and only five data points, the fit is strikingly 
good, and this gives us some reason to hope that our data is consistent with the theoretical 
prediction. Of course, no firm conclusions can be drawn without a more explicit expression 



for the O 



ln(ln(/3/12)) 



correction. However, the ratio of coefficients of the free fit does 



at least confirm that most of the temperature dependence is contained within the leading 



^ ln(/3/12) term. 
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4.1. Interface Behaviour Approaching Criticality 

lii:ii|i:i:l:i' N THE previous chapter, we examined some properties of the Z{3) interface at 
:|||||i|||| high temperatures, (3 ^ 15c ^ 8.175. It was seen that the presence of a twist in 
l:p:::M::::::::N::N::::;::::; Qur simulations is sufficient to put the left- and right-hand longitudinal ends of 
the lattice in different Z{'i) phases, with one physical Z{2>) interface separating the two. At 
these high temperatures, the phases are very well defined, with only small fluctuations seen 
in the value of the order parameter (the Polyakov line) across their extent. The interface 
is similarly well defined, with a fairly rigid shape. Both of these aspects are due to the 
suppression of fiuctuations in phase resulting from the increased energy penalty at high 
temperatures. 

Consider, however, the phase structure at the critical temperature itself, /?c, where the 
deconfinement transition occurs. At and below this temperature, the three ordered, Z{3)- 
broken phases cease to exist, being replaced by a single disordered phase with symmetry 
restored. Obviously, therefore, no Z{3) interface can exist in this regime either. The order- 
disorder, confinement-deconfinement transition itself has been studied in the past. In [8], 
for instance, the surface tensions of order-order and order-disorder interfaces were compared 
in 3 -|- 1 dimensions near the phase transition, and various thermodynamic properties were 
examined. Also, the prediction of complete wetting, mentioned in the introductory chapter, 
was tested: viz. that an order-order {Z{3)) interface consists of two order-disorder ones 
stuck back to back, so that aoo ~ 2aod- A question that has not been addressed sufficiently 
in the past, though, is that of the behaviour of the interface: how does this change from 
the rigid structure we see at high temperatures, as we reduce the temperature towards the 
critical value? This chapter aims to rectify this gap in our knowledge by studying various 
properties of the Z{3) interface at temperatures just above the critical value. 
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So, then, what changes can we predict as the temperature drops? To begin with, a 
decrease in temperature causes a corresponding reduction in the energy penalty suppressing 
fluctuations in phase. It follows that more and larger fluctuations, bubbles of different 
phase, will appear within the two main Z{2>) domains on the lattice, following from the 
reduction in free energy of the bubble walls. Since the phase transition at the critical point 
is second-order, one also expects the expectation values of Polyakov lines in the two regions 
to converge smoothly on zero as the critical point is approached, this being the value taken 
in the disordered phase below (3c- This means that the "height" of the interface (the change 
in the real part of the Polyakov line as one moves across it) will also decrease. These 
two consequences present a double bind, with the interface decreasing in height just as the 
surrounding fluctuations increase. As the interface becomes less and less discernible amidst 
the phase turbulence, one might expect a detailed study of its structure to become well nigh 
impossible near to the critical temperature. 




p=io 




Fig. 4.1: Polyakov line (real part) profiles of the same ^(3) interface at different 
temperatures. The top left picture is identical to that of fig. 3.3, with /3 = 50 and 
a 2 X 24 X 96 lattice. Going from left to right, and then top to bottom, the other 
pictures illustrate an interface between the same phases (< L >ss e^'^*" on the 
left and the < L >?a 1 on the right) for (3 — 10, 9, 8.5, 8.3 and 8.1 respectively. 
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The pictures of Polyakov lines shown in fig. 4.f appear to confirm these fears. One 
particular Z{2>) interface is shown at a series of different temperatures, which drop towards 
the critical value. By the time (3 ~ 10, it is looking much more ragged than at the high 
temperatures of the last chapter, illustrated in the first picture. By the time /3 < 9, the 
interface is difficult to spot at all, much less study in detail. In the last picture, just below 
the critical temperature, all traces of the interface have vanished, as expected, and the 
Polyakov lines are fluctuating about zero. 

The first question to address is that of the qualitative behaviour of the interface as the 
temperature drops. Given that it must vanish at /3c, several possibilities seem to present 
themselves: 

1. The interface maintains the rigid, taut structure seen at high temperatures. It keeps its 
"area" , defined as the Polyakov "height" multiplied by the lattice width, to a minimum, 
by avoiding fluctuations in its shape. Its height shrinks steadily until it reaches zero 
at the critical temperature, where the interface consequently vanishes. Before this, 
though, it becomes invisible amidst the growing fluctuations. 

2. The interface spreads out in the longitudinal direction, following the behaviour of 
correlation lengths (such as the Debye screening length) on the lattice. These diverge as 
the temperature approaches the critical value, meaning that regions of phase separated 
by a large distance on the lattice have an increasing effect on each other. The interface 
thus reveals more and more disordered phase within, in line with the prediction of 
complete wetting, and as its height decreases and its width increases, talk of the 
interface as a physical object becomes less and less meaningful. 

3. In addition to small bubbles of phase, many more entire interfaces form across the 
lattice as the energy penalty drops. These, together with the original interface, interact 
with each other and with the fluctuating bubbles of phase, which are also forming in 
increasing numbers. These interactions, collisions and buddings off, eventually wash out 
all phase structure, leaving disordered phase throughout the lattice, with the Polyakov 
line averaging to zero. We would expect this behaviour to be suppressed by increases 
in the transverse size of the lattice, and use a variety of such sizes later in this chapter. 

4. The interface maintains a relatively narrow width, behaving essentially as a one- 
dimensional object in the spatial dimensions — a string, stretching across the lattice 
— but it suffers increasingly violent transverse fluctuations, contorting its shape and 
thereby increasing its area {i.e. its length, disregarding the time dimension). The 
fluctuations diverge towards the critical temperature, and the interface eventually 
breaks into pieces to leave only disorder. 
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As mentioned above, it is fairly clear from fig. 4.1 that the interface has vanished by the 
final picture {[3 = 8.1), but the manner of its collapse remains ambiguous. The fluctuations 
wash out the interface so much that we cannot really rule out any of the possible behaviours 
at this stage. Clearly, we need a better way to keep track of the interface as the temperature 
drops, and this is the matter addressed in the next section. 



4.2. Tracking the Interface 



/^NE way to deal with the increase in phase fluctuations is to cut out the highest frequency 
modes altogether. This can be achieved most simply by a "box-car" average, where the 
Polyakov lines of each three-by-three spatial array of points are averaged over to give a new 
value, allocated to the central point. Note that the periodic boundary conditions give the 
spatial dimensions the topology of a torus, so we have no edge effects to worry about. The 
expectation value of the Polyakov line in a region is unaffected by this averaging procedure. 



p=io 



-100 





Fig. 4.2: Smoothed Polyakov line profiles, produced by processing the raw pictures 
of fig. 4.1. As in that figure, (3 — 50, 10, 9, 8.5, 8.3 and 8.1 respectively. 
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The effect of smoothing out the fluctuations is shown dramatically in the pictures of 
fig. 4.2. These have been produced simply by applying the box-car average to the pictures 
of fig. 4.1. To make things even clearer, the Polyakov lines have first been re-binned so that 
the 24 X 96 spatial dimensions become 12 x 48, with each two-by-two array being averaged 
to one point. Now the Z{2>) interface can clearly be seen for all but the last picture, i.e. for 
all temperatures above (3c. The fact that the interface can still be seen for temperatures just 
above the critical value gives encouragement to the belief that a study of its structure should 
still be possible quite close to the phase transition, and that the shape of the interface can 
be monitored. 

However, there is a drawback with this smoothing technique: it seems less than ideal 
as a mechanism for study of the interface, as it averages out a whole class of fiuctuations 
which may be important in the interface collapse. Thus, we would like to find a better way 
to keep track of the interface. We shall, though, use smoothing later on in conjunction with 
a more reliable method of tracking the interface, and a comparison of the results will allow 
us to judge how much distortion is introduced by smoothing. 

A more satisfactory method is to produce a contour map from the real part of the 
Polyakov lines on the lattice. At its simplest, consider following a contour whose height is 
mid-way between that of the average Polyakov line at the left-hand end of the lattice and 
that at the right-hand end. This will pick out the mid-height of the interface as it goes 
across the lattice, as well as any bubbles of fiuctuating phase which are large enough to cross 
the contour height. The crucial point to note is that only the Z{'i) interface will cross the 
entire lattice, i.e. the only contour that will wrap once around the lattice in the transverse 
direction is that corresponding to the interface. In a sense, it has unit "winding number" , 
whereas any contour representing a bubble of phase will join up with itself without a net 
crossing of any lattice boundary (winding number zero). This gives us a way to tell the 
interface apart from any fluctuations, and to identify the position of its mid-height. By 
following contours at various heights between the two extremes, we can study more of the 
structure of the interface. Also, the contours allow us to measure any fluctuations in the 
shape of the interface, as well as its width in lattice units. The advantage of using contours 
is that no prior processing of the Polyakov line data is required. However, we can study 
contours of the box-car-averaged data as well if we wish, as this may be useful very close 
to the critical temperature. 



Chapter 4: Interface Properties Near the Critical Temperature 

Contours of Raw Polyakov Lines 



52 



\.' (j<ri 




40 60 

Longitudinal Direction, z 

Contours of Smoothed Polyakov Lines 




40 60 

Longitudinal Direction, z 



Fig. 4.3: Contour maps of raw and processed Polyakov lines for /3 = 10 on a 
2 X 24 X 96 lattice. These maps correspond to the /3 = 10 pictures of fig. 4.1 and 
fig. 4.2 respectively. The contour levels are at —0.25, 0, 0.25, 0.5 and 0.75, darker 
regions corresponding to higher values. 



An example of a contour map is shown in fig. 4.3, for the same gauge configuration 
as featured in the prior two diagrams for /? = 10. Notice the interface in the centre of the 
maps, represented by the only contours going ah the way across the lattice. The first picture 
shows the raw data of fig. 4.1; the second, the smoothed data of fig. 4.2, but without first 
re-binning into a 12 x 48 lattice in this case. 



4.3. Fixing the Interface Position 

("pHE ideas behind contour following also solve another problem for us: how to keep the 
interface fixed in place. In the previous chapter, we allowed it to execute a random 
walk along the lattice, as permitted by translation invariance. This allowed the system to 
tunnel into different ^(3) vacua, producing a different Z{'i) interface each time. However, 
this is unsuitable for our present purposes, as we wish to study the profile of one particular 
Z(3) interface without interference from other vacua. Of course, the properties that we 
find will apply equally to all the Z{2>) interfaces, but our measurement techniques cope best 
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with one particular interface. We need to stop the interface wandering through the twist, 
as this is the process by which the system tunnels. Ideally, the interface should be held in 
place well away from the twist, so that no part of it, even its wildest fluctuations, passes 
through the twist. 

How, then, can we pin the twist down in the longitudinal centre of the lattice? One 
method might be to use a partial anti-twist, i.e. premultiply just one or two plaquettes at 
z = L/2 by the inverse of the factor used in the twist at z = 0. If we inserted a full anti- 
twist, it would just act as another change of variables, taking the Polyakov lines which have 
been flipped by the twist back to their original value, without physical consequence; no Z['i) 
interface would need to form. However, by flipping just one or two plaquettes rather than 
the whole set across the lattice, we make it energetically favourable for the Polyakov lines 
to flip at that position, whilst still forcing the appearance of a physical interface because of 
the untwisted state of most of the plaquettes. 

The problems with this suggestion are twofold. First, an objection of principle: this 
method of pinning is distinctly artificial, breaking the translation invariance of the interface. 
How, then, can we be sure that the behaviour of the interface is not be affected by 
the procedure? The presence of the anti-twisted plaquettes, together with the periodic 
boundary conditions in the transverse direction, fixes the part of the interface at the 
boundary, suppressing the even (cosinusoidal) fluctuations in the interface in favour of the 
odd (sinusoidal) ones. Second, an objection of practice: the anti-twist makes it energetically 
favourable for the interface to form at its location, but we have already seen that energetic 
considerations become less and less important as the temperature drops towards the critical 
value, i.e. the penalty for the interface to wander away from the anti-twist becomes less and 
less. Thus, we need to introduce more and more anti-twisted plaquettes as the temperature 
drops, in order to keep the interface pinned in place. This causes more and more distortion 
in the fluctuations of the interface. In computational trials, we found that this technique 
was simply not effective at the sort of temperatures in which we were interested. 

A far more appealing approach is to locate the position of the interface after each 
sweep of the lattice, and then slide the whole gauge configuration along the lattice until the 
interface is re-centred. Any link variables which are shifted through the twist are multiplied 
by the appropriate factor, so that there is no physical effect from the procedure. This 
is clearly superior to using the anti-twist: there is no physical effect and no violation of 
translation invariance, and it can be used close to the critical temperature, for as long as 
the interface can be located. This, therefore, is the method used in all following simulations, 
unless stated otherwise. 
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Raw Polyakov Profile 



Smoothed Polyakov Profile 




40 60 
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Fig. 4.4: The procedure for locating the interface, illustrated for the (3 — \0 
configuration of previous figures, on a 2 x 24 x 96 lattice. The raw Polyakov line 
profile is smoothed, and then averaged in the transverse direction. 



4.4. Locating &; Profiling the Interface 



'"po find the position of the interface after each sweep of the lattice, we construct a 
smoothed Polyakov profile using box-car averaging, as above, and then average this 
profile in the transverse direction. This process is illustrated in fig. 4.4. The Polyakov 
line averages of this last profile are then averaged separately at opposite ends, to give the 
long dotted lines marking the top and bottom levels of the interface. Three evenly spaced 
contours are then chosen between these two, as illustrated by the short dotted lines crossing 
the interface. Starting from the (longitudinal) centre of the lattice, the algorithm searches 
for the nearest point to pass through the middle contour level in the correct direction 
(increasing from left to right). It is assumed that the Polyakov lines interpolate linearly 
between neighbouring lattice sites. The point found is defined to be the centre of the 
interface, and the configuration is then shifted to reposition this at the centre of the lattice. 
Given that this procedure is repeated after each sweep, it is extremely unlikely that the 
algorithm will be misled by a large fluctuation passing through the middle contour level; it 
will almost certainly find the true interface first, as this will be closest to the centre of the 
lattice. 
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To find a contour at a particular level, we return to the raw Polyakov line data. For 
X = 0, the algorithm follows a procedure similar to that above: starting from the centre, z = 
L2/2, it looks for the nearest point passing through the contour level. Then, it considers the 
square of points formed by the two neighbouring lattice sites at x = and the corresponding 
sites at X = 1. 



(i) 




(ii) 






+ - + 

Fig. 4.5: Options for the contour-following algorithm. Given a point at the 
desired contour level on the lower side, the Polyakov lines at the upper two 
corners determine the progress of the contour in cases (i)-(iii). A plus (minus) 
sign indicates a Polyakov line value greater (less) than the contour level. The 
same corner configurations of (iv) and (v) give two possible contour directions; 
we choose between them by considering the central average of the corner values. 
These pictures exhaust all possibilities for the values of the lower corners given, 
but any general configuration can be rotated and reflected into one of these. 



The algorithm decides to which side of the square the contour should go, as in fig. 4.5. 
Denoting Polyakov line values which are greater or less than the contour level by plus and 
minus signs respectively, the five possible configurations illustrated determine the progress 
of the contour. The last two are ambiguous, with a central average of the four corner 
values used to decide which direction of the two possibilities should be chosen, as shown. 
Once the destination side is determined, the point at which the contour crosses that side 
is found by extrapolating between the two relevant corners, as usual. Then, starting with 
this side, the whole procedure continues. Eventually, the contour must join up with itself. 
If it has wrapped across the lattice exactly once in this process, it must be the contour 
representing the interface; if not, it has merely identified a bubble of phase, and can be 
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Centre of Interface 

Fig. 4.6: An interface contour is illustrated crossing the lattice: the dashed line. 
The dotted line marks the average position of the contour, at the longitudinal 
centre of the lattice. As the contour twists and turns to cross some transverse 
lattice positions more than once, we shall approximate it, where necessary, by 
the solid contour, which crosses each transverse lattice position only once, at the 
closest approach of the dashed contour to the central (dotted) line. 



discarded. Repeated application of this technique for different contour levels determines 
the structure of the interface to the desired precision. The same procedure can be applied 
to the smoothed Polyakov line configurations, to obtain smoothed contour maps. 

In a tiny fraction of cases, the algorithm above may fail. For instance, cases (iv) and 
(v) of fig. 4.5 cannot be distinguished with certainty as Polyakov lines are only defined on 
the lattice sites. Also, extremely close to the phase transition, the interface may break up 
into several pieces, so that no contour goes all the way across the lattice. In the event of 
such failure, a reserve method is used to find the rough positions of the contours; in fact, 
this is the same method as originally used to locate the centre of the interface, but using 
raw rather than smoothed Polyakov lines: for each transverse position, x, the nearest point 
to the centre of the lattice is found at which the Polyakov line passes through the desired 
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contour level in the correct direction (i.e. increasing or decreasing with z, as appropriate 
for the interface). 

Later, we shall need to define the interface just once for each transverse lattice site, 
but at very low temperatures, the contour sometimes twists and turns so much that it 
crosses some transverse values several times. In these cases, we take the interface to be 
at that crossing point which is closest to the average longitudinal position of the interface 
across the lattice. This procedure is illustrated in fig. 4.6, with the solid contour being the 
approximation taken to the dashed contour. 



4.5. Qualitative Results for Interface Behaviour 



'T'O address the question posed in section 4.1 of how the interface behaves near the critical 
temperature, we use the contour techniques described above to follow the changes in 
its structure as the temperature drops. We start by considering some results from our 
simulations on a 2 x 36 x 72 lattice. The length of the lattice, seventy-two sites, was 
chosen after studies with various lengths, to ensure that not even the wildest fluctuations 
of the interface would reach as far as the twist during our full-length runs. The simulations 
consisted, in each case, of 2k heat-bath sweeps and 100k main sweeps. Since we were 
interested in the evolution of the interface, contour measurements were taken after every 
sweep of the lattice. Approximately 100—200 hours of CPU time on a DEC 2100 A500MP 
machine were needed for each value of /3, as shown in the following table (occasionally, parts 
of simulations were performed on DEC 3000 machines; these times have been adjusted to 
compensate): 



Table 4.1: Simulation Times for 2 x 36 x 72 lattice 


(i 


CPU time for 2k+100k Sweeps (hrs) 


8.25 


96 


8.35 


134 


8.50 


152 


8.75 


184 


9.00 


197 



In fig. 4.7, a contour snapshot of the interface is shown for each temperature studied, 
at the same point in each simulation: just over halfway through. Three contours have been 
followed in each case. We can see that the contours stay close together across most of the 
interface, indicating that the interface remains relatively narrow as the temperature drops. 
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(3=8.25 
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(3=8.35 
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p=9 Contour levels: -0.059, 0.119, 0.296 
(3=8.75 Contour levels: -0.056, 0.112, 0.279 
(3=8.5 Contour levels: -0.051, 0.102, 0.256 
(3=8.35 Contour levels: -0.042, 0.097, 0.235 
p=8.25 Contour levels: -0.037, 0.088, 0.214 



Fig. 4.7: Contour maps of the interface on a 2 x 36 x 72 lattice, for decreasing 
temperature. Moving left to right, and then top to bottom, the picture illustrates 
the interface structure around the longitudinal centre of the lattice after 2k heat- 
bath and 60k main sweeps for /3 = 9, 8.75, 8.5, 8.35 and 8.25. 
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The quantitative behaviour of the interface width will be studied more closely later in this 
chapter. However, the shape of the interface appears to fluctuate rather more wildly with 
the falling temperature. Again, a detailed study of the fluctuations will be undertaken later 
in the chapter. 



(3=8-25 




Longitudinal position, z 

Fig. 4.8: A contour map of the whole 2 x 36 x 72 lattice for (3 = 8.25, corresponding 
to the last picture of fig. 4.7, showing that only one interface exists. 



These observations suggest a behaviour of the fourth type in section 4.1, with the 
interface retaining a relatively narrow width as it approaches the critical temperature but 
suffering increasingly violent fluctuations in shape. In fig. 4.7, only the central longitudinal 
portion of the lattice is shown, to show the interface contours more clearly. However, lest 
one wonder whether more interfaces are appearing across the lattice further away from the 
one shown, we show in fig. 4.8 the full structure corresponding to the last picture of fig. 4.7. 
This demonstrates that even at /3 = 8.25, no other interfaces appear. Only bubbles of phase 
form, increasing in size and occurrence as the temperature drops and occasionally budding 
from the main interface or recombining with it. As already mentioned, the interface does 
not spread out noticeably in width, which we can define to be the average difference in 
position between the lowest and highest contour across the interface, and so it appears to 
be relatively unaffected by the divergence of physical correlation lengths such as the Debye 
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screening length. Neither does it remain in the sort of rigid shape of minimal "area" that 
we saw at higher temperatures, following from the reduction in free energy per unit area of 
the interface. 

The snapshots of fig. 4.9 for (3 = 8.5 give an idea of how the interface continually 
changes shape as it evolves in computer time. The initial configuration is shown in the 
first picture: the lattice is set up to be in two definite Z{2>) phases, with a sharp boundary 
between them midway along the lattice, as discussed in section 3.1. The configuration after 
the equilibration sweeps have occurred is shown in the second picture, with the physical 
Z{2>) interface broadening into an equilibrated shape. Pictures after this are spaced at 
intervals of 20k main sweeps, showing some snapshots of the interface as it writhes and 
twists. It is clear that the it maintains its narrow width, but fluctuates a great deal as 
bubbles of phase appear and disappear around it. Note that the interface remains centred, 
and that the fluctuations are not large enough to go beyond our field of view in the centre 
of the lattice. This reassures us that the interface does not bump into the twist at any time. 

We also follow the interface evolution by contouring the smoothed (box-car-averaged) 
configurations after each sweep, as illustrated in fig. 4.10; the pictures shown are those from 
fig. 4.9, but from the smoothed Polyakov lines rather than the raw ones. Notice that we can 
identify the interface a great deal more clearly, at the price of losing a great deal of phase 
structure. We continue to monitor the smoothed configurations, as the greater clarity of 
the interface may allow us to follow its behaviour closer to the critical temperature than 
the raw data alone. Later in this chapter, we shall have the opportunity to check how 
results from the smoothed data compare with those from the raw data, and thus gauge the 
usefulness of this procedure. The increase in computer time to keep track of the smoothed 
configurations is extremely small. 

To provide further evidence for the behaviour that we are suggesting above, we take a 
transverse-averaged Polyakov profile of the lattice after each sweep. Averaging this profile 
over all 100k main sweeps, we obtain an average profile such as that of fig. 4.11. Notice that 
this average one-dimensional profile of the interface fits the expected shape for an instanton 
interpolating between two vacua, namely a tanh function, remarkably well. This is true 
even at (3 = 8.25, though the divergence from the tanh profile is more noticeable at the 
lower temperatures; this is discussed briefly in the next section but one. 

We can determine the Debye mass of the interface from the coefficient of the fitted 
tanh function, since we expect the interface width to be of order the Debye screening 
length, viz. the inverse of the Debye mass. We can also estimate the intrinsic width of the 
interface, characterising the width at any time and point of the interface, rather than the 
Debye screening length, where all fluctuations are washed out to give the width of a "fuzzy" 
average interface. The intrinsic width is estimated by taking the average distance between 
the lowest contour, at 25% of the height of the interface, and the highest, at 75%, all the 
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Fig. 4.9: Contour snapshots of the interface at various stages oi a (3 — 8.5 run on 
a 2 X 36 X 72 lattice, showing just the central longitudinal portion of the lattice. 
Three contour levels were chosen by the method of fig. 4.4, with values as given in 
the first picture. 
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Fig. 4.10: The same contour snapshots as in fig. 4.9, but this time using the 
smoothed Polyakov lines. 
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Averaged Transverse Polyakov Profile 
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Fig. 4.11: Polyakov line profile for /3 = 8.5 on a 2 x 36 x 72 lattice, produced by 
taking an average over all sweeps of the transverse-averaged profile (see text). The 
dotted line represents the best fit to a tanh function, giving the Debye screening 
length (inverse mass); the value shown is actually (2m) ~^, as discussed in section 
2.4. The average interface width, between the lowest (25%) and highest (75%) 
contours, is also shown. 



way across the lattice, after each sweep. The result for this case is shown in fig. 4.11, and 
full results for each type of width are given later in this chapter. 



4.6. Simulation Parameters 

TTAVING established the qualitative features of the behaviour of the interface near the critical 
temperature, we now proceed with a detailed quantitative analysis of our results. Since 
we wish to examine aspects of the interface behaviour as the temperature drops to the 
critical value, we perform simulations for the various values of (3 mentioned previously; 
these are tabulated below for easy reference. It will be seen later that it is difficult for us 
to obtain reliable results closer to the critical temperature than fi = 8.25. 
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Table 4.2: Simulation Temperatures, (3 



l.2b 



^.35 



^.50 



.75 



9.00 



For each value of (3 in table 4.2, we carry out simulations with 2k heat-bath and 100k 
main sweeps on 2 x N^ x 72 lattices with a range of transverse lattice sizes, N^- This will 
allow us to examine the dependence of the interface behaviour on its size, i.e. the lattice 
width. 



Table 4.3: Simulation Lattice Widths, L^, 



18 



24 



30 



36 



42 



48 



54 



4.7. Interface Wetting 



"DEFORE moving on, we pause briefly to consider the question of interface wetting. As 
mentioned earlier, a Z{2>) interface is believed to consist of two order-disorder interfaces 
stuck back to back, enclosing a narrow slice of disordered (low-temperature) phase. Evidence 
to support this theory has been provided by simulations in 3+1 dimensions [8], measuring the 
surface tension of both types of interface at the critical temperature, and it seems reasonable 
that the same behaviour should be seen in 2 + 1 dimensions. As we do not measure the 
properties of order-disorder interfaces in this survey, it is worth considering whether we can 
see any indication of the wetting behaviour in our study of the Z{'i) interface alone. At the 
critical temperature itself, complete wetting would suggest that the interface profile should 
split in two to include an area of disordered (< L >= 0) phase. Although average profiles 
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akin to fig. 4.11 show slightly more divergence from the tanh form at lower /3, it is not 
possible to see such a splitting clearly in our system. However, we would expect any area of 
disordered phase within the interface to remain extremely small until (3 ^ f3c = 8. 175 [44], 
and so we cannot take this as evidence against wetting in 2 + 1 dimensions. 

We cannot obtain reliable average profiles much below (3 = 8.25 because of the difficulty 
of identifying the interface so close to collapse, in order to keep it centred for the average. 
Another approach to the problem which can be followed at lower /3, though, is artificially 
to split the interface into two by changing the initial conditions to include a large area of 
disordered phase between the two ordered domains. The persistence, or otherwise, of this 
interface structure can then be examined at various (3 ~ [3c- In fig. 4.12, instantaneous 
profiles are shown at fifty-sweep intervals on a long lattice at /3 = 8.15 (actually just above 
13c for our finite-sized system). Disordered phase is present only where both left- and 
right-hand plots {i.e. real and imaginary Polyakov profiles) coincide with the dotted line 
(< L >= 0: disordered phase). The disordered region given by the initial conditions is 
seen not to persist, but rather quickly dissolves into many bubbles of Z{3) phase. Again, 
this cannot be taken as strong evidence against complete wetting, because of the imprecise 
nature of the observation; but neither does there seem anything in our observations which 
could be taken as positive evidence of the wetting. 



4.8. Estimates of the Critical Coupling and the 
Magnetisation Exponent 

"\^7E now proceed to more concrete measurements of interface properties. We know that 
the expectation values of the Polyakov lines in the main phases to either side of the 
interface must tend to zero as the temperature descends to its critical value, since the 
disordered phase which exists below that temperature has the Z{3) vacuum symmetry 
restored: < L >= 0. This means that by monitoring the average values of < L > at 
either end of the lattice, which can be read off from graphs such as fig. 4.11, one can obtain 
an estimate for the critical temperature itself, and also the magnetisation exponent, (3m, 
determining the fall in the Polyakov expectation: 

<L>~(/3-/3e)^^, p>pc. 

Our estimates of the Polyakov line expectations at opposite ends of the lattice are given 
in the following tables for each simulation performed, along with an average value thence 
obtained for each value of /3: 
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Fig. 4.12: Polyakov line profiles for /3 = 8.15 on a 2 x 48 x 200 lattice, produced 
by taking instantaneous transverse averages at fifty-sweep intervals. 
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Table 4.4: Estimates of < L >l 


/3 ■•• L^ 


18 


24 


30 


36 


42 


48 


54 


Average 


8.25 


-0.15414 


-0.15975 


-0.15421 


-0.16121 


-0.16205 


-0.16053 


-0.16360 


-0.1594(37) 


8.35 


-0.18135 


-0.18355 


-0.18407 


-0.18655 


-0.18585 


-0.18511 


-0.18647 


-0.1847(19) 


8.50 


-0.20172 


-0.20304 


-0.20348 


-0.20306 


-0.20409 


-0.20379 


-0.20457 


-0.2034(9) 


8.75 


-0.22241 


-0.22277 


-0.22247 


-0.22296 


-0.22306 


-0.22286 


-0.22340 


-0.2228(3) 


9.00 


-0.23644 


-0.23666 


-0.23687 


-0.23685 


-0.23724 


-0.23715 


-0.23711 


-0.2369(3) 



Table 4.5: Estimates of < L >^ 


(3 ■•• L^ 


18 


24 


30 


36 


42 


48 


54 


Average 


8.25 


0.33667 


0.33570 


0.32767 


0.34012 


0.33955 


0.33600 


0.33651 


0.3360(41) 


8.35 


0.37655 


0.37486 


0.37494 


0.37491 


0.37511 


0.37453 


0.37612 


0.3753(7) 


8.50 


0.40925 


0.40898 


0.40797 


0.40816 


0.40895 


0.40908 


0.40855 


0.4087(5) 


8.75 


0.44595 


0.44575 


0.44549 


0.44665 


0.44629 


0.44620 


0.44635 


0.4461(4) 


9.00 


0.47399 


0.47416 


0.47412 


0.47340 


0.47391 


0.47452 


0.47407 


0.4740(3) 



The figures from each simulation come from the transverse-averaged Polyakov profiles, 
averaging over all 100k main sweeps, and then averaging over the six longitudinal points at 
each end of the lattice (0 < z < 5 and 66 < z < 71). The errors from a given simulation are 
O(10~^), decreasing to O(10~^) for the larger (5 values, and are not shown, being negligible 
compared to the error between different simulations, which is included with the average. 

An accurate fit to the form 



<L>= aiP - /?e) 



Pm 



is difficult, since we only have five points and must fit three parameters. Using the gradient 
expansion fitting of IDUs function "CURVEFIT" , we plot the average values from the tables 
in fig. 4.13. The independent estimates from each of the two Z{2>) phases give the following 
results for the critical temperature and magnetisation exponent: 

/3c = 8.170(5), 



/3m = 0.159(3). 



Chapter 4'- 


Interface Properties Near th 




Linear Plot (L.H. Phase) 


-0.14 




-0.16 


: A - 


^ -0.18 

V 

§, -0-20 


■ A-, '- 

■ A,, : 


-0.22 


■ "A-,, : 


-0.24 


■ ^~VA 


8.2 8.4 8.6 8.8 9.0 



68 



-0.60 - 

7- -0.65 : 

J 

V 

^ -0-70: 

I -0.75 : 

-0.80 L 



Log Plot (L.H. Phase) 

Best Fit: |3 =8.176+0.001 
|3m=0.163+0.001 



-1.2 -1.0 -0.8 -0.6 -0.4 -0.2 -0.0 
log,o(P-Pc) 



Linear Plot (R.H. Phase) 



0.50 r 



_□ 



0.45: 

A 

V 0.40- 
0.35- 



a- 



-0.30 r 



,« -0.35 - 



B^ 



Q 



A 

J 

V 



o 



Log Plot (R.H. Phase) 

Best Fit: (3=8.163+0.005 
|3m=0.156+0.003 



0.30 1 



8.2 



8.4 



8.6 



9.0 



-0.40 : 

-0.45 : 

-0.50 L_ 
-1.2 




-1.0 -0.8 -0.6 -0.4 
logio(P-Pc) 



-0.2 -0.0 



Fig. 4.13: A three-parameter fit for the critical temperature, /3c, and the 
magnetisation exponent, /3m, from the expectation value of the Polyakov line deep 
in the phases at the right- and left-hand ends of the lattice (denoted by "R.H." 
and "L.H." respectively), as found from averaged profiles such as fig. 4.11. 



These are some distance from the values of Pc ~ 8.14 and Pm ^ 0.13 obtained in [45], 
the estimate for Pc being closer to the continuum estimate of 8.175 than to the value of 8.14 
which broadly corresponds to our lattice sizes. As a check, fig. 4.14 shows our estimates for 
Pm if we assume that Pc = 8.14; the estimates are little changed. 



4.9. Quantitative Results 



TTAVING established that the interface remains as a thin, localised object even very close to 
the critical temperature, we now wish to examine its behaviour in more quantitative 
detail. To do this, we need a way to quantify the fluctuations in the interface. The phase 
transition is second-order, with the Polyakov lines (and the interface) vanishing smoothly 
at Pc', thus, we expect the fluctuations to diverge as P ^ Pc, and we wish to quantify this 
divergence. 
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Fig. 4.14: A two-parameter fit for tfie magnetisation exponent, /3m, using /3c 
8.14. 



Illustrated in fig. 4.15 is a contour going across the lattice, with displacements from its 
average position at transverse location x given by (t){x). To characterise the fluctuations, we 
monitor the fluctuation "moments", where the nth moment is defined to be < (j){x)^ >. We 
track the first six moments (1 < n < 6) in our simulations; higher moments are subject to 
unacceptably large errors. We take a separate average for each transverse position across the 
lattice, and, because of the translation invariance of the interface, we expect the moments 
to be roughly constant across the lattice after a large number of sweeps. 

First of all, we need to examine the qualitative behaviour of the fluctuation moments. 
As mentioned above, we expect the moments to be essentially constant across the lattice, 
owing to the periodic boundary conditions in the transverse direction and the translation 
invariance of the interface. Of course, we would expect the moments to approach constancy 
more closely with increasing simulation length. 

The following graphs show some moments for (3 = 8.5, averaged over the full 100k 
sweeps. Separate sets of data are shown for each of the three contours followed, at 25%, 50% 
and 75% of the height of the interface respectively. The qualitative behaviour is roughly the 
same for each contour. The upper and lower contours are subject to rather more fluctuation 
than the middle one, as we expect since they are closer to the outside of the interface and 
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Fig. 4.15: An interface contour is illustrated crossing the lattice: the bold line. 
The dotted line marks the average position of the contour, at the longitudinal 
centre of the lattice. We quantify the fluctuations in the interface shape by the 
displacement, 4>{x), from the average position at each point x across the lattice. 



therefore more susceptible to bubbles of phase forming nearby and distorting the edge of the 
interface. The even moments are fairly constant across the lattice, whilst the odd moments 
are very much smaller and show more variation. Note that the first moment must average 
to zero across the lattice, by definition, and we likewise expect all odd moments to be 
comparatively close to zero, becoming closer with even longer runs. 

In fig. 4.17, we show the moment data from the smoothed contours measured in the 
same simulation as fig. 4.2. One can see that the behaviour is very similar to the raw data, 
both qualitatively and quantitatively. Since this may indicate that the smoothed data has 
not lost a significant amount of information about the fiuctuations, we shall present our 
data later on for both raw and smoothed moments, in order to compare the two. 
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Fig. 4.16: The first six moments for [3 — 8.5 on a 2 
separately for each of the three contour levels followed. 
100k sweeps for the raw Polyakov data. 
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Fig. 4.17: The first six moments for /? — 8.5 on a 2 x 36 x 72 lattice, sliown 
separately for each of the three contour levels followed. These are averages over 
100k sweeps for the smoothed Polyakov data corresponding to fig. 4.16. 
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4.10. A Toy Model: A Scalar Field Theory with 
Associated Feynman Rules 

QlNCE the interface has been shown to behave as a fluctuating, narrow object, winding 
across the lattice, it may be helpful to think about it as a one-dimensional object, viz. a 
string. At high temperatures, the string stretches tightly across the lattice, its length being 
equal to the transverse size of the lattice. At low temperatures, though, if coils and twists 
ever more violently, and loops form which break off and recombine with it, representing 
bubbles of Z{'i) phase. If we take this picture further, we can consider the possibility of 
describing the string in terms of a one-dimensional field theory. The scalar field in this case 
is simply the displacement of the string from equilibrium at a position x along its length: 
(j){x). At high temperatures, we essentially have a free field theory; at low temperatures, 
we need an interaction term to account for the breaking off and recombination of loops of 
string. 

To preserve the translation invariance of the interface, the Lagrangian cannot depend 
directly on (j){x), but instead must depend only on the derivative dx4'{x). We keep the 
normal kinetic term for a scalar field theory, and add a quartic interaction term; a cubic 
interaction is ruled out by "parity" , as we expect the interface to be invariant under (j) ^ —cp. 
Thus, we have a Lagrangian of the form 

£ = i7(5.<A)' + ^a#, (4.1) 

where 7 and A are both unknown functions of (3. By examining the dependence of the 
fluctuations on the temperature, we shall be able to determine these functions. Since we 
expect the fluctuations to diverge as /3 ^ /?c, this implies that 7 and A will both be functions 

of(/3-/3c). 

To test how well this toy model describes the actual behaviour of the Z{'i) interface, 
we need to calculate the correlation functions of 0, which we can then compare with the 
behaviour seen. The n-point connected vacuum correlation function at point x, < 0(x)" >, 
corresponds to the nth Wick-subtracted fluctuation moment at the same point. The Wick 
subtracted moments, denoted by the subscript "Ws", are given, as is well known, by the 
following relations: 



2nd Moment: 


< (t? >Ws =< <f >, 




4th Moment: 


< (j) >Ws =< 4> > " 


-3 < 0^ >^ 


6th Moment: 


< (j)^ >Ws =<</>*">■ 


-15 < 02 >3 



■15 < (t? X 0'* >ws ■ 

For the Lagrangian above, the odd correlation functions are zero, so we expect all the odd 
moments to be zero. To calculate the even moments, we need to calculate the appropriate 
vacuum diagrams. 
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Fig. 4.18: The vacuum diagrams corresponding to the 2nd, 4th and 6th connected 
vacuum correlation functions respectively. The open circle simply denotes point 
X, the endpoint of each external leg; the solid circle denotes an interaction vertex, 
A; and the momentum travelling along a leg is denoted by k. 



These diagrams are shown in fig. 4.18, with the momenta denoted by the /c's, and the 
four-way vertices marked by A, corresponding to the interaction term in the Lagrangian. 
The diagrams shown are the lowest-order Wick-subtracted diagrams with two, four and 
six legs respectively leaving the point x. To perform the calculations, we need to derive 
Feynman rules for a {dxCJ))^ theory, rather than the usual (j)'^, and we do this now. 

The path integral, for our one spatial dimension only, is 

W[J] =N f -p^pj M^i-P^9.'P)+J'P) ^ 



where N is the normalisation and C is given by (4.1), written in the form 

In the absence of the interaction, 

W[J] = Wo[J] = Arel/^^'^^"'(^)^^(^-^)-^(^), 

with the Feynman propagator for the scalar field, 

AF{x-y)= / ^^^^. 
J 27r7 p^ 

To include the interaction, we expand W[J] as a power series in A: 

W[J] = e'^''^'-'^^'^Wo[J] 

6 



1 + 
Wo[J] + 



A 



4! 



dx 



6J{x 
d 6 



+ ... 



dx 6J{x) 



Wo[J] 
Wo[J] + .... 
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Now, 

^^ = j dyd^[^F{x - y)]J{y)Wo[J], 

and further differentiation of this yields 



W^[J] = {l + ^/ 



dx ?,{dl^F{^)f 

+6Jdyidy2dl[AF{0)d:,[AF{x-yi)d,[AFix-y2)]]]J{yi)J{y2) 
+ / dyidy2dy3dy4dx[AF{x - yi)9^[Air(x - 2/2) 

X a^[Ai7(x - y3)d^[AF{x - ^4)]]]] J(yi)^(y2) ^(^3)^(^4) 



+ ... 



Wo[J]. 



The free-field Green's functions are defined by 



G^'^Hx, 



, . . . , Xji 



,5J(xi)...5J(x„)i-^=o' 
giving the propagator 

G^'^\xi,X2) = Air(xi - X2). 

The connected Green's functions are defined by 



<5"X[J] 

• • ) Xn) — YTT 

dJ{2 

giving interaction term 



'^""(^■- ■ ■ ■•'•'^ = ^J(:,,)..^V) '^-°- ^■'" = '° "'■''■■ 



G('^)(xi,X2,X3,X4) = -^ / (ia;5^[Ai7(x - Xi)9^[Air(x - X2)9a;[Ai7(x - X3)d^[AF{x -X4)]]]] 

= A / dx{-d:ci [Af(x - xi)])(-aa;2 [Af(x - X2)]) 

X (-5^3[Af(x - X3)]){-d^^[AFix - X4)]). 
Converting to momentum space, 

G'-^\pi,P2,P3,P4){'^-^)^Hpi +P2 +P3+P4) 

= /"dxidx2dx3dx4G(^Ha;i,X2,X3,X4)e-^(Pi^i+P^^^+P=^^+P^^^\ 

and, integrating by parts, 

dxid^^Apixi - x)e-*Pi^i = ipi j Af{xi - x)e-'P^''\ 

Thus, each vertex picks up a factor of 

G^'^\pi,P2,P3,P4:) = >^PiP2P3Pi- where Pi + P2 + P3 + Pi = 0. 
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4.11. Field Theoretic Predictions for Fluctuation 
Moment Behaviour 

TN the continuum, therefore, the Feynman diagrams of fig. 4.18 are given by the foHowing 
formulae: 



27r7 J k^ ' 



<<P[x) >ws = j^j^^,J ^...jr^ 



X 2'Kkik2k^k4^5{ki + k2 + k-^ — ^4) 
A 1 /" dki dk2 dks 1 



\6 



4174 (27r)3 J fci fc2 ^3 {ki +k2 + ks) ' 

X {2it)'^ kik2ksk4k5kQk'^5{ki — k2 — k^ — /c7)(5(/c4 + k^ — k^ + k-j) 
A^ 1 /" dk2 dks dk^ dk^ ,,1 1 

~ 2!^(27r)5 J "fc7"fc7 ^4 ^5 ^ (^2 + fcs + fc?) (fc4 + ^5 + kj) ■ 

Note that there is no exponential in the first formula, as the propagator starts and 
finishes at the same point, x. Also, the symmetry factor in the last formula is from ^ x 
(il)^ X (4 X 3 X 2)^ = 21, from the two vertices, the four legs per vertex, the number of legs 
connected to x, and the number connected within the loop respectively. 

Since we are working on a finite space of extent L^, we have to replace the integrals 
over k with sums over integers n, putting k = 2Tm/Lx'- 

27rJ L, ^ 

n=l 

Using the numerically-calculated results 

E— = — T^ 111 1 -649(1) 



I'a •• • "■c 



v? 6 

n 

and y —- — — , ^\ w ^\ , ^60(1), 

/ ■' ria rib ric rid ("a+n-b+n-c) (na+nd+n^) \ -" 

ria-.-rie 

we obtain the following expressions: 



< cj){xf >ws = (6.49(1)) ^ ^^ 

< 4>ixf >ws = (60(1)) 



4! (27r)474 ' 

2r(2^ov' 
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Note that the above expressions for the correlation functions do not depend on the 
position x, so we now write (p{x) simply as (p. This confirms our original assertion that the 
fluctuation moments should be constant across the lattice. In future, we shall work with 
the average of a fiuctuation moment across the lattice. Now, since the action, S = J Cdx, 
must be a dimensionless quantity, and assuming 7 to be dimensionless, we find that must 
have space dimension 2, and coupling A dimension 1. Thus, as the only fundamental length 
in our space is L^, this implies that A ~ L^, or A = XqLx- Substituting this dependence 
into the previous expressions, we obtain the following predictions: 

< 0^ >~ Lj., < 4>'^ >vKs~ LI and < (j)'^ >ws^ L^. 



Checking the La;-dependence of our results will thus help to verify whether the scalar field 
theory really is a good model for the interface fluctuations. However, the above calculations 
have all been for the continuum, so we flrst need to calculate some corrections. 

We have neglected the fact that we are working on a lattice with site spacing a and 
periodic boundary conditions. To take this into account, we need to replace k^ with 
4 sin {ak/2). This is equivalent to the following change in our sums over n: 

A;2 \2TmJ Asin^iirna/Lx) 4sin^(7rn/iV^) ' 

We now calculate the corrections to our previous expressions which result from this change. 
First, we tackle the two-point function. To do this, we calculate the sum from an 
integral, working with an expansion in the small parameter ^■. 






1 



sin X 



Consider the n 



1 term: 

TV 



dx 



1 



1 



sin X 



dx 



Nx 1 
2 -" 

n=l 



1 



sm 



2 f nn 



+ 



d^. 



'0 sm'(ii+0 
Similarly, for the n > 1 terms, 
-1 



d^ 



vr/iV, 



sin^(7r/iV^ 



^ + Ol{n/N,n 



2 

E 

n=2 



1 



sin2(^+e)" 



n=2 

+ 3? 



CSC 






. vrn. 



2Ccot( — jcsc- 






cot^ 



TT 



2 



E ^«^'( 



n=2 



vrn 

1%' 



I CSC 



.vrn. 



+ -CSC 



, vrn 



vr 



2 2 

E 

ra=2 



1 
cot( 



+ 



dc 



vrn. 



IT 



3 2 ^ 

E 

n=2 



COt^ 



, vrn. 



CSC 



, vrn. 



H — CSC 
3 



I CSC 



, vrn. 



, vrn. 



+ ... 



JVx 



-1 



vr 



E 

n=2 



CSC 



, vrn. 



+ 



N^. 



vr 



77,3 ^ 77,4 



+ 0[{7r/N, 



Chapter 4: Interface Properties Near the Critical Temperature 78 

Now, 

n=2 n=2 

Thus, our result for the original integral, dropping 0{{t: /Nx)"^] terms and above, is 

!hL-l j^ _2 

n=l n=l 

This series expansion in N~'^ must be equal to the value of the original integral, 

[-COtxl „"^ =2cot-^ =2(^ -TTTT + •••), 

so our final result is that 

1 ^2 ^^ ^^ X 2 X 2 
> CSC ss — 1 . 

The correction to our original prediction, that < (jp' >^ N^, is negligible for our purposes 
(about 0.3% for A^^ = 18). We assume that the corrections to the higher moments are 
similarly negligible, as the calculation of these lattice corrections is considerably more 
difficult. 



4.12. Observed Behaviour of Fluctuation Moments 

T^T^ now compare these predictions with the results that we obtain from our computer 
simulations. Tables 4.6-8 list a sample set of results for a contour halfway up the 
interface, using raw (i.e. unprocessed) data. An estimate of each even moment is presented 
for each value of f3 and L^, together with the calculated gradients of log-log plots of the 
moments against (/3 — (3c) or L^ as appropriate, which give the power-law dependence 
of the moments on these quantities. To examine the behaviour of the interface more 
fully, we tracked fluctuations in a total of three contours, defined to be equally spaced 
between the top and bottom of the interface. For each of these three contours, we also 
measured the fiuctuations in the interface after smoothing, to see whether the procedure 
would affect results. Appendix B lists complete results for all three contours used, and 
both raw and smoothed Polyakov data. In this chapter, tables 4.9-11 summarise the power- 
law dependence of the moments from each dataset on L^, while tables 4.12-14 give the 
equivalent dependence on (3 — (3c- 
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Table 4.6: Estimates oi < (j)^ > from Raw Binned Data, Middle (50%) Contour 




P'--L^ 


18 


24 


30 


36 


42 


48 


54 


Slope 


8.25 


6.871°;?^ 


8.27;°;J1 


12.31°:^ 


14.41°:! 


i7.8i°:l 


20.81°:^ 


23.11°:? 


1.20+0.03 


8.35 


5.02t°;^? 


6.7ltH'5 


9.14tH^ 


11.41°:! 


i3.4i°:i 


i6.5i°:l 


i8.2i°:i 


1.21+0.01 


8.50 


q oQ+0.03 
"J-y'J-0.03 


c 1Q+0.03 
O-J-y-0.04 


7 17+0.05 
'•-L'-0.05 


8.89in^ 


io.ol°:l 


ii.9l°:l 


i3.8l°:l 


1.15+0.01 


8.75 


o Q7+0.02 
'^•'-"-0.02 


A 91+0.03 

^•^-^-0.03 


c- 09+0.04 
"-••"^^-0.03 


6.34lH^ 


7.63lH^ 


oc:+0.07 
»-»O_0.05 


9.571°:°^ 


1.05+0.01 


9.00 


9 fiO+O.Ol 
■^•0^-0.02 


44+0.02 
'^•^^-0.02 


A 00+0. 03 
^••30-0.02 


c 97+0.03 
"-••^'-0.03 


fi 90+0.05 

"•^^-0.03 


7 07+0.05 
'•"^'-0.05 


7 sq+0.04 


1.03+0.01 


Slope 


-0.41+0.01 


-0.37+0.02 


-0.44+0.02 


-0.43+0.03 


-0.44+0.02 


-0.45+0.03 


-0.46+0.04 




1 





Table 4.7: Estimates of < cf)'^ >Ws from Raw Binned Data, Middle (50%) Contour 




I3'--L,, 


18 


24 


30 


36 


42 


48 


54 


Slope 


8.25 


2270lt«° 


3341!^^ 


mtlf 


463l^f 


4921^1 


585_i3 


690l?o 


-0.54+0.27 


8.35 


775l?f 


63.11?:^ 


111 + 2 


1661° 


1991° 


344124 


2991" 


1.31+0.49 


8.50 


19.51°:^ 


28.3_o:8 


59.71?:^ 


09 cr+3.3 


88.01^6 


1261? 


16311 


1.95+0.05 


8.15 


9.92l°:f7 


16.71°:^ 


22.3l°:« 


30.31^:^7 


43.211:^ 


48.6li:3 


49.51- 


1.54+0.04 


9.00 


5.85_o:i7 


8.80_o:23 


10 1+0.5 

J- "^•-'--0.3 


15.61°:° 


23.31°:^ 


28.61^:^ 


28.611:^ 


1.54+0.03 


Slope 


-2.40+0.44 


-1.33+0.12 


-1.53+0.17 


-1.47+0.08 


-1.29+0.05 


-1.26+0.12 


-1.32+0.11 




1 





Table 4.8: Estimates of < (/)° >ws from Raw Binned Data, Middle (50%) Contour 




^'■■L,. 


18 


24 


30 


36 


42 


48 


54 


Slope 


8.25 


5720000 

+ 1370000 
-510000 


391000 

+441000 
+45000 


1280000 

+720000 
-180000 


297000 

+320000 
-16000 


199000 

+295000 
-23000 


24900 

+ 13600 
-3800 


71700 

+112000 

-7300 


-4.40+0.39 


8.35 


1880000 

+970000 
-20000 


1830 

+260 
-170 


3710 

+640 
-330 


5150 

+700 
-560 


6180 

+ 1300 
-980 


42800 

+41100 
-2700 


2640 

+2570 
-1400 


-2.87+1.05 


8.50 


2891^^ 


380lt? 


12601-° 


21601?;,° 


9391^1^ 


1660t^^° 


968l°«° 


1.82+0.26 


8.75 


1091^2 


1671^4 


noli 


98.011^7 


2851-^ 


77 c; + 161.1 
' '-^-61.2 


-6251^1^ 


0.44+0.25 


9.00 


27.ll|^ 


48.21^:^9 


39.411?:° 


-19.71?°:^ 


-65.213^:7 


-2131^^ 


-1251^1 


1.06+0.14 


Slope 


-5.57+0.98 


-3.16+0.55 


-3.95+0.47 


-3.33+0.37 


-3.15+0.23 


-2.27+0.82 


-2.98+0.29 




1 
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Table 4.9: Power-Law Dependence of < 0^ > on L^. 




Contour 


Data 


f3 = 8.25 


(3 = 8.35 


(3 = 8.50 


P = 8.75 


P = 9.00 
















Lower 


Raw 


1.30±0.02 


1.22±0.01 


1.13±0.01 


1.04±0.01 


1.02±0.01 


Lower 


Smooth 


1.21±0.02 


1.16±0.02 


1.14±0.02 


1.09±0.01 


l.lOiO.Ol 
















Middle 


Raw 


1.20±0.03 


1.21±0.01 


1.15±0.01 


1.05±0.01 


1.03±0.01 


Middle 


Smooth 


1.26±0.01 


1.23±0.01 


1.17±0.01 


l.lliO.Ol 


1.13±0.01 
















Higher 


Raw 


0.91±0.03 


1.08±0.02 


1.04±0.01 


0.92±0.01 


0.89±0.01 


Higher 


Smooth 


1.23±0.01 


1.12±0.01 


1.03±0.01 


0.94±0.01 


0.95±0.01 








Table 4.10: Power-Law Dependence of < 0'* >ws on L^ 




Contour 


Data 


(3 = 8.25 


(3 = 8.35 


(3 = 8.50 


13 = 8.75 


P = 9.00 
















Lower 


Raw 


0.27±0.25 


1.84±0.47 


1.99±0.05 


1.58±0.03 


1.51±0.03 


Lower 


Smooth 


-0.18±0.15 


-0.11±0.38 


1.65±0.07 


1.21±0.05 


1.19±0.05 
















Middle 


Raw 


-0.54±0.27 


1.31±0.49 


1.95±0.05 


1.54±0.04 


1.54±0.03 


Middle 


Smooth 


-0.56±0.17 


-0.02±0.31 


1.68±0.07 


1.17±0.06 


1.27±0.05 
















Higher 


Raw 


-1. Olio. 24 


0.88±0.54 


1.62±0.04 


1.32±0.03 


1.28±0.02 


Higher 


Smooth 


0.08±0.16 


0.24±0.33 


1.36±0.05 


1.06±0.04 


1.02±0.03 








Table 4.11: Power-Law Dependence of < (f)^ >Ws on Lr^ 




Contour 


Data 


(3 = 8.25 


(3 = 8.35 


P = 8.50 


P = 8.75 


P = 9.00 
















Lower 


Raw 


-4.38±0.24 


0.04±1.18 


1.60±0.29 


0.73±0.17 


-0.52±0.11 


Lower 


Smooth 


-1.62±0.20 


-5.38±1.17 


1.19±0.32 


1.02±0.09 


0.90±0.08 
















Middle 


Raw 


-4.40±0.39 


-2.87±1.05 


1.82±0.26 


0.44±0.25 


1.06±0.14 


Middle 


Smooth 


-2.55±0.24 


-5.29±1.10 


1.66±0.22 


0.92±0.08 


?±? 
















Higher 


Raw 


-4.40±0.39 


-3.69±0.99 


1.54±0.27 


0.94±0.26 


0.60±0.19 


Higher 


Smooth 


-2.60±0.19 


-4.14±0.86 


1.56±0.18 


1.21±0.10 


1.03±0.07 
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Fig. 4.19: Graphs of even moments against transverse size, Lx, for /3 = 8.5. Data 
are shown for both raw (unprocessed) and smoothed Polyakov configurations. For 
the 4th and 6th moments, the left-hand graphs are as measured initially; the right- 
hand ones are Wick-subtracted. 
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1.2 



1.0 



a 0.8 



0.6 



0.4 



Slope(U)= 1.15+ 0.01 
Slope(S)= 1.17+0.01 




1.20 1.30 1.40 1.50 1.60 1.70 1.80 
L^-Averaged 4th Moment (Wick-Subtracted) 



V 1 



2.4 


Slope(U)= 1.95+ 0.05 


- 


2.2 


Slope(S)= 1.68+ 0.07 




2.0 


y 


4 : 
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y 


- 
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- 
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^ j^ 


_ 




*■ 


- 
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1 |- 
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Unprocessed 
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Fig. 4.20: Log plots corresponding to the Wick-subtracted data of fig. 4.19, with 
best fits to a power-law dependence on the lattice width. 



A typical set of results is shown in fig. 4.19, showing the dependence of < (/)(x)'^ > 
on the lattice width (interface length), L^, for a particular value of /?. Both raw and 
smoothed data are shown, and these are in good qualitative agreement, though there is 
a small quantitative difference between them. In fig. 4.20, we obtain best estimates for 
the power-law dependences of the even Wick-subtracted moments on the lattice width. 
These are displayed in the graphs, matching the values given in the earlier tables, and show 
reasonable agreement with our theoretical predictions for the 2nd and 4th moments. The 
agreement is distinctly less good for the 6th, but the nature of this Wick-subtracted moment, 
the small difference between two large numbers derived from powers of a fluctuation, may go 
a long way towards explaining this discrepancy. Our results are similar for the other values 
of /3, as shown in tables 4.12-14. Though not in perfect agreement with our predictions for 
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the higher moments, these measurements encourage us to beheve that our toy Lagrangian 
is a reasonable model for the interface behaviour. 

It should further be pointed out that some of the results at low (3 and small L^ seem 
to suffer from anomalous large "spikes" in the data, which can skew the estimates of the 
higher moments severely, and with them the slopes. This problem is addressed at the end 
of the next chapter, and oddities such as the negative slopes of the higher moments against 
Lx for small (3 then disappear from the corrected results given there. 

Having obtained reasonable results for our model, we now turn to the question of the 
unknown functions, 7 and A. To tell anything about these functions, we need to produce 
similar plots to fig. 4.19 and fig. 4.20, but showing the dependence on (3 instead of L^. In 
fig. 4.21, we show the divergence of the non- Wick-subtracted even moments as /3 ^ (3c- The 
data shown are for lattices with L^ = 36. To find the degree of divergence, we need to fit 
this behaviour to a power law dependence on {(3 — f3c)- 

In fig. 4.22, we fit our results to just such a power-law divergence, taking [3c = 8.175. 
We obtain good fits even for the sixth moment, giving another indication that our model 
is reasonable. For the data shown in the graphs of this chapter (L^ = 36, middle contour, 
raw data), the fitted slopes for the second, fourth and sixth moments are: 

< Ct? >~ (/3-/?,)-0-43(3)^ < ^4 >^^^ (^_^J-1.47(8)^ < ^6 >^^^ (/3 -/?,) -3-33(37) _ 

Referring back to our formulae for the moments, we see that the first two imply the following 
for our unknown functions of (3: 

These two together predict, from our previous formulae, that 

within roughly one standard deviation of the observed behaviour. Our other sets of data 
give similar results for 7 and A. 

It can be seen from the tables that the different contours give similar results, as hoped. 
For this reason, we use the middle contour, as it should be less susceptible to the influence 
of nearby bubbles of phase. Also, the results after smoothing are close to those from the raw 
data, but the discrepancies increase for the higher moments. Smoothing appears to be a 
useful technique for keeping track of the interface at low /?, indeed lower (3 than is possible 
in its absence. It appears not to alter the results signiflcantly for the lower moments, 
though the discrepancies increase fast enough for the results to disagree somewhat with the 
unprocessed results over the higher moments. 
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Table 4.12: Power-Law Dependence oi < (p^ > on (/3 — 13c) 




Contour 


Data 


L^ = 18 


L. = 24 


L^ = 30 


L^ = m 


L^ = 42 


L^ = 48 


L^ =54 




















Lower 


Raw 


-0.33±0.02 


-0.32±0.02 


-0.39±0.02 


-0.40±0.03 


-0.41±0.01 


-0.42±0.02 


-0.44±0.03 


Lower 


Smooth 


-0.51±0.03 


-0.46±0.03 


-0.51±0.03 


-0.50±0.03 


-0.50±0.02 


-0.51±0.02 


-0.53±0.03 




















Middle 


Raw 


-0.41±0.01 


-0.37±0.02 


-0.44±0.02 


-0.43±0.03 


-0.44±0.02 


-0.45±0.03 


-0.46±0.04 


Middle 


Smooth 


-0.50±0.03 


-0.48±0.03 


-0.54±0.03 


-0.52±0.04 


-0.52±0.02 


-0.53±0.03 


-0.56±0.03 




















Higher 


Raw 


-0.41±0.03 


-0.37±0.02 


-0.42±0.02 


-0.39±0.03 


-0.40±0.02 


-0.41±0.03 


-0.42±0.04 


Higher 


Smooth 


-0.40±0.03 


-0.42±0.03 


-0.47±0.03 


-0.48±0.03 


-0.49±0.02 


-0.51±0.02 


-0.54±0.03 








Table 4.13: Power-Law Dependence of < cj)'^ >ws on (/? — (5c) 




Contour 


Data 


L^ = 18 


L, = 24 


L^ = 30 


i^. = 36 


L. = 42 


L^ = 48 


L^ =54 




















Lower 


Raw 


-2.00±0.43 


-1.25±0.10 


-1.40±0.06 


-1.23±0.08 


-1.19±0.04 


-1.24±0.10 


-1.29±0.11 


Lower 


Smooth 


-2.69±0.34 


-1.63±0.15 


-1.84±0.14 


-1.67±0.10 


-1.54±0.06 


-1.65±0.12 


-1.76±0.11 




















Middle 


Raw 


-2.40±0.44 


-1.33±0.12 


-1.53±0.17 


-1.47±0.08 


-1.29±0.05 


-1.26±0.12 


-1.32±0.11 


Middle 


Smooth 


-2.76±0.25 


-1.80±0.17 


-1.88±0.17 


-1.74±0.10 


-1.61±0.04 


-1.62±0.13 


-1.82±0.11 




















Higher 


Raw 


-2.33±0.52 


-1.20±0.25 


-1.33±0.24 


-1.21±0.07 


-1.11±0.05 


-1.18±0.07 


-1.26±0.06 


Higher 


Smooth 


-2.28±0.34 


-1.44±0.13 


-1.38±0.15 


-1.43±0.05 


-1.34±0.04 


-1.45±0.06 


-1.50±0.02 








Table 4.14: Power-Law Dependence of < (f)^ >Ws on (/3 — (ic) 




Contour 


Data 


^:. = 18 


L. = 24 


L^ = 30 


L^ = 36 


Lx = 42 


L, = 48 


i^. = 54 




















Lower 


Raw 


-4.99±1.05 


-3.49±0.57 


-4.08±0.37 


-3.11±0.38 


-2.75±0.12 


-1.98±0.47 


-3.00±0.32 


Lower 


Smooth 


-5.40±0.94 


-3.30±0.47 


-3.96±0.37 


-3.25±0.48 


-3.45±0.33 


-3.47±0.51 


-3.55±0.22 




















Middle 


Raw 


-5.57±0.98 


-3.16±0.55 


-3.95±0.47 


-3.33±0.37 


-3.15±0.23 


-2.27±0.82 


-2.98±0.29 


Middle 


Smooth 


-5.78±0.92 


-3.35±0.59 


-3.79±0.52 


-3.24±0.32 


-3.39±0.25 


?±? 


-3.71±0.17 




















Higher 


Raw 


-5.36±1.01 


-2.22±0.67 


-2.58±0.70 


-2.39±0.24 


-1.99±0.07 


-2.36±0.23 


-2.85±0.11 


Higher 


Smooth 


-5.25±0.99 


-2.74±0.43 


-2.43±0.57 


-2.44±0.22 


-2.29±0.10 


-2.46±0.17 


-2.36±0.19 
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Fig. 4.21: Graphs of even moments against /3, for Lx — 36. Data are shown for 
both raw (unprocessed) and smoothed Polyakov configurations. For the 4th and 
6th moments, the first graphs are as measured by the computer, and the second 
graphs are Wick-subtracted. 
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Fig. 4.22: Log plots corresponding to the Wick-subtracted data of fig. 4.21, with 
best fits to a power- law dependence on (/3 — /3c). 



4.13. Interface Width: Intrinsic &; Screening 



T N addition to our survey of fluctuation moments, we also measure the Debye electric 
screening mass and interface width for each value of (3 and L^, as illustrated in fig. 4.11. 
As touched on previously, we can measure the "intrinsic" width of the interface by taking 
the distance between the lower and upper contours at every point across the lattice, and 
then averaging this quantity across the lattice, and then over all sweeps. By contrast, 
we obtain the screening mass by averaging the transverse-averaged Polyakov profile of the 
interface, obtained after each sweep, over all sweeps, and fitting the resulting function to a 
tanh with coefficient proportional to the Debye mass. The inverse of this mass is the Debye 
screening length, giving the screening width of the interface, in contrast to the intrinsic 
width mentioned before. 
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Fig. 4.23: Measurements of the average interface width in units of the lattice 
spacing. For reasons of accuracy, twice the half-width (the average separation of 
the 25% and 50% contours) is plotted rather than the width (25% — 75%). This 
remains roughly constant as the width of the lattice increases, and rises somewhat 
as the temperature drops. 



We can return to our picture of the interface as a string. The intrinsic width is the 
width of the string itself, which should be independent of the length of the string (L^;), and 
which we would not expect to change much with temperature. To visualise the screening 
width, however, we need to imagine letting the string fluctuate very fast while we look at 
it, so that the different shapes blur into one another; at very high temperatures, where 
the fluctuations are suppressed, this blurred string should look only a little wider than the 
stationary string, but as the temperature drops and the fluctuations diverge, we expect 
the blurred string to appear wider and wider, also diverging in width towards the critical 
temperature. Since the mean square fluctuations (the second moments) grow like L^,, we 
might expect the screening width of the interface to grow like \/Z^. 

The results for the intrinsic width, shown in fig. 4.23, show that the width does 
indeed remain constant as L^ increases, and that it increases only moderately as /? ^ fie- 
Presumably, this increase is due to the less rigid shape of the interface as the energy penalty 
drops. 

Our results for the Debye screening length are shown in fig. 4.24 and fig. 4.25. Rather 
surprisingly, the best fit is for a linear dependence on L^, with a constant term to which 
it tends at very high temperatures when the interface is rigid. The length also appears to 
diverge somewhat as /3 ^ /?c, as expected. The best fit of this divergence is to a curve of 
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Fig. 4.24: Measurements of the Debye screening length for each value of (/3, L^), 
showing a roughly linear dependence on L^ ■ 
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Fig. 4.25: Measurements of the Debye screening length for each value of (/3, Lx), 
showing a mild divergence as /3 ^ /3c, which we fit to a curve of the form a + 
61n(/3 — /3c) for each Lx- 



the form a + feln(/3 — Pc) for each L^, as predicted for completely wet Z{3) interfaces in 
3 + 1 dimensions [44]. 
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Asa curiosity, we can also measure the speed of the random walk of the interface near 
the critical temperature if we turn off the mechanism that constantly re-centres the 
interface. As the temperature drops and more and larger bubbles of phase froth within the 
main phase domains, the interface will combine with, and bud off, more and more of these 
bubbles. Therefore, we would expect it to wander more and more rapidly along the lattice, 
moving in bursts as it merges with bubbles in "front" of it or buds off bubbles "behind" it. 
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Fig. 4.26: This shows the random wandering of the interface, for temperatures 
near the critical temperature, over 40k sweeps after an initial 400 heat-bath sweeps. 



In fig. 4.26, we do indeed see an increase in speed as the temperature drops. In fact, in 
common with many other parameters, the speed appears to diverge as /3 ^ f3c- Of course, 
it could be argued that this is simply an artefact of the Monte-Carlo heatbath algorithm, 
with no relevance to a "real" interface at all. In the absence of a physical "time" to measure 
and use as a yardstick in the Monte-Carlo procedure, the best way to test this hypothesis 
would seem to be to to examine the behaviour of the interface under another Monte-Carlo 
scheme, with different autocorrelation behaviour as (3 drops; we do not do this here. 

In fig. 4.27, we fit the speed to a function of the form 






W-Pc 



Our results suggest values of 



f3, = 8.14(3), f3nw = -0.71(6). 
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Fig. 4.27: This uses the divergence of the speed of the random walk in fig. 4.26 
to estimate the critical temperature and exponent of divergence. 



This value for (3c is certainly consistent with our previous estimates, and the dynamical 
exponent governing the critical acceleration of the random walk is given by (inw- 



4.15. Interface Equilibration 

T T is interesting to note that studies of discrete models [46] and continuum growth 

equations [47] in condensed matter make predictions for the behaviour of the width 

of an interface. The "width" here is usually defined to be the square-root of our second 

moment, 

W{L^,t) = 



—2 



where t is the time since formation of the interface (initial width zero) . 

Models generally predict that for times much smaller than some critical time, tx{Lx)^ 
determined by the interface size, an exponent /3g will govern the growth of the width: 



VF(t)~t^«, 



t^tx{L^)- 



For times much larger than this, as might be expected, the width is predicted to saturate at 
a value whose relationship to the interface size is governed by ac, the roughness exponent: 



W,,t(L,)~L°«. 
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Fig. 4.28: An example of the growth in size of an interface for (3 = 8.5, Lx = 36, 
and a growth profile averaged over all L^ and (3. In each case, (3g is estimated 
from the slope of the dotted line. 
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For the Z{2>) interface, we have already seen behaviour corresponding to the latter of 
these two, since 

< (/>2 >~ L^-^3<^) ^ OG = 0.56(1). 

It can be seen again in fig. 4.28, an example of the growth of a Z{'i) interface. At high t, 
the size saturates as usual, but at low t, we see growth governed by an exponent, as in the 
condensed matter models. For our interface, the average profile yields 

/?G = 0.61(2). 

This, again, could be dismissed as an artefact of the simulation technique, but the manner of 
formation simulated by the Monte-Carlo technique seems eminently plausible for a "real" 
interface, and it is suggestive that the interface shares many characteristics with more 
familiar interfaces in condensed matter physics. 



Cnapter 5. 

Statistical Discussion: 
Bootstrap lecnniques 



5.1. The Bootstrap 



::K|fel|::- INCE IT is not clear that the Wick-subtracted moments will be normally 
::i|||iiM:: distributed, we need a reliable method to construct confidence intervals for them 
::;:::::::::::::::;:!::;::::; whlch docs uot rely ou our knowiug their dlstrlbutlon at the start. A suitable 
method is given by the "bootstrap" principle [48], and we now discuss some aspects of it as 
applied to our analysis. 

Let X = (xi,X2, ■ ■ ■ ,Xn) be a vector of samples from a normal distribution with 
unknown mean and variance. Further, sample the random variable T from the t-distribution 
with n — 1 degrees of freedom, and define the quartile tn-i,a, shortened in future to t^, by 
Prob(T < ta) = a. If x„ and s^ are the sample mean and variance, then 

{Xn + Snto.975 i Xn — Snto.025) 

is the shortest unbiased exact 95% confidence interval for the mean. To obtain the standard 
deviation (a 68% confidence interval for the normal distribution), we simply change the 
subscripts on t to 0.84 and 0.18. 

How, then, can we construct confidence intervals for an unknown distribution? This is 
where the bootstrap principle comes to our aid, as it enables us, as its name suggests, to 
estimate confidence intervals just from our observed data, with no prior knowledge of the 
distribution. Let us denote the population parameter, in our case corresponding to one of 
the moments, by 9. Given a random sample from the population, x, we should like to be 
able to construct an estimate, @ = 0(x), for 6. The idea of the bootstrap is to learn about 
the relationship between 9 and 0(x) by examining that between 9 = 9{xobs) and ^(x*), 
where x* is a resample with replacement from the observed data, Xots (our results for the 
particular moment). This is known as the "nonparametric bootstrap". 

93 
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In the simplest case, we choose to consider the statistic W = Q — 9^ invoking the 
bootstrap principle to estimate its quantiles by those of W* = Q* — 9. We follow the 
procedure given here: 

• Fix 6 at the value given by ^ = 9{xobs)', 

• Sample with replacement from Xof,s to obtain x*; 

• From this, make an estimate of 9{xobs), denoted 9* = 9{x.*); 

• Repeat the previous steps {m — 1) times to obtain m estimates for ^(xq^^); 

• Sort the 9* into order, denoting the pth ordered estimate by 0f -,; 

• The 100(1 — a)th percentile of the distribution of W* is now consistently estimated by 
9l_^ — 9, where we define 9l_^ = 9J •., choosing m so that p = m{l — a) is an integer. 

It then follows that the approximate 95% confidence interval so obtained is 



(29 - ^Q^975 , 29 - 9*y 



Q2b) 



There are two approximations behind this. The first is that ^* — ^ is a consistent estimate 
of ttJ* , where Prob(VF* < ti;* ) = a, and can be reduced as desired by increasing the number 
of resamplings, m. The second relies on W having the same distribution as W* ^ and 
will be true if the distribution of W is pivotal, i.e. depends on no unknown parameters. 
Unfortunately, a sizeable error may be introduced if this is not the case. 



5.2. The Bootstrap t 

'"po overcome this problem, we need to make the statistic W more accurately pivotal. 
Such modification is suggested by considering a sample x = (xi, . . . , Xn) from normal 
distribution with mean 9 and variance o"^. The distribution of x — 0, our previous choice 
of VF, depends on cj^, and is hence non-pivotal. However, introducing the sample variance, 
s^ = -^—^ S"=i(2^i~^)^) we can see that the distribution of y/n{x — 9)/s is pivotal, following 
a t distribution with (n — 1) degrees of freedom. Thus, we "studentise" our previous W: we 
divide by an estimate of the standard error of 0, giving a new statistic: 

nn 0-^ 

T = — ^. 

.(9) 

Using T, we proceed exactly as for W , invoking the bootstrap principle to estimate the 
quantiles by those of 



s(e*) 

and finding the 95% confidence interval to be 

{9-s{9)il,,, , 9-s{9)il 



Q2h) 
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The approximations are as before, but the error associated with the pivotal quahty of the 
statistic should be much reduced. It has been shown, in fact [49], that the coverage error 
of the interval reduces from 0{n~^-''^) to 0{n~'^) for many cases. If T is truly pivotal, the 
error will reduce to zero. 



5.3. The Non-Parametric Delta Method 

TpOR our purposes, there is an added complication: whilst we can calculate the sample 
variance for the second moment, cfp', in the usual manner, the form of the Wick 



subtraction in (t)\Yg and (pws makes it impossible to use a similar formula for them. 
Fortunately, a variance approximation exists for many general statistics T which can take 
the representation 9 = t{F). In this case, F is the cumulative distribution function F (the 
integral of the probability density function) for a distribution X. If the sample data is 
denoted by the usual vector x, then a simple example would be the mean of X, given by 
t{F) = JxdF{x). 

The method works by applying aa Taylor expansion to t(F): 



t{G) = t{F)+ I Lt{x-F)dG{x). 



Lt gives the first derivative of i at F in direction G, and is sometimes called the "influence 
function" : 

L,{.;F)=lim ^ - ^)F + .H.) - t(F) _ m((l - e)F ^ .H, 



e^o e de 



e=0 



Here, Hx{y) is the Heaviside step function, stepping from to 1 at y = x. The empirical 

approximations given by the data values, li = Lt{xi,F), are the "empirical influence values" . 

The "nonparametric delta method" relies on applying this approximation with G = F, so 

that 

1 " 

t{F)=t{F) + -Y,Lt{x,). (5.1) 



n 



Now, the difference, {T — 6), between t{F) and t{F) will follow a normal distribution, by 
the central limit theorem, with mean zero and variance, 

MF) = - I Ll{x)dF{x). 
n J 

We approximate i'l{F) by the sample version, 

1 " 
^^iP) = -2Y.^l (5-2) 
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the "nonparametric delta method variance estimate" . Note that subsitution of F for F in 
(5.1) imphes that ^^ li = 0. For our purposes, we also note that we can use the chain rule 
for statistics which are functions of others, t{F) = T{ti{f), . . . ,tm{F)}, yielding 

Ltix) = ^m—Lt^{x). 

i=l * 

Let us consider the three even moments, denoting the population means of the 
unsubtracted moments by H2, /i4 and hq respectively. First, the second moment: 

t=^^ Ltix) = a,[(l - e)fi2 + ex%=o 

^ l{xi) = {x1 — x^). 
Similarly, the unsubtracted fourth and sixth moments are: 

t = ^^/(xi) = {x1-x^) 

t=^^ l[Xi) = {xf - X6). 

Using the chain rule then allows us to calculate the empirical influence values for the 
subtracted moments: 

t = ^ - 3x^^ 

=^ l{xi) = {xf — x^) — 6x'^{xf — x^) 
t = x^ — 15x^ — 15x^(x'* — 3x^ ) 

^ l{xi) = (x^ - x6) + 90x2^ (xf - x2) - 15x2 (xf - ^) - 15^(xf - ^). 

In each case, the variance is estimated by (5.2); also, as expected, YliK^i) = 0. 

Finally, then, we are able to construct confidence intervals for the subtracted moments, 
and this is how we obtain the standard errors for the fluctuation moments given in tables 
and graphs. For a full data vector, x, of 100,000 values, approximately 2,000 resamplings, 
X*, are necessary. 



5.4. Variance Stabilisation 

A check on the standard errors obtained in this way is provided by using a numerical 
approach rather than (5.2). To do this, we perform a second-level bootstrap analysis 
on each resampled vector, x*, resampling each vector itself (<C 2,000 times) to obtain 
X**. The sample variance of the second- level estimates of the mean of the subtracted 
moment gives the standard deviation for the "parent" first-level estimate. To check for any 
correlation between the mean and deviation of the first-level estimates, which might skew 
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Fig. 5.1: Variance stabilisation in numerical second- level bootstrap procedures, 
illustrated for (/? — 8.75, Lx — 36). The left-hand plots show the second-level 
standard deviations cr(6), against the first-level estimates for the average moment, 
0, together with a fit to a smooth function. The right-hand plots show g{0), the 
integral of the function a(9)~^. In this case, 200 second-level bootstraps are used. 
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the confidence intervals, we then plot one against the other. If a trend is seen, we introduce 
"variance stabilisation" [50] , as follows. A smooth curve is fitted to the data, denoted by 
cr{9). We then construct a compensating function by numerically integrating the inverse: 

as illustrated in fig. 5.1 for data from a (/3 = 8.75, L^ = 36) run. We rescale g{6) so that its 
mean and variance are the same as that of the original data vector x, order it, and plot it 
against the estimates 9. Then, we simply read off the values for 9 which correspond to the 
desired percentiles in g{9), giving the desired confidence interval (quoted for 95%): 



{29 - 9Q^gj^ , 29 - 9*y 



025 J 



The best way to compare the validity of the different methods is to construct "Q-Q 
plots" , where the ordered bootstrap vectors are plotted against quantiles of the standard 
normal distribution, as illustrated in fig. 5.2. The dotted lines on these graphs show what 
would be obtained by a normal distribution with mean and standard deviation given by the 
relevant bootstrap vector. Thus, the most reliable data is that which most closely matches 
the dotted line over the desired confidence interval. The upper three and middle three sets 
of plots show respectively that the bootstrap and bootstrap t distributions are close to, 
but not quite, normal. The variance stabilisation technique is designed to transform to a 
scale where the variance is constant relative to the mean, and will therefore compensate 
for any divergence from normality in the bootstrap t distribution; the lower three plots 
are included merely to show that it does indeed produce curves which compensate for the 
deviations from normality of the bootstrap t curves above. 

As just mentioned, the intervals from variance stabilisation are the most reliable, being 
produced from accurately normal bootstrap t distributions, but they differ little from the 
bootstrap t estimates for our data. As an example of this, the following table lists the 
95% and 1-a confidence intervals for /3 = 8.75 and L^ = 36, matching the data of fig. 5.2. 
Since the variance stabilisation technique requires computationally intensive second-level 
bootstrap analysis, we therefore use straight bootstrap t intervals for our results. 
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Fig. 5.2: Q-Q plots comparing the different methods for obtaining confidence 
intervals, using data for {/3 — 8.75, Lx — 36). Horizontally, the plots refer to the 
second, fourth and sixth moments respectively. Vertically, the upper plots are for 
the basic bootstrap, the central ones for the bootstrap t, and the lower ones for 
variance stabilisation. The dotted lines are those which would be obtained for 
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Table 5.1: Estimates of Confidence Intervals 


Moment 


Technique 


95% (Lower) 


a (Lower) 


a (Upper) 


95% (Upper) 


2nd 


Bootstrap 


6.269 


6.313 


6.396 


6.436 


2nd 


Bootstrap t 


6.276 


6.314 


6.397 


6.442 


2nd 


V. Stabilisation 


6.270 


6.316 


6.395 


6.437 


4th 


Bootstrap 


28.43 


29.38 


31.30 


32.32 


4th 


Bootstrap t 


28.47 


29.43 


31.34 


32.44 


4th 


V. Stabilisation 


28.63 


29.46 


31.44 


32.38 


6th 


Bootstrap 


-5.3 


42.1 


155.2 


214.5 


6th 


Bootstrap t 


12.4 


48.8 


171.6 


257.8 


6th 


V. Stabilisation 


15.7 


51.1 


167.1 


251.3 
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Fig. 5.3: Correlations are given between interface contour shapes separated by 
a certain number of sweeps, for Lx — 36. A particular correlation is shown on 
the right, whilst approximately 4,000 such graphs are averaged on the left, for 
four values of /3. It is clear that the correlation function drops to one half within 
approximately ten sweeps, as illustrated on the left. 



5.5. Correlations and Data Bins 



TN our simulations, we record the full 100k measurements for each raw, even moment of the 
central (50%) contour, but bin all data from other contours, and all smoothed data, in 
groups of one thousand. Thus, we must consider how to ensure that the standard errors are 
not distorted by the binning used. In fig. 5.3, the correlation between two interface contour 
shapes is given as a function of the number of sweeps separating them, and it is seen to drop 
to less than one half after about ten sweeps. Therefore, bins containing far fewer than ten 
samples may give results tainted by unacceptably large correlations, leading to artificially 
small errors, whereas those containing far more than ten will lose some information, giving 
errors which are too large. Thus, we bin our raw data for the central contour into bins of 
size ten before calculating bootstrap errors. 

Illustrated in fig. 5.4 are errors from the non-studentised bootstrap analysis, and from 
the bootstrap t, with more differences apparent between them for the higher moments. The 
graphs show the effect of rebinning the raw moment data into bins of different size. We take 
the full 100k measurements for each moment, and bin the results into groups of between one 
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Fig. 5.4: The bootstrap errors, both studentised ("bootstrap f ) and non- 
studentised, are shown here for (/? — 8.75, L^ — 36). On the left, the ztcr bounds 
(68% confidence intervals) are shown for different binnings of the data; on the 
right, a logarithmic plot of interval sizes against number of bins shows a linear 
dependence for large numbers of bins, whose gradient gives the corresponding 
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and ten thousand. Whereas a random variable with normal distribution would have obey 
a relationship of the form a ~ (number of bins)"*''^, the moment errors flatten off for bins 
containing more than a few hundred measurements. For bins with fewer items than this, 
an estimated power law dependence of a on the number of bins is shown, and this enables 
us to trivially adjust the errors obtained from our bins of size one thousand to those that 
would be obtained from bins of size ten. 



5.6. Statistical Problems Near Criticality 

/^UR results for (3 > 8.35 show reasonable statistical distributions for the moments over the 
100k sweeps. However, those for (3 = 8.25 and small L^ occasionally contain "spikes", 
observations for the moments which are far larger than we would expect, sometimes by 
several magnitudes. This can be seen by comparing the distributions in the following two 
sets of figures: the first, where L^ = 54 and no spikes are apparent; and the second, where 
L^ = 18 and spikes have given each distributions an artificially long, flat upper tail. The 
spikes are unlikely to be caused by wild fluctuations passing through the twist, since similar 
spikes are not seen for wider lattices where fluctuations are larger. It is most likely that 
they are an artefact of the greater energetic instability of the smaller interfaces at very low 
temperatures: if parts of the interface collapse during interaction with nearby bubbles of 
phase, the interface contour may briefly undergo a drastic distortion or partial collapse, 
even appearing to pass through the twist because of the change in Polyakov values there. 

It is difficult to prevent the spikes occurring, but they have a massively disproportionate 
effect on the overall averages, especially for the higher moments. Thus, we would like to be 
able to examine the data with the spikes removed. To do this, we impose an upper cut-off 
on the statistical distributions, as illustrated by the dashed lines in fig. 5.5 and fig. 5.6. 
To ensure that only the spikes are removed, we need to define the cut-off in a statistically 
sensible manner. We take it, for each distribution, to be the median plus 2.5 times the 
interquartile range of the distribution of logarithmic values of data, where the interquartile 
range is defined to run between the 25th and 75th percentiles of the distribution. This 
ensures that the cut-off corresponds only to the upper tip (to be precise, the top 0.04%) of 
a normal distribution overlaid on the right-hand plots of fig. 5.5 and fig. 5.6, ensuring that 
only data sets containing spikes are affected by the procedure. The exact multiplier "2.5" 
can, of course, be varied, but we find that a lower number {e.g. 2.0, cutting off 0.3% of 
the distribution), has too great an effect on the normal contribution to the distributions, 
whereas a higher number {e.g. 3.0, cutting off only 0.002% of the normal distribution) 
permits too much influence from the spikes. 
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Fig. 5.5: Illustrated on the left are the statistical distributions of unsubtracted 
moments for (/3 = 8.25, L^ — 54). The distributions have long upper tails, as 
expected for moments of a roughly normal distribution. To the right are the 
corresponding distributions of the logarithm of the moments, each resembling a 
normal distribution. The dashed vertical lines represent a suggested upper cut-off, 
as discussed in the text. 
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Fig. 5.6: The partners to fig. 5.5, with (/3 — 8.25, L^ = 54), show unexpectedly 
long, flat upper tails in the right-hand plots. The effects of these tails can be 
eliminated by the use of a cut-off defined in the same way as that in fig. 5.5. 
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5.7. Revised Moment Averages (Without Spikes) 

T N this section, we give revised results for moment averages after the removal of spikes as 
discussed, and then some plots for comparison with those of the previous chapter. We 
quote revised results only for the middle (50%) contour and raw Polyakov data: 





Table 5.2: Estimates of < 0^ > from Raw Data (Corrected), Middle (50%) Contour 




(i'--L^ 


18 


24 


30 


36 


42 


48 


54 


Slope 


8.25 


r 44+0.05 


8.09iro'6 


11.9;°;1 


14.3;°;1 


17.8;°;1 


20.8+°:2 


23.11°:^ 


1.33+0.02 


8.35 


A 52+0-03 
^•^^-0.03 


6.69+ro^ 


9.13tH^ 


11.4t°:l 


13.4t°;l 


16.5+°;1 


i8.2+°:i 


1.27+0.01 


8.50 


3.93lH'3 


5.19+ro'3 


7 17+0.05 
'•-•-'-0.05 


8-89tH^ 


10.0t°;l 


11.9l°;l 


i3.8+°:i 


1.15+0.01 


8.75 


3.06+H^ 


4 9-1+0.03 
^•^-■--0.03 


c 09+0.04 
"-••"^^-0.03 


"•<34-0.04 


7.62tH^5 


8.85+ro^6 


q 57+0.06 

^•^'-0.06 


1.05+0.01 


9.00 


2 fi2+°-°2 

^••3^-0.02 


o 44+0.02 


A 00+0. 03 
^•"-•"-0.03 


c 97+0.03 
"-"•^'-0.03 


fi 29+°-°'^ 

"•^^-0.04 


7 07+0.05 

'•"^'-0.05 


7 00+0.05 


1.03+0.01 


Slope 


-0.31+0.03 


-0.36+0.03 


-0.42+0.03 


-0.42+0.03 


-0.44+0.02 


-0.45+0.03 
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Table 5.3: Estimates of < (t)"^ >Ws from Raw Data (Corrected), Middle (50%) Contour 




13' --L^ 
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24 


30 


36 


42 


48 


5^ 


Slope 


8.25 


Q2,.2tll 


llhtl 


273tlt 


326tl 


425+11 


580tlt 


Q58tll 


2.12+0.07 


8.35 


30.6ti;? 


Qi.rtll 


1081^ 


166tl 


1971^ 


322tl 


299t^o 


2.17+0.06 


8.50 


19.2t°j 


27.8l°J 


59.0t?;° 


82.0^2^ 


87.8+2^ 


126+1 


163t^ 


1.97+0.05 


8.15 


9.68+^3^ 


i6.3+r6 


22.3t°;^ 


30.3tl;J 


43.0+lg 


48.2tl;9 


49.7tlj 


1.55+0.04 


9.00 


5.84+°:!^ 


8.76+°:i 


13.0t°;l 


15.6t°j 


90 0+0.8 


28.6tl;J 


28.5tl;i^ 
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-1.25+0.11 
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Table 5.4: Estimates of < (f)^ >Ws from Raw Data (Corrected), Middle (50%) Contour 




P'--L^ 
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24 


30 


36 


42 
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54 
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8.25 


3530 
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-400 


6330 
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23500 
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14500 
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15200 
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-2900 


1.46+0.22 


8.35 


712 
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-55 


1700 

+210 
-190 


2860 

+440 
-390 


5180 

+670 
-580 


4790 

+ 1050 
-840 


10200 

+2100 

-1700 


2640 

+2300 
-1680 


2.38+0.29 


8.50 


259ti 


3501?^ 


1090+lfo 


1710+1° 


91211^0 


1580t^^° 


W20ttl'o 


1.84+0.22 


8.75 


on c+9.9 


1421^° 


1021^0 


61.51^1° 


24411^^ 


-17 0+101.8 

-■-'•"^-78.0 


-634tl\' 


?±? 


9.00 


24.711? 


34.41^-1 


97 0+15.8 

^'•"J-13.9 


05 7+19.5 
•-*^-'-16.2 


f-r. 9+34.2 


914+47 
^J^*-39 


-139+1^ 


0.49+0.11 


Slope 


-1.98+0.11 


-2.10+0.15 


-2.58+0.25 


-1.94+0.60 


-2.09+0.12 
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-1.87+0.07 
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Fig. 5.7: Graphs of even moments against transverse size, Lx, for /3 = 8.5. Data 
are shown for corrected raw Polyakov configurations. 
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It can be seen by comparison with tables 4.6-8 that revision of estimates beyond 
the original errors generally occurs only for low (3 and L^. However, this can have an 
appreciable effect on the calculated power-law dependence of moments on L^ and {(3 — (3c)- 
For comparison, figures corresponding to fig. 4.19, fig. 4.20, fig. 4.21 and fig. 4.22 conclude 
this chapter, showing the effects of removing data spikes. 
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Fig. 5.8: Log plots corresponding to the Wick-subtracted data of fig. 5.7, with 
best fits to a power-law dependence on the lattice width. The previous fits from 
fig. 5.7 are given for ease of comparison. 



For the data shown in the graphs of this chapter {Lr^ = 36, middle contour, raw data), 
the fitted slopes for the revised second, fourth and sixth moments are: 



<02>^(^_^^)-O.42(3)^ <^4>^^^(^_^^)-1.25(ll)^ < / >H..~ (/3 - /^e)" 



-1.94(60) 
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Fig. 5.9: Graphs of even moments against /3, for Lx = 36. Data are shown for 
corrected raw Polyakov configurations. 



Chapter 5: Statistical Discussion: Bootstrap Techniques 
Lj.- Averaged 2nd Moment 



110 



1.20 










' ^■■. 


Slope=-0.42± 0.03 


: 


1.10 


- 


(Old: -0.43+ 0.03) 


- 






'■-.^ 


; 


A- 1.00 


r 




-_ 


5 0.90 


- 


'■'.^ 


'-_ 


0.80 


r 


^, 


-_ 


0.70 






^ ; 



-1.2 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 

L^-Averaged 4th Moment (Wick-Subtracted) 
Slope=- 1.25+ 0.11 



2.5 



"v 2.0 



1.5 



1.0 -, 



(Old: -1.47+ 0.08) 



-1.2 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 



Middle (50%) Contour, Width=36 



L^-Averaged 6th Moment (Wick-Subtracted) 

c p ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' 

Slope=- 1.94+ 0.60 

^.^ (Old: -3.33+ 0.37) 



V 3 



1-1 



<> 



-1.2 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 



Fig. 5.10: Log plots corresponding to the Wick-subtracted data of fig. 5.9, with 
best fits to a power-law dependence on the lattice width. The previous fits from 
fig. 5.9 are given for ease of comparison. 



Referring back to our formulae for the moments, we see that the first two imply the following 
for our unknown functions of (3: 



^^{P-P,f-^m^ A~(/3-/3er^3(i6), 



These two together predict, from our previous formulae, that 



^ ^6 ^ (o o \-2. 08(24) 

< 9 >VKs~ [P - Pc) ^ ', 



extremely close to the observed behaviour. The other sets of revised data now give very 
consistent results for 7 and A. 




6.1. Conclusions 



I 



N THIS thesis, we have examined the properties of the Z{3) interface in two very 
different regimes: 



At high temperatures, we have measured the interface tension on the lattice, and 
extrapolated our results to the infinite-temperature limit, where we have found good 
agreement with the prediction of an instanton calculation. This reinforces the argument 
for the existence of Z(3)-breaking phases in the Euclidean functional integral, although 
their physical existence in the early universe remains in question, and suggests that the 
predictions of perturbation theory are indeed reliable in this regime. We have also examined 
the Debye profile of the interface, but have not been able to compare our results accurately 
with predictions at the temperatures concerned. 

Just above the critical temperature, we have seen that it is possible to follow the 
fluctuations in profile of the interface even though it appears, at first sight, to disappear in 
the phase turbulence. We have observed that the interface maintains a narrow width right 
down to the critical temperature, but undergoes more and more severe fluctuations in shape. 
We have characterised these fluctuations using moments, and have seen that these diverge 
at the critical temperature; and we have constructed a toy scalar field theory model which 
reproduces the observed behaviour to a reasonable approximation. We have also seen that 
whilst the intrinsic width of the interface increases rather slowly near the critical point, the 
Debye screening width appears to diverge logarithmically, in line with the predictions for a 
completely wet interface; but we have not seen direct evidence for this wetting. Finally, we 
have noted that the formation and movement of the interface on the lattice shares many of 
the characteristics of interface behaviour in condensed matter systems. 



Ill 



lUCUNBI ACTI LABOMES 

Cicero 



pltlh! 
(Klingon saying) 





Table A.l: Values of F(2r, 3) for r < 30 


r 


F(2r,3) 


1 


2.9257441147538591e-02 


2 


2.8588133947184931e-03 


3 


7.7625872330950078e-04 


4 


3.1782999811761687e-04 


5 


1.5992165454570213e-04 


6 


9.1150725927608572e-05 


7 


5.6469174319874464e-05 


8 


3.7115784724567322e-05 


9 


2.5484489828066094e-05 


10 


1.8087564731975696e-05 


11 


1.3170447587156097e-05 


12 


9.7840479099121027e-06 


13 


7.3840421441212328e-06 


14 


5.6428031018154392e-06 


15 


4.3549604790025167e-06 


16 


3.3872387617351618e-06 


17 


2.6505179822474595e-06 


18 


2.0836256901432198e-06 


19 


1.6435963519300661e-06 


20 


1.2996333087139130e-06 


21 


1.0292617591260539e-06 


22 


8.1581683728676361e-07 


23 


6.4676674563294586e-07 


24 


5.1257056950486260e-07 


25 


4.0588577520678718e-07 


26 


3.2100885546335727e-07 


27 


2.5347419875900234e-07 


28 


1.9976210969186481e-07 


29 


1.5708328861159867e-07 


30 


1.2321764912912764e-07 


1 
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Table B.l: Estimates of < (j)^ > from Raw Binned Data, Lower (25%) Contour 




P'-'L:, 


18 


24 


30 


36 


42 


^5 


54 


Slope 


8.25 


5 ftQ+0.08 


7 5n+o-04 


io.9+°:l 


i3.5+°:l 


16.6t°;l 


19.7t°i 


01 7+0.2 


1.30+0.02 


8.35 


4^1+0.11 
*-OJ--o.oi 


"•"-•^-0.04 


8.67t°;°« 


10.6+°:1 


12.4t°;l 


15.5t°;l 


17.2t°;l 


1.22+0.01 


8.50 


q 09 + 0.02 

"3-O^-0.03 


5 nn+0.03 
^•'-"-'-0.03 


6.85t°;°^, 


8.63tH^ 


9.57tH^ 


ii.5tS:J 


i3.itS;l 


1.13+0.01 


8.75 


04+0.02 

''•'-'^-0.02 


4 11+0.03 
^•-^^-0.03 


c 91+0.04 
"-"•^^-0.03 


"•^^-0.03 


7 4fi+0.05 
'•*'3-0.04 


cq+0.07 
O-O'3-0.04 


q Q4+0.07 

^•"3^-0. 06 


1.04+0.01 


9.00 


9f.i+0.01 
^•"J--0.02 


40+0.02 
"-'•^"-'-0.02 


4 Q4+0.03 
^•"3^-0. 02 


c 91+0.03 
"-"•^^-0.03 


fi iq+0.05 

"•-'^^-0.04 


7 97+0.05 
'•^'-0.05 


7 77+0.04 
'•''-0.05 


1.02+0.01 


Slope 


-0.33+0.02 


-0.32+0.02 


-0.39+0.02 


-0.40+0.03 


-0.41+0.01 


-0.42+0.02 


-0.44+0.03 




1 





Table B.2: Estimates of < 4>^ > from Raw Binned Data, Middle (50%) Contour 




13'--L^ 


18 


24 


30 


36 


42 


48 


54 


Slope 


8.25 


(^ 07+0.22 


97+0.12 
°-^'-0.03 


19 0+0.2 
1Z.CS_Q2 


14.4+°:1 


17.8t°;l 


20.8t°i 


23.lt°;? 


1.20+0.03 


8.35 


c- P9+0.20 
"-'•'-'^-0.01 


fi 71+0-05 

O-'-L-0.05 


q 14+0.09 

^•-•-^-0.05 


11.41°:! 


13.4;°;1 


16.5;°;1 


18.2;°;1 


1.21+0.01 


8.50 


q Qq + 0.03 

o.yo_Q Q3 


5 iq+0.03 

"3--'-y-0.04 


7 17+0.05 
'•-■-'-0.05 


8.89;n^ 


10.0;°;1 


11.9;°;! 


13.8;°;1 


1.15+0.01 


8.75 


07+0.02 

"3-'-" -0.02 


A 91+0.03 

*-^J--0.03 


c q9+0.04 
"-'•"3^-0.03 


fi 04+0.05 

'J-<J'±-0.04 


7 co+0.06 
'•"'3-0.05 


oc:+0.07 
O.»0_QQ5 


9.57tH^ 


1.05+0.01 


9.00 


9 r.^+om 

z:.u^_Q Q2 


44+0.02 
"-••^^-0.02 


A qO + O.OS 
4:.00_Q Q2 


c 97+0.03 
"-"•^'-0.03 


fi 29+°-°^ 

"•^^-0.03 


7 q7+0.05 
'•"3' -0.05 


7 OQ+0.04 

'•oy_QQg 


1.03+0.01 


Slope 


-0.41+0.01 


-0.37+0.02 


-0.44+0.02 


-0.43+0.03 


-0.44+0.02 


-0.45+0.03 


-0.46+0.04 




1 





Table B.3: Estimates oi < (f)^ > from Raw Binned Data, Higher (75%) Contour 




I3'--L^ 


18 


24 


30 


36 


42 


48 


54 


Slope 


8.25 


q 91+0.50 
^•^-■--0.30 


q ftq+0.48 

y-O'3-0.08 


13.41°:^ 


14.6l°:l 


17.6t°;l 


20.lt°;l 


22.0t°;? 


0.91+0.03 


8.35 


fi iq+0.36 

"•-'^^-0.04 


7 qi+0.04 

' -"3-^-0.05 


9.7lt°;°^ 


11.9+°:1 


13.9t°:l 


16.5t°:l 


18.0+°:1 


1.08+0.02 


8.50 


4 44+0.02 

*-*^-0.02 


5 qf;+0.03 
"3-^0-0.04 


7.86tH5 


9.48tH^ 


io.6t°:l 


12.3t°:l 


i4.i+°:i 


1.04+0.01 


8.75 


cr+om 

'3-'JJ-0.02 


A SS+0.03 
4-OO-0.03 


5 qfi+0.04 


6.94;n^ 


s.iTint 


9.29;°;°^ 


9.981^^ 


0.92+0.01 


9.00 


q 17+0.01 
"3--L'_0.02 


4 nzL+0-02 


4 04+0.03 
4-^4-0.02 


5 7q+o.03 

"-••'^-0.03 


f- 77+0.04 

O-' '-0.03 


7 7c:+0.05 
'•'^-0.05 


qp, + 0.04 

o-"3U_o.05 


0.89+0.01 


Slope 


-0.41 + 0.03 


-0.37+0.02 


-0.42+0.02 


-0.39+0.03 


-0.40+0.02 


-0.41+0.03 


-0.42+0.04 




1 
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Table B.4: Estimates of < 0'* >ws from Raw Binned Data, Lower (25%) Contour 




I3'--L^ 


18 


24 


30 


36 


42 


^5 


54 


Slope 


8.25 


1150t?^° 


or;i+81 


uitlt 


331+^^ 


414lfe 


568+11 


642lfo 


0.27+0.25 


8.35 


438+«f 


55.51M 


102+^ 


144+^ 


166+^ 


313+1^ 


286+? 


1.84±0.47 


8.50 


18.3+°:^ 


26.7+r8 


54.6+?;° 


79.5+fj 


85.0t|° 


121+1 


160l« 


1.99+0.05 


8.75 


9.48;°;1 


16.lir4 


22.0+^4 


m.?.tll 


41.211:6 


49.2;2;3 


49.2;iJ 


1.58+0.03 


9.00 


6.29tS;l^ 


a QQ+0.25 
O.»»_o.i9 


13.4+°i 


17 r, + 0.6 


23.8t°;^ 


29.2tl;^ 


28.7tJJ 


1.51±0.03 


Slope 


-2.00±0.43 


-1.25±0.10 


-1.40±0.06 


-1.23±0.08 


-1.19±0.04 


-1.24±0.10 


-1.29±0.11 




1 





Table B.5: Estimates of < (j)'^ >Ws from Raw Binned Data, Middle (50%) Contour 




I3'--L,, 


18 


24 


30 


36 


42 


48 


5^ 


Slope 


8.25 


2270t^«° 


334^125 


mtlf 


463+^f 


492t^| 


585_i3 


690+?o 


-0.54±0.27 


8.35 


775t?f 


63.1+?:^ 


111 + 2 


166+^ 


199t^ 


344t?l 


299t" 


1.31+0.49 


8.50 


19.51°;^ 


28.3_o;8 


59.7;?;! 


00 r+3.3 

O'^-t'-i.e 


SS.Olle 


l2Qtl 


16311 


1.95+0.05 


8.15 


9.92+°:f7 


16.7+°:^ 


22.3t°J 


30.3tJ;27 


43.2tlj 


48.6+1:3 


49.5+?-i 


1.54±0.04 


9.00 


r oc+0.14 

t).»t>_o 17 


8.80+°:i 


10 1+0.5 

^--'•-'--0.3 


15.6tn 


23.3t°;^ 


28.6+J:2 


28.611:^ 


1.54±0.03 


Slope 


-2.40+0.44 


-1.33+0.12 


-1.53+0.17 


-1.47+0.08 


-1.29+0.05 


-1.26+0.12 


-1.32±0.11 




1 





Table B.6: Estimates of < cj)'^ >Ws from Raw Binned Data, Higher (75%) Contour 




13' --L^ 


18 


24 


30 


36 


42 


48 


54 


Slope 


8.25 


3820_3oo 


000 + 246 
oOo_22 


1370+fo°o 


483+?? 


5l6tll 


705tl? 


9081^^ 


-1.01+0.24 


8.35 


1750t^r 


89.2+?:^ 


139+^ 


210+^ 


268t^ 


375^1^ 


3711^1 


0.88+0.54 


8.50 


30.8+°;^ 


53.4+?:^ 


80.7tl;^ 


114+2 


117t^ 


153t? 


191^4 


1.62±0.04 


8.15 


18.5l°;l 


30.11°:^ 


38.7;[!j 


48.9;^:? 


64.41?:° 


71.811:? 


76.81?:? 


1.32±0.03 


9.00 


12.3l°j 


18.0+°:! 


23.8;°;I 


27.81°;^ 


39.0;^;° 


45.2lJ:^ 


47.ll?:? 


1.28+0.02 


Slope 


-2.33+0.52 


-1.20+0.25 


-1.33+0.24 


-1.21±0.07 


-1.11±0.05 


-1.18+0.07 


-1.26+0.06 




1 
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Table B.7: Estimates of < 0^ >ws from Raw Binned Data, Lower (25%) Contour 




/3'--^. 


18 


24 


30 


36 


42 


48 


54 


Slope 


8.25 


3450000 

+730000 
-310000 


540000 

+699000 
+ 105000 


654000 

+ 166000 
-91000 


84100 

+94600 
-3900 


49200 

+5100 
-2900 


34500 

+7200 
-3500 


43000 

+7900 
-2900 


-4.38+0.24 


8.35 


1370000 

+7100000 
-20000 


1580 

+420 
-160 


3400 

+ 1050 
-300 


3760 

+600 
-470 


3210 

+660 
-640 


19200 

+6700 
-900 


1700 

+ 1580 
-880 


0.04+1.18 


8.50 


278lg 


OOO QQ 


1090l?l° 


2550tgr 


971+foI 


-|oin+1590 
101U_,_47Q 


71 c:+539 


1.60+0.29 


8.15 


ib.2t]y 


135;t« 


7o.9;?°:i 


lOllif 


19llf, 


07 c: + 156.6 
■-"•^-60.7 


-mitll' 


0.73+0.17 


9.00 


39.6+1^ 


34.2t^j 


33.4+^^1^ 


19.61?^:? 


64 2+25-5 


-157+^? 


-149+^9 


-0.52+0.11 


Slope 


-4.99+1.05 


-3.49+0.57 


-4.08+0.37 


-3.11+0.38 


-2.75+0.12 


-1.98+0.47 


-3.00+0.32 




1 





Table B.8: Estimates of < (j)^ >ws from Raw Binned Data, Middle (50%) Contour 




/3'--L, 


18 


24 


30 


36 


42 


48 


54 


Slope 


8.25 


5720000 

+ 1370000 
-510000 


391000 

+441000 
+45000 


1280000 

+720000 
-180000 


297000 

+320000 
-16000 


199000 

+295000 
-23000 


24900 

+ 13600 
-3800 


71700 

+112000 

-7300 


-4.40+0.39 


8.35 


1880000 

+970000 
-20000 


1830 

+260 
-170 


3710 

+640 
-330 


5150 

+700 
-560 


6180 

+ 1300 
-980 


42800 

+41100 
-2700 


2640 

+2570 
-1400 


-2.87+1.05 


8.50 


289+^^ 


380+^7 


1260t?l° 


2160+?;,° 


939+^1^ 


i&mXf^ 


968t^«^ 


1.82+0.26 


8.15 


109+1^ 


167+^4 


iioti 


98.0+1^7 


285+-^ 


77.5+^?y 


-625+^^^ 


0.44+0.25 


9.00 


21.ltll 


48.2+^:^9 


39.411?:° 


-I9.7tf,l 


-65.2+34;7 


-213t^^ 


-125+^^ 


1.06+0.14 


Slope 


-5.57+0.98 


-3.16+0.55 


-3.95+0.47 


-3.33+0.37 


-3.15+0.23 


-2.27+0.82 


-2.98+0.29 




1 





Table B.9: Estimates of < cj)^ >Ws from Raw Binned Data, Higher (75%) Contour 




/3'--i. 


18 


24 


30 


36 


42 


48 


54 


Slope 


8.25 


7170000 

+ 1170000 
-620000 


653000 

+821000 
-43000 


1580000 

+690000 
-180000 


257000 

+323000 
-12000 


59300 

+5600 
-2800 


73900 

+9500 
-6300 


157000 

+46000 
-9000 


-4.40+0.39 


8.35 


3880000 

+500000 
-40000 


3750 

+390 
-260 


5410 

+390 
-270 


9560 

+ 1130 
-800 


13700 

+ 1700 
-1500 


31000 

+3400 
-100 


16000 

+7100 
-2300 


-3.69+0.99 


8.50 


749tf8 


oqon+26650 


2180t?°° 


3900lf8° 


2650ti° 


3650+^[;0 


2190+^fo 


1.54+0.27 


8.15 


334;f2 


6601^^ 


792+_ll 


935l^f 


lOTOliJO 


10201^;,° 


338llg 


0.94+0.26 


9.00 


17311^ 


315lp 


2>7^tf, 


3411^^ 


b'iQtll 


274;^^ 


2071^^ 


0.60+0.19 


Slope 


-5.36+1.01 


-2.22+0.67 


-2.58+0.70 


-2.39+0.24 


-1.99+0.07 


-2.36+0.23 


-2.85+0.11 




1 
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Table B.IO: Estimates of < > from Smooth Binned Data, Lower (25%) Contour 




I3'--L^ 


18 


24 


30 


5^ 


^^ 


48 


54 


Slope 


8.25 


n 14 + 0.23 


10.1^°:? 


u.^tli 


17.8tr2 


91 5+0-2 


25.6^^2 


28.6+°:^ 


1.21+0.02 


8.35 


6.02+°;- 


7 71 +0.06 

'•'-■--o.oe 


10.6t°;l 


13.0t°;l 


i5.oi°:l 


18.8t°;? 


20.6t°;l 


1.16+0.02 


8.50 


4.38+H! 


tr 7C + 0.04 
^•'^-0.05 


s.ootH^ 


10.lt°;l 


11.0+°;1 


i3.i;S:l 


i5.it^;l 


1.14+0.02 


8.75 


q 1/1+0.03 

"-'•-'-^-0.03 


4 40+0.04 
^•4'J-0.04 


5 fi7+0-05 
"J-'J'-0.04 


6.82;n^ 


8.19+ro^ 


9.5liro^ 


10.3;°;! 


1.09+0.01 


9.00 


9 c/i+0.02 
^•^^-0.02 


4«+0.02 
"^•^"-0.02 


4 cn+0.04 


c 4fi+0.04 


6.54+ro^ 


7 71 +0.06 

'•'-•^-o.oa 


8.26t°;°^ 


1.10+0.01 


Slope 


-0.51+0.03 


-0.46+0.03 


-0.51+0.03 


-0.50+0.03 


-0.50+0.02 


-0.51+0.02 


-0.53+0.03 




1 





Table B.ll: Estimates of < 0^ > from Smooth Binned Data, Middle (50%) Contour 




^'■■L,. 


18 


24 


30 


36 


42 


48 


54 


Slope 


8.25 


7qn+0-i9 

'•^'^-0.12 


10.6+°:^ 


15.4t°J 


18.31°;? 


22.51°:^ 


26.61°:^ 


30.31°:^ 


1.26+0.01 


8.35 


5 fi7+0-05 
"-••"' -0.02 


8.00iro? 


11.0t°;l 


13.6l°;l 


i5.8i°:l 


19.61°:? 


21.51°:? 


1.23+0.01 


8.50 


4 qR+0.04 


5.79iro^ 


8.171°:°? 


10.ll°;l 


11.41°:! 


13 5+0-1 
1J.D_04 


i5.6i°:l 


1.17+0.01 


8.15 


q 07+0. 03 
"^•"-"-0.03 


4 /in+o-03 


5 fif^+0-05 
'-'•""-'-0.04 


6.801°;°^ 


Q 99+0.06 
"•^^-0.06 


50+0-08 
^•"3^-0. 05 


io.3l°:l 


1.11+0.01 


9.00 


9 47+0.01 
^•^'-0.02 


07+0.02 

-^•"^'-0.02 


4 /in+0-04 
4-40_o.o3 


5 40+0-04 

"'•*'-'-0.04 


fi 59+0-05 
"•"J^-0.04 


7 fiq+0-06 

'^•0^-0.06 


n 9C-+0.05 

°-^'J-0.06 


1.13+0.01 


Slope 


-0.50+0.03 


-0.48+0.03 


-0.54+0.03 


-0.52+0.04 


-0.52+0.02 


-0.53+0.03 


-0.56+0.03 




1 





Table B.12: Estimates oi < cj)^ > from Smooth Binned Data, Higher (75%) Contour 




/3 ' ■ • ix 


18 


24 


30 


36 


42 


48 


54 


Slope 


8.25 


r^ qc-+0.15 

"•"J^-o.og 


11 q+0-2 


16.11°:^ 


19.81°:^ 


24.01°:^ 


90 q+0.2 

Zo.O_Q 2 


q9 1+0.4 


1.23+0.01 


8.35 


6.88_o:o3 


9 29+°-°l 


i2.ii°:i 


i5.olS:l 


i7.5l°:l 


21.01°:? 


22.8l°:l 


1.12+0.01 


8.50 


r q9+0.03 
"-'•"-'^-0.03 


7 19+0.04 
'•-■-^-0.05 


9.50l°:°o^6 


11.41°:! 


i2.6i°:l 


i4.5i°:l 


i6.6i°:l 


1.03+0.01 


8.75 


4 07+0.03 

^•'-''-0.03 


c C4+O.O4 
"-'•"-'^-0.04 


6.79l°:°oi 


7 Q9+0.05 
'^•^^-0.04 


9.351°:°^ 


io.6i°:J 


ii.4i°:i 


0.94+0.01 


9.00 


q q/1+0.01 
"^•"J^-0.03 


4 qq + 0.03 

^•"J'J-0.02 


5 40+0.04 

"3-4'-'-0.03 


(^ q7+0.03 

O-3'-0.04 


7 50+0-05 
'^•"3 '-'-0.04 


8.631°:°^ 


9.231°:°^ 


0.95+0.01 


Slope 


-0.40+0.03 


-0.42+0.03 


-0.47+0.03 


-0.48+0.03 


-0.49+0.02 


-0.51+0.02 


-0.54+0.03 




1 
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Table B.13: Estimates of < 0^ >Ws from Smooth Binned Data, Lower (25%) Contour 




13' --L^ 


18 


24 


30 


36 


42 


48 


54 


Slope 


8.25 


1730t-° 


6281^° 


134011^° 


QS1+478 
9^1-65 


936t^^5 


1050t^° 


1130t^60 


-0.18+0.15 


8.35 


729tf 


71.8+1^ 


132+^ 


imtl 


182^^ 


456t^^^ 


29lt^i 


-0.11+0.38 


8.50 


24.1+r8 


31.3+J;^ 


62.lt|^ 


m.2+_f,i 


82.51^-^ 


112+1^ 


148+1 


1.65+0.07 


8.75 


10.4+°-i 


16.11SJ 


20.61^;° 


2Q.?.+_ll 


34.9+1:8 


38.21?:^ 


33.7+?:2 


1.21+0.05 


9.00 


Q.Mtlf, 


7.91+°i^ 


ll-5t^J 


13.5t°i 


18.7ir6 


20.5+1:^ 


I8.O+1:? 


1.19+0.05 


Slope 


-2.69+0.34 


-1.63+0.15 


-1.84+0.14 


-1.67+0.10 


-1.54+0.06 


-1.65+0.12 


-1.76+0.11 




1 





Table B.14: Estimates of < 0'* >ws from Smooth Binned Data, Middle (50%) Contour 




I3'--L,, 


18 


24 


30 


36 


42 


48 


54 


Slope 


8.25 


2230t?^° 


861+lf 


1620l?l° 


7i5;^r 


SlQt}^ 


912li 


umtlf 


-0.56+0.17 


8.35 


497tf 


79.8+^:° 


1321^ 


1871^ 


22'itl 


4221^^ 


30011^ 


-0.02+0.31 


8.50 


24.0+Ji 


31.111:3 


U.Atll 


87.71^:^ 


82.3l|^ 


11811° 


1481^ 


1.68+0.07 


8.15 


10.3+°:^ 


-"-o-i-o.e 


2o.o;j:° 


or; 0+1.3 


35.91?:^ 


36.5lf] 


32.21?:° 


1.17+0.06 


9.00 


5.04+°-° 


7-1-1+0.34 
'•-'--'--0.25 


10.1;°:^ 


J--'--"J-0.5 


16.81°:? 


19.111:^0 


17.311:^ 


1.27+0.05 


Slope 


-2.76+0.25 


-1.80+0.17 


-1.88+0.17 


-1.74+0.10 


-1.61+0.04 


-1.62+0.13 


-1.82+0.11 
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Table B.15: Estimates of < cj)'^ >Ws from Smooth Binned Data, Higher (75%) Contour 




/^■••i. 


18 


24 


30 


36 


42 


48 


54 


Slope 


8.25 


21801?!° 


909lf 


1280lij^° 


967li 


98911^ 


13001^° 


16201^° 


0.08+0.16 


8.35 


10601^° 


138l| 


1921^ 


2841^ 


36611J 


5221?^ 


444+1? 


0.24+0.33 


8.50 


43.0lJ:^ 


68.8l2:o 


1021^ 


1431^ 


1311^ 


1691? 


1961^ 


1.36+0.05 


8.15 


24.31°:^ 


37.71J:°8 


45.91U 


56.01^:^ 


71.51?:^ 


74.51?:° 


74.011^ 


1.06+0.04 


9.00 


15.41°:! 


21.61°:^ 


27.11°:^ 


3o.oi[;:f 


41^41^:^ 


43.61^:^ 


45.011:^ 


1.02+0.03 


Slope 


-2.28+0.34 


-1.44+0.13 


-1.38+0.15 


-1.43+0.05 


-1.34+0.04 


-1.45+0.06 


-1.50+0.02 
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Table B.16: Estimates of < (j)^ >Ws from Smooth Binned Data, Lower (25%) Contour 




/3'--^. 


18 


24 


30 


36 


42 


48 


54 


Slope 


8.25 


3450000 

+330000 
-180000 


540000 

+ 170000 
-54000 


654000 

+260000 
-150000 


84100 

+507000 
-38000 


49200 

+67000 
-39000 


34500 

+ 166000 
-57000 


43000 

+359000 
-30000 


-1.62+0.20 


8.35 


1370000 

+51000 
-6000 


1580 

+630 
-270 


3400 

+9720 
-1340 


3760 

+820 
-640 


3210 

+ 1130 
-1020 


19200 

+606000 
-13000 


1700 

+2940 
-1650 


-5.38+1.17 


8.50 


562+^^ 


804+^^° 


1840lf9° 


5580+^^,^ 


i7mtf,l 


2020+^^^° 


lOlOtfyO 


1.19+0.32 


8.15 


1561?^ 


2^ltlt 


222lf^ 


303l^f 


4851^^ 


338l^f 


n^ 7+116.7 

O'-'-86.0 


1.02+0.09 


9.00 


67.0+^^4° 


97.5tl- 


119+^3 


lUtt 


121+^^ 


mt% 


208+^? 


0.90+0.08 


Slope 


-5.40+0.94 


-3.30+0.47 


-3.96+0.37 


-3.25+0.48 


-3.45+0.33 


-3.47+0.51 


-3.55+0.22 




1 





Table B.17: Estimates of < 0^ >Ws from Smooth Binned Data, Middle (50%) Contour 




13'--L,, 


18 


24 


30 


36 


42 


48 


54 


Slope 


8.25 


5720000 

+660000 
-370000 


391000 

+910000 
-80000 


1280000 

+850000 
-290000 


297000 

+452000 
-22000 


199000 

+406000 
-51000 


24900 

+250000 
-40000 


71700 

+610000 
-39000 


-2.55+0.24 


8.35 


1880000 

+28000 
-8000 


1830 

+560 
-360 


3710 

+940 
-500 


5150 

+960 
-790 


6180 

+2600 
-1900 


42800 

+ 109400 
-5700 


2640 

+5970 
-2450 


-5.29+1.10 


8.50 


546+S 


775+^°^ 


nmtfw 


4000+280° 


1850_29o 


3380+^r 


1830t^J° 


1.66+0.22 


8.15 


207+^^ 


286+^^ 


295+- 


313t^« 


574t^f 


484t^f 


-157+^f 


0.92+0.08 


9.00 


49.7+^;0 


76.41^^/ 


9o.2t?^:^ 


7Q Q+29.8 


53.7_3^;3 


19.3+^°:^ 


1781^^ 


n/a 


Slope 


-5.78+0.92 


-3.35+0.59 


-3.79+0.52 


-3.24+0.32 


-3.39+0.25 


n/a 


-3.71+0.17 
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Table B.18: Estimates of < (p^ >Ws from Smooth Binned Data, Higher (75%) Contour 




/3'--^. 


18 


24 


30 


36 


42 


48 


54 


Slope 


8.25 


7170000 

+750000 
-420000 


653000 

+710000 
-40000 


1580000 

+670000 
-240000 


257000 

+322000 
-26000 


59300 

+21000 
-13000 


73900 

+57000 
-32000 


157000 

+65000 
-20000 


-2.60+0.19 


8.35 


3880000 

+ 110000 
-30000 


3750 

+2700 
-1500 


5410 

+800 
-600 


9560 

+3400 
-2200 


13700 

+4800 
-3500 


31000 

+ 10500 
-2800 


16000 

+5700 
-3200 


-4.14+0.86 


8.50 


14601^^0 


36101^^^° 


5230lifo° 


9QmtfiQ 


70101^^° 


7880;^f° 


57001^°^° 


1.56+0.18 


8.15 


707+f3 


1380+-° 


1860+lfo 


2290t?^o 


2640t?°° 


3080+?^° 


2140^250 


1.21+0.10 


9.00 


4111^6 


7Utf, 


812+f6 


1070t^^ 


1390ti?° 


1060t^J° 


1420+1^0 


1.03+0.07 


Slope 


-5.25+0.99 


-2.74+0.43 


-2.43+0.57 


-2.44+0.22 


-2.29+0.10 


-2.46+0.17 


-2.36+0.19 




1 



Rei( 



erences 



[1] E.g. G. Forgacs, R. Lipowsky and Th.M. Nieuwenhuizen in "Phase Transitions and 
Critical Phenomena", vol. 14, Academic Press Ltd. (London) 1991 (ed. C. Domb and 
J. Lebowitz), and references therein. 

[2] A. Polyakov, Phys. Lett. B72 (1978) 477. 

[3] L. Susskind, Phys. Rev. D20 (1979) 2610. 

[4] H.J. Rothe, "Lattice Gauge Theories" (ch. 19—20), World Scientific (1992). 

[5] V. Dixit and M. C. Ogilivie, Phys. Lett. 269 (1991) 353. 

[6] J. Ignatius, K. Kajantie and R. Rummukainen, Phys. Rev. Lett. 68 (1992) 737. 

[7] Z. Frei and A. Patkos, Phys. Lett. B229 (1989) 102. 

[8] K. Kajantie, L. Karkkainen and K. Rummukainen, Nucl. Phys. B357 (1991) 693. 

[9] R. Brower, S. Huang, J. Potvin, C. Rebbi and J. Ross, Phys. Rev. D46 (1992) 4736; 
B. Grossmann, M. L. Laursen, T. Trappenberg and U.-J. Wiese, Nucl. Phys. B396 
(1993) 584. 

[10] V. M. Belyaev, I. I. Kogan, G. W. Semenoff and N. Weiss, Phys. Lett. B277 (1992) 
331. 

[11] N. Weiss, UBCTP93-23 (liep-ph/9311233). 

[12] I. I. Kogan, Phys. Rev. D49 (1994) 6799. 

[13] A. V. Smilga, Ann. Phys. 234 (1994) 1. 

[14] C. P. Korthals Altes, private communication. 

[15] A. Roberge and N. Weiss, Nucl. Phys. B275 (1986) 734. 

[16] C. Bernard, Phys. Rev. D9 (1974) 3312; for a more recent exposition, see e.g. 
J. I. Kapusta, "Finite- Temperature Field Theory", CUP, Cambridge (1989). 

[17] L. D. McLerran and B. Svetitsky, Phys. Rev. D24 (1981) 450. 

[18] B. Svetitsky and L. G. Yaffe, Nucl. Phys. B210[FS6] (1982) 423. 

[19] B. Svetitsky, Phys. Rep. 132 (1986) 1. 

[20] N. Weiss, Phys. Rev. D24 (1981) 475. 

[21] T. Bhattacharya, A. Gocksch, C.P. Korthals Altes and R. D. Pisarski, Nucl. Phys. B383 
(1992) 497. 

[22] B. DeWitt, Phys. Rev. 162 (1967) 1195. 

[23] L. F. Abbott, Nucl. Phys. B185 (1981) 189. 

123 



References 124 

[24] T. Bhattacharya, A. Gocksch, C. Korthals Altes and R. D. Pisarski, Phys. Rev. Lett. 66 

(1991) 998. 

[25] V. M. Belyaev and V. L. Eletsky, Zeit. Phys. C45 (1990) 355. 

[26] J. B. Kogut, Rev. Mod. Phys. 55 (1983) 775. 

[27] Michael Creutz (ed.), "Quantum Fields on the Computer", World Scientific, Singapore 
(1992). 

[28] C. P. Korthals Altes, A. Michels, M. Stephanov and M. Teper, OUTP-96-10P. 

[29] E. D'Hoker, Nucl. Phys. B200[FS4] (1982) 517. 

[30] E. D'Hoker, Nucl. Phys. B201 (1982) 401. 

[31] Following L. D. Landau and E. M. Lifschitz, "Statistical Physics" (Part I, Ed. 3), §78: 
"Thermodynamic Quantities for a Classical Plasma". 

[32] M. Gross and J. F. Wheater, Nucl. Phys. B240 (1984) 253. 

[33] J. Christensen, G. Thorleifsson, P. H. Damgaard and J. F. Wheater, Phys. Lett. 276B 

(1992) 472. 

[34] R. V. Gavai, F. Karsch and B. Petersson, Nucl. Phys. B322 (1989) 738; M. Fukugita, 
M. Okawa and U. Ukawa, Phys. Rev. Lett. 63 (1989) 1768; Nucl. Phys. B337 (1990) 
181; N.A. Alves, B.A. Berg and S. Sanielevici, Phys. Rev. Lett. 64 (1990) 3107. 

[35] M. Creutz, Phys. Rev. D21 (1980) 2308. 

[36] A. D. Kennedy and B. J. Pendleton, Phys. Lett. 156B (1985) 393. 

[37] N. Metropohs, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller and E. Teller, 
J. Chem. Phys. 21 (1953) 1087. 

[38] N. Cabibbo and E. Marinari, Phys. Lett. 119B (1982) 387. 

[39] S. L. Adler, Phys. Rev. D23 (1981) 2901. 

[40] F. R. Brown and T. J. Woch, Phys. Rev. Lett. 58 (1987) 2394. 

[41] M. Creutz, Phys. Rev. D36 (1987) 515. 

[42] K. Akemi et al, Nucl. Phys. B (Proc. Suppl.) 30 (1993) 253. 

[43] C.P. Korthals Altes, A. Michels, M. Stephanov and M. Teper, Nucl. Phys. B (Proc. 
Suppl.) 42 (1995) 517. 

[44] T. Trappenberg and U.-J. Wiese, Nucl. Phys. B372 (1992) 703. 

[45] J. Christensen, G. Thorleifsson, P. H. Damgaard and J. F. Wheater, Nucl. Phys. B374 
(1992) 225. 

[46] F. Family and T. Vicsek, J. Phys. A18 (1985) L75; F. Family, J. Phys. A19 (1986) 
L441; P. Meakin, P. Ramanlal, L. M. Sander and R. C. Ball, Phys. Rev. A34 (1986) 
5091; J. M. Kim and J. M. Kosterlitz, Phys. Rev. Lett. 62 (1989) 2289. 



References 125 

[47] S. F. Edwards and D. R. Wilkinson, Proc. R. Soc., London, A381 (1982) 17; M. Kardar, 
G. Parisi and Y.-C. Zhang, Phys. Rev. Lett. 56 (1986) 889. 

[48] B. Efron and R. J. Tibshirani, "An Introduction to the Bootstrap", Chapman and Hall, 
London (1993); J. R. Carpenter, "Simulated Confidence Intervals for Parameters in 
Epidemiological Models", D. Phil thesis. University of Oxford (1996). 

[49] P. Hall, "The Bootstrap and Edgeworth Expansion", Springer- Ver lag, London (1992). 

[50] A. C. Davison and D. V. Hinkley, "Bootstrap Methods", CUP, Cambridge (1996). 



